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Development of High-Speed Water-Based 
Aircraft 


ERNEST G. STOUT* 
Consolidated Vultee Aircraft Corporation 


SUMMARY 


Outstanding progress has been made in recent years in the 
field of seaplane hydrodynamics which has done much to bring 
the faltering seaplane up to modern standards of aerodynamic 
performance and design efficiency In fact, extensive research 
programs, utilizing the dynamically similar mode] technique, 
promise to place consideration of water-based aircraft high on the 
program of current defense planning. 

This paper outlines the salient features involved in the develop- 
ment of the dynamic model and associated research techniques 
into a practical, accurate research tool for designing water-based 
aircraft that demand little, if any, compromise with contemporary 
aerodynamic design. New seaplane design criteria are discussed, 
and their application to high-performance water-based aircraft 
isanalyzed. With the high-speed propeller-driven seaplane now 
an actuality, attention is drawn to the solution of the supersonic 
water-based problem 

It is concluded that adequate design criteria and technological 
experience are available to meet satisfactorily or to exceed any 
aircraft requirement with a suitable water-based configuration 


INTRODUCTION 


Fe ITS INCEPTION, the technical development of 
However, in recent 


the airplane has been rapid. 
years, it has been virtually land-bound, unless one was 
willing to accept the large performance penalties that 
accompanied adaptation to water operation. Because 
of these performance deficiencies, the seaplane for 
many years was restricted to only those functions where 
water-based operation was mandatory and outweighed 
left to the 


land-based airplane the principal task of securing and 


other performance considerations. This 


maintaining the mastery of the air—even over the 
seas—in spite of the hazards of forced landings and the 
restricted freedom of action imposed by the necessity 
for elaborately prepared terminal and auxiliary landing 
facilities. 

Presented at the Aircraft Design Session, Eighteenth Annual 
Meeting, I.A.S.. New York, January 23-26, 1950. 

* Assistant to the Chief Engineer. 


This situation was not always true. Dating from 
the first practical seaplane developed by Glenn Curtiss 
in 1911, seaplane design progressed rapidly through the 
first World War. As a matter of fact, water-based 
development outstripped contemporary landplane de- 
sign during that period to such an extent that, by 1938, 
all long-range transport was still accomplished by flying 
boats and the world’s speed record was held by a sea- 
plane. From a hydrodynamic point of view, this early 
development period culminated in the design and con- 
struction of the famous NC flying boats in 1918. The 
hull of the NC was 45 ft. long with a beam of 10 ft., 
and it was this hull that influenced basic seaplane de- 
sign and hydrodynamic criteria for the following 20 
years, during which period little research was directed 
toward the improvement of hull design. 

During this static period, seaplane development 
centered about a few fundamental design criteria that 
had gradually evolved from trial and error and long 
experience. While it was recognized that basic improve- 
ments in hydrodynamic design were no doubt possible, 
research had lagged to the point where aircraft manu 
facturers were not justified in deviating, nor could they 
afford to deviate, far from established trends. Limited 
progress continued, but this was largely due to costly 
flight-test programs and the incorporation of improved 
air-frame and engine installations. This situation re- 
sulted from the fact that there was no satisfactory re- 
search procedure available for the rapid and accurate 
determination of the complex discontinuous variables 
involved during simultaneous operation in two media 
i.e., air and water. 

The introduction in this country of the dynamically 
similar model technique of research (by the author in 
July, 1938') for the first time presented the airplane de- 
signer with a practical research tool for designing the 
hull of a seaplane which was complementary to, and as 
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Original configuration of the XPB2Y-1 which 
(Bottom). 


Fic. 1 (top). 
gave an operational gross weight limit of 40,000 Ibs. 
The revised X PB2Y-1 hull which resulted in an operational gross 
weight of 86,000 Ibs. 


modern in concept as, his science of aerodynamics. 
With this versatile scientific tool available, extensive 
hydrodynamic research programs were initiated during 
the past war to meet the ever increasing demands for 
increased seaplane efficiency and reliability. For ex- 
ample, the Coronado, a large four-engined 
bomber flying boat, was increased from a maximum 
operational gross weight of 40,000 Ibs. to a weight of 
86,000 Ibs.—over double—solely through the applica- 
tion of an intensive dynamic model research program. 
Fig. | shows the modification that eventually was found 
necessary in this case, and it is obvious that such an 


patrol 


extensive revision of lines would not have been dis- 
covered by the costly time-consuming procedure of full- 
scale cut and try. The model tests, on the other hand, 
indicated a complete solution to the problem in 3 weeks 
of testing. In a like manner, the performance and 
efficiency of the famous Convair Catalina, Martin 
Mariner, Boeing Clipper, and other World War II 
seaplanes were greatly increased at a most critical 
period of the war effort. 

With this outstanding confirmation of the effective- 
ness of the dynamic model in rapidly solving complex 
hydrodynamic problems and in spite of the urgency 
attached to improving existing seaplanes during the 
war, considerable effort was directed toward applying 
this new research technique to the fundamental solu- 
tion of basic hydrodynamic criteria. The war was 
demanding that the range and load-carrying ability of 
seaplanes be improved. The only known method at 
the time was to overload continually these aircraft to 
the point where seaworthiness and water-handling 
characteristics were greatly impaired, attempt to refine 
the existing design still further, and then overload 
again. It was imperative, therefore, that some other 
fundamental means be found whereby the desired sea- 
worthiness could be attained without impairing the 
load-carrying ability while at the same time allowing 
marked improvements in performance, which at best 


was dangerously low. 


CAL SCIENCES—AUGUST, 1950 

In searching for a fundamental concept that would 
allow radical improvement in hull performance, ay 
N.A.C.A. Technical Memorandum? was recalled which 
briefly abstracted some experiments conducted py 
Sottorf in Germany. This report presented meager 
resistance and spray information on a series of high 
length-beam ratio hulls which, when carefully analyzeq, 
held considerable promise for a solution to the problem 
just discussed. The brief summary presented jn the 
N.A.C.A. publication was insufficient, however, for g 
thorough analysis, and the complete reports were 
DVL in Berlin. These reports, which 

received after outbreak 


ordered from 
were actually 


Europe, were completely translated by Convair! and 


of war in 
were found to contain a complete account of the Ger 
man experiments. 

More thorough analysis of these later data confirmed 
preliminary theoretical studies that high length-beam 
ratio hulls held out a definite promise for marked hydro- 
dynamic improvement. However, it was also realized 
that these studies were based on a series of resistance- 
type tank models and that many more data, particu- 
larly on the stability and relationship of secondary de- 
sign parameters, were required before utilization of the 
basic principles embodied in these hulls could be 
realized. 

A Navy research authorization was granted to Con- 
solidated Vultee, and, under the direction of Con- 
vair hydrodynamicists, carefully planned and inte- 
grated dynamic model research programs were organ- 
ized at the N.A.C.A. Towing Tank, Stevens Institute 
Towing Tank, and Convair’s own Hydrodynamics Re- 
This resulted in a new family of 
It was through 


search Laboratory. 
hull forms of superior performance. 
the application of basic design criteria obtained from 
this intensive research program that the XP5Y-1, the 
Navy's newest flying boat, was developed. This high- 
speed turboprop seaplane, shown in Fig. 2, represents 
the first water-based design in this country to deviate 
radically from the pattern established by the outmoded 
NC flying boats of World War I. 

Like the NC boat, the XP5Y-1 has a beam of 10 ft 


Fully loaded, the NC had a gross weight of 28,000 Ibs. 


whereas the XP5Y-1 demonstrates excellent hydro- 





Fic. 2. The Navy’s new XP5Y-1_ turbopropeller-powere 
high-speed flying boat incorporating a high length-beam rat 
hull. 


| 
| 
i 





dyn 
crec 
hull 
inst 
brou 
moc 
plar 
the ’ 
flyin 
rati 
and 
prob 
and 
resul 
decr 
such 
integ 
auto 
plane 
prop 
the ¢ 
takes 
tage 
tiona 
force: 
If s 
did a 
succe 
seapli 
steps 
into ¢ 
tually 
onic ¢ 
sacrifi 
perf | 
desigt 
appea 
appar 
refine 
would 
withot 
liness 1 
therefi 
tion a 
search 
applic: 
result | 
norma: 
form. 


In |i 
the ba 
smooth 
section 
floatati 
of bod) 
persont 


of thes 


uld 


an 
lich 
by 
ger 
igh 
zed, 
lem 
the 
ra 
vere 
ich 
in 
and 


3eT- 


ned 
am 
lro- 
ized 
nce- 
ICU- 
de- 
the 


be 


‘on- 
‘on- 
nte- 
yan- 
tute 
Re- 
y of 
ugh 
rom 
the 
igh- 
ents 
jate 


ded 


() ft. 
lbs., 
dro- 








HIGH-SPEED 


dynamic characteristics up to 150,000 Ibs. The in- 
crease in efficiency of the new high length-beam ratio 


Is does not, however, stop with load. Hydrodynamic 


l 
ne om such as porpoising and skipping, has been 
brought under control through the use of dynamic 
models and no longer restricts the operation of new sea- 
plane designs to a narrow band of stable limits; hence, 
the versatility, safety, and seaworthiness of the modern 
flying boat were greatly increased. High hull fineness 
ratio, which allows the incorporation of easy buttocks 
and sharp waterlines, has done much to minimize the 
problem of spray and greatly increase seaworthiness 
and rough-water performance. This same feature has 
resulted in hulls of low frontal area, representing large 
decreases in aerodynamic drag. In conjunction with 
such design innovations as multicellular watertight 
integrity, propeller-turbine engine developments, and 
automatic mooring and docking techniques, the sea- 
plane designer is now able to put a high-performance 
propeller-driven seaplane into the air which fully meets 
the exacting requirements of modern design and yet 
takes full advantage of the inherent strategic advan- 


tage of water-based aircraft from the  opera- 
tional standpoint of mobility and = dispersion of 
forces. 


If seaplane development were not to falter again, as it 
did after the first World War, it was obvious that the 
successful development of high-speed propeller-driven 
seaplanes would provide only temporary respite unless 
steps were taken to develop the water-based aircraft 
into an efficient and dependable transonic, and even- 
tually supersonic, aircraft. The adaptation of trans- 
onic aircraft design to water-based operation, without 
sacrificing any of the high standards of hydrodynamic 
performance which have been introduced into seaplane 
design during the past few years, posed what at first 
appeared to be an insurmountable problem. It was 
apparent from the outset that extensive aerodynamic 
refinement of known hydrodynamic forms of high quality 
would merely result in poor hydrodynamic performance 
without ever fully obtaining the aerodynamic clean- 
liness required for transonic flight. It appeared logical, 
therefore, to start with an ideal aerodynamic configura- 
tion and, by means of extensive dynamic model re- 
search, develop new hydrodynamic principles and 
applications that, though different in concept, would 
result in the same degree of seaworthiness and stability 
normally associated with the best hulls of conventional 
form. 


In line with the above reasoning, it was obvious that 
“in-flight”’ 

smooth, unbroken form, utilizing high critical Mach 

sections throughout. 


the basic configuration must consist of a 
In order to secure elementary 
floatation, it was necessary to supply sufficient volume 
of body so that engines, air inlets, jet exhausts, and 
personnel would not be inundated. The combination 


of these two fundamental requirements led logically to 


WATER-BASED 


AtRCRAPT 159 





Fic. 3 


A dynamically similar flying model of the first jet- 
propelled, blended-hull configuration to be developed for a trans- 
onic water-based aircraft 


the blended-wing-hull configuration as the most promis- 
ing approach to the basic problem. 

By proceeding in this logical manner, a large number 
of transonic hydrodynamic configurations were de- 
signed and tested as free-body dynamic models—some 
were catapulted, others were towed, and many were 
actually flown, self-propelled with their own jet engines. 
It became increasingly apparent that a true aerody- 
namic form could meet adequately, or in some cases ex- 
ceed, the established design criteria for satisfactory 
The 


series of transonic water-based configurations is illus 


hydrodynamic performance. first in this long 


trated in Fig. 3. Continual refinement has radically 
altered the appearance of these high-speed aircraft over 
that shown in this figure. However, the new research 
philosophy out of which it was conceived is firmly estab- 
lished, and the little dynamic model shown in Fig. 3 
may well represent, in the years to come, the real 
rebirth of sea power. For, with the transonic water- 
based aircraft a practical reality from a technical stand 
point, the revitalized hydrodynamic research program 
at Convair is already well into the problems of water 
based supersonic flight. 

The following sections of this paper will outline the 
development of the dynamically similar model and its 
Within 


security classification, the more important hydrody- 


associated research techniques. the limits of 
namic design criteria that have been obtained from the 


above research will be discussed. 


DYNAMIC SIMILARITY 


Fluid Resistance 

Before discussing the various research techniques, 
equipment, and design criteria that have been de- 
veloped around the use of dynamically similar models, 
it is believed desirable to review briefly the basic theo 
retical background upon which this versatile applied 
science is based. Classical hydrodynamics considers 
the perfect fluid that is defined as incompressible and 
i.e., there is no resistance to shearing be 
Water is, all 
has viscosity and other 


nonviscous 


tween the elements. for purposes, in 


compressible; however, it 
characteristics such as surface tension which require 
special attention and cannot be neglected. The general 


resistance equation for the force acting on a body in 
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motion, partly or wholly immersed in an imperfect 
fluid such as water, can be expressed in the following 
form: 


Vv? pVL Viv VEL 


’ , ,’ , » ,’ (1) 
co pp @ & 4h y Ob 


R = (pV°L2)f | ¢, 


where the variables are density of the fluid, p; velocity 
of the body, V; linear size of the body, Z; trim or 
angle of attack, 7; gravity, g; coefficient of viscosity, 
“; compressibility of the fluid, ’/a; surface roughness, 
1/L; texture of the fluid flow or turbulence, v/V; sur- 
face tension, y; and fineness or aspect ratio, L/d. 

The term p/°L? is the main term of the expression 
and is recognized as the force due to density. We are 
accustomed to the use of a single nondimensional coef- 
ficient of resistance instead of the complex function in 
the brackets. It must be borne in mind, however, that 
such a single coefficient, when used, varies with, and 
attempts to be, the equivalent of the terms listed. 
Some of these terms may be negligible, but all are sub- 
ject, to some extent, only to experimental determina- 
tion. 

In hydrodynamic experimental testing, it is neces- 
sary to determine f;(7) or the variation of resistance co- 
efficient with trim. Because of the gravity effects as- 
sociated with wave-making, the term fo(V’*/gL), or 
Froude Number, must be met. Because of the depend- 
ence of frictional resistance on viscosity, f3(o VL/), or 
Reynolds Number, must be calculable or an allowance 
made. As we are dealing with incompressible water, 
the compressibility term f;(l’/a) may be neglected. If 
tests on two models do not agree, the surface roughness 
fs(//L) may not bear the same relation in the two 
eases. This function has a strong bearing on frictional 
resistance as does fs(v/ V), which is the stream turbu- 
lence as represented by the ratio of the average lateral 
turbulence velocity to the measured axial velocity. 
Surface tension represented by f7(V?L/y) can, in most 
cases, be neglected, except possibly in accurate studies 
of spray formation. Experiments have been conducted 
where this term is eliminated by the introduction of 
Aerosol, or some other commercial wetting agent, to 
the water. Finally, it is apparent in aerodynamics, and 
is becoming increasingly so in hydrodynamics, that 
fineness or aspect ratio, represented by fs(L/b), cannot 
be neglected. In hydrodynamics, this becomes the ra- 
tio of the planing bottom fineness or length to beam 
ratio. These are by no means all of the dimensionless 
combinations that may be written, but they com- 
prise what are usually considered the most impor- 
tant. 

Using this typical example of the numerous factors 
influencing fluid resistance, three important facts can 
be derived which apply equally well, though varying 
in detail, to all other quantities for which scale repro- 
duction is desirable: 


(1) Theoretically, no model test can completely rep- 
resent a full-scale condition unless all of the dimension. 
less ratios are held constant. 

(2) It is impossible to hold all of these ratios constant 
at the full-scale value in a model test, for some of them 
are contradictory. 

(3) The dimensionless ratios may not have equal 
weight, but theory does not show it. 

Experience tells us, however, that the angle of attack 
or trim is the most important of all the terms given and 
that for hydrodynamic work the Froude Number should 
be held constant. It is obvious from this discussion 
that experimental research on models, in an attempt to 
determine or predict full-scale behavior, is not a cyt. 
and-dried procedure but is a series of compromises 
based on judicious judgment and experience. Even 
the meticulous matching of Froude Number does not 
assure similarity. For the lower speed régimes where 
wave-making predominates, it will give close approxi- 
mations, but, as the speed increases to high-speed plan- 
ing on the water, the function of Reynolds Number be. 
comes more and more predominant. While the above 
discussion appears to present a pessimistic picture of 
experimental model testing, it should be pointed out 
that, by comparison, the theoretical or analytical solu- 
tion to hydrodynamic problems, particularly those 
concerned with stability, is practically hopeless, with 
no prospect of improvement. On the other hand, 
rapid progress has been made toward the solution of 
the problems associated with dynamic similarity, and 
it appears that the full application of this relatively in- 
expensive and direct approach to the solution of dis- 
continuous complex variables is virtually unlimited. 
Later sections of this paper will discuss in more detail 
the specific deviations from theoretical similitude and 
the remarkable degree of accuracy that can be ob- 
tained from the modern dynamic model when it is prop- 
erly designed and operated. 


Froude’s Law of Comparison 

In the previous section on fluid resistance, it was 
stated that Froude Number, |/gl, must be held con- 
stant in all cases of experimental hydrodynamics when 
anything but the specific full-scale configuration itself 
is being investigated. It was further stated that, with 
calculated reservation, the holding of this function 
constant would most nearly result in model exper 
ments producing quantitative data directly commen- 
surate with full scale. As these statements contain the 
heart of all hydrodynamic research and analysis, it 1s 
advisable to become more familiar with the character 
and limitations of this function. 

In 1867, William Froude had constructed what 1s 
generally conceded to be the first modern towing tank, 
and, by 1879, when he died, he had contributed vol- 
umes to the science of naval architecture. Of greatest 
significance was his famous Greyhound tests, where 4 
172-ft. boat named the “Greyhound” was towed by 
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By duplicating the full-scale tests in 


H.M.S. ‘‘Active.”’ 
owing tank with a '/16-scale model and by applying 
appropriate coeflicients of friction to both full scale and 
rom data previously determined in the tank, he 
able to separate that portion of the total resistance 
By 


analysis of these data, Froude determined that, if any 


his t 


model, f 
was 
due to gravitational wave formation and eddies. 
full-scale speed was divided by V 16 (the square root of 
the model scale), the wave and eddy resistance of the 
model at that speed was equal to the corresponding re- 
sistance of the full-scale ship divided by 16% (the model 
This relationship Froude called the 
By referring to Eq. (1), we 


scale cubed). 
“Law of Comparison.”’ 
see that the second function in the bracket, I’*/gZ, pre- 
viously called Froude Number, demands that V’ vary as 
the square root of the linear dimension Z if model com- 
parison is to be valid. 

As pointed out, this relationship applies only to re- 
siduary resistances and wave formation. If we were 
limited to ideal fluids, Froude’s Law of Comparison 
would be sufficient. However, the functions of vis- 
cosity, surface tension, and other properties of water are 
such that any attempt to bring them into a law of sim- 
ilarity requires the introduction of Reynolds Number 
as well. The conditions where Reynolds and Froude 
agree exist only at full scale. If the frictional compo- 
nent is small, as is usually the case in the larger tank 
models of conventional speed seaplanes, it is possible 
to neglect its effect and proceed under Froude’s com- 
parison and scale total resistance directly to the scale 
cubed. In most cases this will result in a degree of full- 
scale conservatism which will cover items of surface 
roughness and protuberances not usually present in the 
model tested. If, on the other hand, we are dealing 
with small models of extremely high-speed planing 
may be small 
the fric- 
at the 
theory. 


craft, gravitational resistance or 


negligible, and large errors will 
not 


result if 
tional resistance is separately treated 
proper Reynolds Number as in aerodynamic 
Fortunately, the majority of critical items of hydro- 
dynamic performance, such as stability, spray forma- 
tion, seaworthiness, etc., occur at low speeds in the 
range of transition from displacement to planing and 
thereby closely follow Froude’s Law. For emphasis it 
should be pointed out again that one will find all of the 
standard coefficients and parameters of seaplane design 
based on Froude’s Law of Comparison and, unless the 
research engineer is fully aware of the limitations and 
treatment of this function, that serious inconsistencies 
will result, particularly as the trend for higher speed 
water-based aircraft continues. 


Dimensional Analysis 


Inasmuch as experimental research is the foundation 
of all applied hydrodynamics for the seaplane designer, 
it is necessary that the engineer understand the dimen- 
sional test equivalents for all factors encountered full 
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scale. 
the components of hydrodynamic resistance and how 
In this section, with 


In the two previous sections, we have discussed 


this force varies with linear scale. 
Froude’s Law as a basis, we will develop the relation- 
ship of many other interrelated physical functions of 
dynamics and mechanics to linear scale. Keeping in 
mind the inherent, but calculable, discrepancy of 
Froude’s fundamental Law of Comparison, we can es- 
tablish mathematically how every other physical 
quantity must vary to keep the total dimensionally 
correct. 

In the early days of seaplane research and ship de- 
sign, the designer was primarily occupied with the re- 
sistance of his hull forms, and, consequently, the basic 
law of comparison, as stated by Froude, dealing with 
the relation of resistance to speed with varying scale 
was adequate. However, with increasing knowledge of 
hull form and the rapid development of powerful en- 
gines, resistance was gradually, but steadily, subordin- 
ated to the more critical studies of dynamic stability 
on the water and the factors affecting spray formation 
and seaworthiness. As stated previously, the most di- 
rect approach to these problems is through the use of 
dynamically similar scale models. To obtain accurate 
results for these complex dynamic problems, it is nec- 
essary for the models to be dimensionally correct in all 
respects to the full-scale prototype—i.e., they must not 
only be geometrically to scale, such as a wind-tunnel 
model or resistance-tank model, but it is also required 
that they have, among other things, gross weight, in- 
ertia, power, accelerations, and all aerodynamic forces 
and moments to scale. It is obvious, therefore, that, 
if this can be accomplished throughout for every factor, 
we have in effect a flying miniature of the full-scale air- 
plane that will perform every maneuver of the full-scale 
aircraft and at a rate of movement directly to scale. 
In effect, the dynamic model becomes a complex integrat- 
ing mechanism that automatically picks up every known 
or unsuspected force, in the proper magnitude, point of 
application, direction, and sequence; integrates all these 
reactions instantaneously; and provides the observer with 
the resultant motion and rate. Even if there were no 
unknown transient forces, the task of integrating all 
known forces in a complex dynamic reaction by analy- 
tical means, for just one speed point, is enormous. It is 
this goal of tremendous simplification of integrating all 
forces on a free body which has made the problems as- 
sociated with the attainment of such a model seem in- 
consequential by comparison. 

Because the dynamically similar flying model repre- 
sents the ultimate in experimental research at reduced 
scale, we shall develop the principal scale relationships 
involved in its design and analysis. It is apparent 
that, if sufficient basic factors are dimensionally cor- 
rect, other minor dependent variables will automatically 
follow. Therefore, the derivations and relationships 
that follow are the principal functions involved and 
should be sufficient for most analyses. The engineer 
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will find it a simple matter to derive certain other func- 
tions not specifically listed in this compilation. 

In the following derivations we will consider that the 
symbol for scale, A, represents the whole number 
if \ = S then the linear scale is '/gor\~'. As the linear 
or geometric scale, \~', is usually given, it is desired to 
get all other physical relations in terms of this one value. 
It is obvious, therefore, that, if the linear scale is a ratio 
of lengths L, then Z varies directly as 1/A or A~', 
written L « X~!. Following this procedure, it is plain 
that an area is made up of a length times a length or L?. 
Therefore, area, or L’, varies as the linear scale squared 

i.e., area « X~*. In a similar manner, it may be 
reasoned that a volume or mass is an area, L?, multi- 
plied by thickness or height, Z, giving L*. Hence, 
volume, weight, or, as Froude determined, force « 
A~*. As the moment of a force is that force multiplied 
by an arm, L, we may extend our reasoning to show that 
Similarly, the moment of inertia 


i.€., 


moment « Le orx-* 
is a mass multiplied by the arm squared, or moment of 
inertia «)\~, 

Since our system of similarity is based on Froude’s 
Law of Comparison and since we have seen from Eq. 
(1) that this expression depends upon the speed varying 
with the square root of the linear dimension, it follows 
that velocity «VL or i As distance is a linear 
dimension, L, and velocity varies as the V L, then time, 
which is distance divided by velocity, or L/V L, must 
also vary with VL. Hence, time « ) Now that 
we have the basic variations of mass, length, time, and 
velocity, it is a simple matter to substitute in the ex- 
pressions for any physical function and derive its varia- 
tions with linear scale. For instance, revolutions per 
minute is revolutions, which are nondimensional, di- 
vided by time, or | VL, which gives us the relation 

« Likewise, acceleration is feet per second 
squared, or L (V L)?, giving acceleration « 1, or unity; 
This means that 


r.p.m. 


hence, it becomes nondimensional. 
all model linear accelerations will be identical to the 
full-scale accelerations in magnitude. On the other 
hand, if we consider angular motion, we find that angu- 
lar velocity is nondimensional radians divided by sec- 
onds, or 1/V L, giving, as in the case of r.p.m., w « d'/*, 
As angular acceleration is radians per second squared, or 
1/L, we have simply a«X for this quantity. Whereas 
we found that linear acceleration is identical to both 
model and full scale, we see that angular acceleration 
will be the whole number ) times as great in the model 
as the full-scale value. The fact that linear accelera- 
tions do not vary with scale is fortunate, inasmuch as 
the value g, which is the acceleration due to gravity, is a 
constant over which we have no control and corrections 
would be extremely difficult at best, if not impossible. 
Continuing our derivations into the more complex 
functions, we recall that power is defined as the work 
accomplished per unit of time, where work is the prod- 
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uct of a force times the distance through which it acts 
Following the previous line of thought, we can consider 
that force, L*, times the distance of action, L, Causes 
work to vary with the fourth power of linear scale 

i.e., work « \~‘4 and power will therefore be L4/y L 
which gives L’” or power «X That these rela. 
tionships are dimensionally correct can be quickly 
checked by substituting these derived values into any 
formula defining some nondimensional coefficient and 
thereby demonstrate that the numerical value of the 
coefficient does not change with scale. For this dem. 
onstration we can pick such an expression as the fo. 
lowing, which defines the well-known nondimensiona| 
power coefficient, Cp: 


5 X& 10” X b.hp. 


Cp 
(VN)? xX (D)5 


(full scale) 
where b.hp. is the brake horsepower, JN is the propeller 
revolutions per minute, and D is the propeller diameter 
in feet, all full-scale values. From our previous dis- 
cussion we know that the Cp for the model would be 


5 X 10" X (b.hp. XK A~’’) 
(N X d”)® X (D X dN) 


C P(model) 
and, solving for the variation of A, we find 


: ,-" nN 
C; a = =] (4 


(model) he 4 d ) d 


Therefore, the ratio of Cp fut seatey tO Cr¢modet) iS unity. 
For convenience, Table 1 is presented which sum- 
marizes the principal relationships in condensed form, 
giving, in addition, a typical set of values for an as- 
sumed value of \ = 8 giving a linear scale of !/s. 


TABLE 1 
Dimensional Conversion for Linear Scale 


General 1/5 Scale, 
Unit Conversion 1=8 
Linear dimensions 7! 1/8 
Area 2 1/64 
Volume, mass, force \ 1/512 
Moment r 1/4,096 
Moment of inertia 75 1 /32,768 
Linear velocity N77? 1/2.83 
Linear acceleration Constant l 
Angular velocity ee 2.83 
Angular acceleration 8 
Time \ 1/2.83 
R.p.m. x 2.83 
Work r 1 /4,096 
Power r 1/1,446 
Wing loading d 1/8 
r 2.83 


Power loading 


The interesting fact to note in closing this brief dis- 
cussion of dimensional analysis is that all of the factors 
listed in Table 1 have been experimentally checked on 
numerous occasions by constructing and testing models 
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xisting airplanes for the purpose of positive correla- 


of e 
For instance, in the case of the Navy's XP4Y-1, 


tion. 
\/,-scale, radio-controlled, dynamically similar model 


a 
4) was constructed and thoroughly correlated with 


(Fig. 
iain airplane, as well as wind-tunnel and tow- 
ing-basin, tests. With accurate scale propellers set at 
the actual full-scale blade angle and with the r.p.m. 
adjusted to 2.83 times the full-scale value (see Table 1), 
the engine power was measured on a dynamometer and 
was found to be 1/1,446 the full-scale value, and the 
thrust developed was 1/512. With this power and a 
model weight of 1/512 full scale, the model was found 
to become air-borne at a time and speed equal to 1 /2.83 
that observed during flight tests. Linear acceleration 
at the hump and getaway was found to be the same on 
model and full size. Additional description and results 
of these correlation studies, particularly with regard to 
the more complex functions involved in the stability 
derivatives, are covered in more detail in reference 4. 


RESEARCH TECHNIQUES 


Towing Tank Testing 


The oldest and most usual form of dynamic model 
testing isin a towing tank. Prior to the introduction of 
dynamic models, these towing tanks were used primar- 
ily to tow solid resistance models through the water, 
much as William Froude did in the late 1800's. These 
resistance models were attached to a towing carriage 
through a dynamometer linkage such that resistance 
and moments at fixed trims or free to trim resistance at 
various speeds could be recorded. The results were 
similar to those obtained in a wind tunnel, inasmuch as 
forces and moments due to form only were the princi- 
pal products of the testing. 

With the introduction of dynamically similar scale 
models in 1938, it became necessary to revise all the old 
procedures and equipment of tank testing in order to 
take advantage of the unlimited opportunities afforded 
to investigate the dynamic effects of hydrodynamic sta- 
bility and its interrelated influence on the associated 
aerodynamic parameters of the aircraft as a whole. 
The towing tank became a free-stream turbulence wind 
tunnel, as well as a water channel, and the models had 
to be given as many degrees of freedom as physically 
possible in order that resultant motion would be unre- 
strained. Inasmuch as pitching and heaving motions 
associated with the cyclic hydrodynamic stability prob- 
lem known as porpoising were of immediate concern and 
because of the fact that the towing channels were nar- 
row, precluding lateral deviations from a course down 
the center of the tank, it was customary to mount the 
models so that they had freedom in pitch and rise but 
were restrained in yaw, roll, and freedom along the 
longitudinal axis. 

This is accomplished in a manner illustrated in Fig. 
5, which shows the '/1 scale powered model of the 
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The first completely free-body, self-propelled dynami- 
A '/s 


Fic. 4. 
cally similar model to be flown by remote radio cortrol 
scale model of the XP4Y-1. 


XP5Y-1 attached to the towing carriage of N.A.C.A. 
Tank No. | at Langley Air Force Base, Va. The rec- 
tangular towing staff is pivoted at the center of gravity 
of the model and allows freedom in pitch while restrain- 
ing the model in yawand roll. The staff is machined on 
all faces and rides in a roller cage attached to the car- 
rier. The roller cage restrains the staff in yaw and 
holds it to a vertical position but allows complete free- 
dom in rise. It will be further noted in Fig. 5 that the 
roller cage is mounted to the carriage truss through a 
dynamometer linkage so that resistance and/or thrust 
may be measured during the run. 

Early experiments were conducted with unpowered 
models, and the towing force or thrust was transmitted 
to the model through the towing staff to the model cen- 
ter of gravity. While this procedure was providing 
many results theretofore unattainable, true dynamic 
similarity was not obtained, and every effort was made, 
in conjuncticn with the N.A.C.A. staff, to improve con- 
tinually the efficiency of dynamic models and their 
technique of testing. In 1940, it became obvious that 
the models should incorporate scale power and running 
propellers in order to further enhance and refine the ex- 





10 scale XP5Y-1 dynamically similar powered 
model in the N.A.C.A. towing tank at Langley Air Force Base. 


Fic. 5. The 
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The '/s scale Coronado powered dynamically similar 
towing tank model. 


Fic. 6. 


cellent data being obtained. Convair was given a con- 
tract by the Navy to design and install variable fre- 
quency power equipment in the N.A.C.A towing tank, 
along with associated development of lightweight elec- 
tric motors, adjustable-pitch metal propellers, and 
other equipment involved in going to power operation. 
A 4/s scale dynamic model of the Convair Coronado, 
shown in Fig. 6, became the first model to be tested with 
this equipment. 

The absence of complete freedom of motion still re- 
stricted the scope of testing and precluded any study of 
directional problems. Various attempts were made to 
incorporate additional degrees of freedom of a limited 
nature, mostly without success. One important result 
of these studies was the development of a rail system 
for the roller cage which allowed an appreciable amount 
of fore and aft freedom. With careful control of model 
and carriage power, it became possible to operate the 
model truly self-propelled for short periods of time. 
However, except for certain special constant-speed runs, 
the basic towing tank technique consists of towing a 
powered model, rigidly attached to the towing carriage 
through the center of gravity pivot and having freedom 
in only pitch and rise. 

There are four basic testing techniques utilized in 
tank-testing of dynamic models which will be briefly dis- 
cussed : 

(1) Constant Speed Run. 
not a dynamically similar function inasmuch as the im- 
Never- 


The constant speed run is 


portant acceleration parameter is neglected. 
theless, this technique is valuable for establishing cer- 
tain fundamental hydrodynamic criteria such as trim 
limits of stability, resistance, underwater flow photo- 
graphs, and spray patterns. Because of the ease in du- 
plicating test conditions and the time available for ob- 
servation, the influence of extremely small variables 
may be detected. 

(2) Accelerated Run.—The 
the model is accelerated at a constant rate from stand- 


accelerated run, where 
still to getaway, most nearly simulates full-scale opera- 
tion. While the actual take-off run of a seaplane is not 
made at a constant acceleration rate, it is not feasible 
to vary this factor in the towing tank because of the 
high inertia and complicated speed control of the tow- 
ing carriage. It is customary to set a constant rate 
equal to the acceleration of the aircraft during the criti- 
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cal speed range at the hump and maintain this rate gop, 
stant to getaway. An even more serious drawback t) 
accelerated runs in the towing tank is the inability of 
any present-day equipment to match the acceleration. 
being studied for modern high-speed water-based air 
The ac- 
celerated run technique ts used in all cases for evaluating 
the stability range, aerodynamic control, and wm 


craft, which are approached !/3g in some cases. 


characteristics of a specific design configuration. 

(3) Landing Run.—This technique involves getting 
the model up to stabilized flight as rapidly as possible 
in order that a landing from the flight attitude may ly 
made. The majority of these landings are made with 
the model attached to the towing carriage, which is de. 
celerated from flight speed at a constant rate. Here 
again, in most cases, it is not feasible to decelerate the 
heavy towing carriage at a rate commensurate with ful] 
scale, and, being rigidly attached, the predetermined 
rate of the carriage is impressed upon the model regard. 
less of the varying decelerating forces it may experience 
In certain cases where the model is small, it has been 
found practical to launch the model from the carriage in 
free flight. 
namically similar and, consequently, are much more de- 


Landings under these conditions are dy- 


sirable in that they overcome the above objecticns to 
the rigid attachment. Obviously, control and instru. 
mentation is complicated considerably by this latter 
technique. 

(4) Generalized Similarity. 
procedure is commonly called the “Stevens Method,’ 


This last towing tank 


inasmuch as it represents the unique testing technique f 


employed by the Experimental Towing Tank at Stey 
ens Institute of Technology. Fig. 7 illustrates the 
Stevens towing carriage with a dynamic model at 
tached. Here it will be noted that the aerodynam« 
parameters and derivatives are introduced through 
calibrated weights, springs, and dash pots. While all 
of the full-scale parameters must be known in order t 
set up their equivalents on the model system for specific 
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Fic. 7. 
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tests, this method does prove versatile for generalized 
problems where it is desirable to introduce the aerody- 
namic derivatives as independent variables. It is 
dificult to determine a specific evaluation by this 
method, such as is possible with the accelerated run 
technique; however, for rapid comparative study of 
many variables, this technique is superlative. 

Steady progress has been made in the continual ef- 
fort to minimize or eliminate inaccuracies in the scale- 
testing of dynamically similar models in the towing 
tank. Obviously, the ideal situation is to have a fully 
instrumented, self-propelled free body under the posi- 
tive control of an operator. Such a model would in- 
herently produce the proper variation of acceleration 
and unrestrained motion in space. In spite of the 
physical limitations and compromises of the towing 
tank, this equipment produces accurate results in the 
hands of experienced analysts who fully recognize, and 
allow for, the deviations from dynamic similitude. 


Free-Body Testing 


During World War II, the facilities at the N.A.C.A. 
and other towing tanks in this country became crowded, 
and many high priority developments were suffering 
long delays. To alleviate this condition and at the 
same time obtain for the first time the optimum condi- 
tion previously mentioned (that of having complete 
free-body operation), Convair started to develop a sys- 
tem of free-body dynamic model research. The heart 
of this development is the large, natural outdoor towing 
basin available in the form of the U.S. Naval Training 
Station Estuary, illustrated in Fig. 8. This ideal test- 
ing facility is an arm of San Diego Bay, and it is com- 
pletely contained in the center of a military reservation, 
which provides excellent security protection. Located 
adjacent to Convair property, its smooth sand beaches 
and calm protected water allow dynamic models to be 
operated by the contractor at his plant, thus ensuring 
uninterrupted research and development. 

As discussed previously, the problems associated 
with complete free-body operation are numerous. The 
models must be propelled and be under positive control, 
and, to be of any value, technical test data must be ac- 
curately and rapidly available. This has resulted in 
the development of many new pieces of test equipment 
and of new research techniques. However, as in the 
early development history of the dynamic model itself, 
the potential rewards in the form of versatility, speed of 
results, and accuracy outweigh the purely mechanical 
problems associated with the attainment of this goal. 
For instance, unlike the inherent restrictions of the tow- 
ing tank, the free-body model produces unrestrained re- 
sultant motion about all axes, allowing complete stud- 
ies to be made of directional stability, water-looping, 
low-speed maneuvering and control, drifting in seas, 
course-keeping, and other specified problems, in addi- 
tion to the standard straightaway run. Take-off time 
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Fic. 8. The San Diego Naval Training Station Estuary used 
by the Convair Hydrodynamics Laboratory for open-water free- 
body testing of dynamically similar models. 


and distance under all conditions of varying seas, wind, 
and loading can be determined immediately without re- 
course to the complicated take-off integral. Perhaps 
of greatest value, however, is the ability to operate in 
an environment as natural as nature itself. The tow- 
ing tank has always presented a problem when evalua- 
tion or study of the influence of rough or mixed sea con- 
ditions is desired. Inasmuch as an actual seaplane 
rarely operates in a flat calm, the entire field of typical 
service operational problems has received but cursory 
examination. 
always meet the model head-on, which is the most 
hazardous and least recommended procedure for land- 
To make mat- 


The waves propagated in a towing tank 


ing a full-scale seaplane in rough water. 
ters more difficult, the assisting stiff head wind invari- 
ably present in head-on seas of the type generated in a 
towing tank is not duplicated, which allows the model 
to encounter these -heavy seas with extremely excessive 
speeds. On the other hand, at the outdoor estuary 
the gradually rising mid-day and afternoon breeze 
provides a daily testing period from 3 to 4 hours of flat 
calm in the morning, for precision laboratory type of 
testing, to scale head winds of 30 to 40 knots in the early 
afternoon. By judicious selection of time of day, any 
degree of full-scale operational conditions may be ob- 
tained. Likewise, heading and sea condition may be 
varied at will. With the desired wind and heading se- 
lected, artificial seas of any desired pattern and com- 
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The '/s scale X P4Y-1 radio-controlled model and the 
seven-channel proportional control transmitter developed by the 
Convair Hydrodynamics Laboratory. 


Fic. 9. 


plexity can be generated in the path of the model by ex- 
ternal sources. With the models under remote control 
and free of all restraint, valuable information can be 
obtained concerning characteristics most difficult and 
hazardous to secure full scale and impossible to secure 
in the towing tank. 

There are four basic open-water, free-body research 
techniques currently being used to develop high-speed 
water-based aircraft; these will be discussed as follows: 

(1) Radio Control.—As stated previously, to take full 
advantage of the great potentialities of the free-body 
research technique, it is necessary to have a self-pro- 
pelled dynamically similar model under the positive 
control of an operator with an accurate system for re- 
cording test data. The precise development of any one 
of these requirements is a major undertaking when the 
small available weight and space allowances are consid- 
ered. 

In June of 1943, development was started of a !/s 
scale dynamically similar model of the Convair XP4Y-1 
twin-engined flying boat to be remotely controlled by 
positioning, multichannel radio. Scale requirements 
called for a gasoline power plant of 1'/2 hp. at 4,250 
r.p.m. and not to exceed a weight of 4 Ibs. Ohlsson & 
Rice Manufacturing Company, of Los Angeles, pro- 
duced this engine in 4 months after date of contract, and 
the first problem, that of self-propulsion with scale 
power, was solved. 

To attain precision positive control, the Convair Ra- 
dio Laboratory developed a seven-channel, positioning 
radio transmitter, lightweight receiver, and associated 
servomechanisms. The detail technical development 
of this interesting and elaborate system is fully de- 
scribed in reference 5 and will not be repeated here. 
With self-propulsion and a precise remote control sys- 
tem available, attention was directed toward developing 
an accurate phototheodolite tracking and film analyzer 
system for recording the motions and accelerations of 
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the free body in space. The first completely instry. | 
mented, radio-controlled free-flight dynamic model was 
successfully demonstrated by Convair in August, 1944 
just 14 months after initiation of the program, This 
original model is shown in Fig. 9 and a complete descrip. 
tion of its development and correlation with full-scale 
flight tests and the towing tank is given in reference 4 
Following the successful demonstration of precision 





free-body control and analysis, a concerted program oj 
research was initiated in an effort to overcome the ge. 
rious seaplane problems discussed in the introduction tg 
this paper and, once again, to place efficient, high-speed 
water-based aircraft in the forefront of defense plan 
ning. <A total 
trolled research model configurations were tested to de 


of 27 dynamically similar radio-cop. 
termine the fundamental factors affecting hydrody. 
namic efficiency. Out of this intensive study came the 
full development of the remarkable high length-beam 
ratio hull, which resulted in the Navy's new turbopro- 





10. The '/j seale XP5Y-1 radio-controlled free-flight 


model sitting in front of the full-scale mock-up. 


Fic. 





The XP5Y-1 free-flight model just prior to beiti 
released for a self-propelled radio-controlled test flight at the ¢ 
tuary. 


Fic. 11. 
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peller-powered XP5Y-1 previously discussed and shown 
in Fig. 2. 

The '/1 scale radio-controlled dynamically similar 
model of the XP5Y-1 is shown in Figs. 10 and 11. This 
model has made over 2,000 test runs during the past 3 
years and has completed the full hydrodynamic flight- 
test program scheduled for the full-scale aircraft. In 
addition to performing the functions normally expected 
of the full-scale aircraft, the model has thoroughly 
explored régimes, attitudes, and conditions considered 
far too hazardous to risk human life and costly equip- 
ment to investigate during flight testing. Many of 
these conditions, once investigated by the model, have 
proved to be satisfactory and thereby serve to enhance 
the operational safety and utility of the aircraft. Those 
that fail result in relative minor costs to repair the 
model, no one is injured, and valuable experience is 
gained. ‘ 

Since the first radio-control system was developed in 
1944, there have been many improvements in electronic 
control, recording equipment, and operating technique. 
For all precise, laboratory-type analysis of straightaway 
control and stability in calm water, it is now customary 
to utilize a fast-acting bang-bang type of radio control. 
As it is desired to accelerate to getaway on a straight 
course in smooth water with the principal controls 
treated as independent variables, the elevators, rudder, 
flaps, and throttles are fixed at their test settings and 
the bang-bang control is connected to only the ailerons 
and ignition so that slight corrections to maintain 
wings level and the “blipping”’ of port or starboard en- 
gines to maintain course are the only variables involved. 
Just after getaway, all ignition is cut and the run is 
terminated. For the elaborate evaluations of 
seaworthiness, maneuverability, rough-water take-offs, 
and flight, the original seven-channel proportional con- 
trol system is used which provides simultaneous preci- 


more 


sion control of the elevator, rudder, ailerons, and inde- 
pendent throttles, in addition to selective ignition cut- 
off. 

(2) Catapult Launching. 
powered model is used for the majority of conventional 
hydrodynamic tests, it is not desirable or necessary to 
employ these elaborate models for routine investiga- 
Through the 


While the radio-controlled 


tions of landing stability and impact. 
use of a catapult and open-water conditions, where the 
heading and air speed relative to the water may be se- 
lected to match the conditions under study, the engines 
and radio equipment can be removed and their weight 
replaced with more elaborate recording equipment. 
For these tests the aerodynamic controls are preset 
prior to each run. Through the use of an accurate 
velometer to record wind velocity and direction, of 
calibrated catapult spring tension, and of a little exper- 
ience in presetting the controls, it is possible to execute 
precise landings at any attitude and rate of sink within 
the capabilities of the aircraft configuration being 
tested. 


WATER 
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The XP5Y-1 free-flight catapult model during landing 
stability tests at the estuary. 


Fic. 12. 


Fig. 12 shows a typical catapult launching of the ' \o 
scale XP5Y-1 catapult model. The equipment shown 
is a standard 30-ft. A-2 target drone catapult, which is 
capable of launching a 100-Ib. model at a flight speed of 
60 ft. per sec. The catapult is portable, and its head- 
ing and height above the water may be varied at will. 
For these tests, the models carry elaborate internal in- 
strumentation, in the phototheodolyte 
tracking equipment previously mentioned. When the 
lanyard is pulled, an inertia switch automatically turns 


addition to 


on a 16-mm. gun camera mounted within the model. 
This camera records a continuous record of trim, air 
speed, relative wind, bow, step and sternpost contact 
) 


with the water, and impact acceleration. Fig. 13 is a 


picture of the internal instrumentation installation. 

The free-body catapult technique is used to obtain 
smooth- and rough-water landing stability, aerody- 
namic ground effect, landing impact, and water-looping 
characteristics. Here again, particularly in the study 
of rough-water impact and water-looping, attitudes and 
conditions that would never be attempted full scale can 
be thoroughly studied at no risk. Only in this way can 
the full capabilities of a new design be determined. 

(3) Bridle Tow. 
good advantage in all those cases where it is desired to 


The bridle tow technique is used to 


obtain resistance data, position the model close to the 
observer for detail study, determine preliminary sta- 
bility and spray information from small survey models 
incapable of carrying power plants or radio, and con- 
stant speed running under controlled conditions. This 
procedure most nearly approaches the towing-tank 
method except for the fact that selection of heading, 
wind, and sea is still available to the operator. 

Fig. 14 illustrates a typical bridle tow installation for 
a small jet-propelled research model. The launch 
tows the model through the thrust line from the end of 
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Fic. 13. Internal instrumentation in the X P5Y-1 dynamically 
similar free-flight model. 


a towing boom that contains a resistance dynamometer. 
The instrumentation box on the launch contains the 
resistance and water-speed indicators, which are re- 
corded by a camera. Just below the recording instru- 
ments, the model being towed is framed in an aperture 
that allows the image of the model, the shore line, re- 
sistance, and water speed to be recorded simultaneously. 
In this manner, complete trim tracks with a record of 
spray and resistance at any speed may be determined. 





Fic. 14. The bridle tow technique developed by the Con 
vair Hydrodynamics Laboratory for the testing of small un 
powered survey models of supersonic water-based aircraft prior 


to self-propelled tests 


Also, the condition of the sea at the moment is perma. 
nently recorded for future reference or evaluation, 
This system is extremely versatile and is used extens 
ively during early survey investigations of a radical new 
design prior to constructing the large and more elabo- 
rate powered models. 

(4) Free Launch.—The testing procedures discussed 
thus far have dealt primarily with gasoline engine pro 
peller-driven models or models containing no power 
plants at all, such as the catapult and bridle tow tech 
niques. To date, no large jet-propelled radio-controlled 
models have been tested comparable to the large pro 
peller-driven types previously discussed. This is dye 
to several causes: (1) suitable jet units for producing 
sufficient scale thrust have not been available; (2) the 
study of jet-propelled water-based aircraft is so new and 
the problems so diversified that emphasis has been put 
on a large number of small survey models in an effort 
to establish basic design criteria; and (3) the inherently 
greater speeds of the jet-propelled aircraft have posed 
the tremendous problem of speeding up pilot and servo 
reaction times beyond the point that is now critical, 
It will be recalled that, in dynamic similarity, time var 
ies as the square root of the scale, so that on a '/;9 scale 
model events occur roughly three times faster than the 
equivalent full-scale event. Operation in confined 
areas with extremely fast aircraft rapidly multiplies the 
problems of control reaction. 

It was inevitable, however, that, if water-based air 
craft were to keep pace in the aviation field, means to 
obtain powered dynamically similar research data, 
comparable to that being secured on conventional-type 
seaplanes, must be developed. The first step in this 
direction was the availability of an efficient, lightweight 
pulse-jet engine. Aeromarine Company of Vandalia, 
Ohio, produces such engines in 4'/2 and 30 Ibs. of thrust 
per unit ratings. For additional versatility, the 4'/.-lb 
thrust units may be manifolded together to produce an 
S'/5-lb. thrust unit. 

With suitable power plants available, detailed stud 
ies of the many mounting, manifolding, cooling, and as 
sociated problems were started. Obviously, it was de 
sirable to check out these many jet-power development 
problems with as small and economical a model as pos 
sible in view of the attrition rate anticipated. It was 
at this point that the free launch technique was evolved 
It was reasoned that excellent results are attained from 
catapult launchings with preset controls and that with 
the accelerations anticipated the take-off run would 
closely approximate a catapult launching. In_ the 
original attempts at free launch, the fuel was metered 


' 





— 





closely, and progressively longer bursts were permitted 


as the trim and balance of the model proved satisfac 
tory. This technique resulted in some exciting, inad 
vertent uncontrolled take-offs and climbs to high alti 
tudes when combinations of factors became ideal; and, 
as a consequence, a light, single-channel, radio-con- 
trolled fuel shutoff had to be developed. This shutoll 
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echanism was so light that it could be easily installed 
all 15-Ib. models utilizing a pair of 4'/2-lb. thrust 
This complete installation, in- 


m 
in sm 
units for propulsion. 
cluding radio, fuel tank and cutoff, and a pair of en- 
gines, is shown in Fig. 15. 

Subsequent testing with the single-channel control of 
fuel shutoff has proved so successful that the free 
launch technique has become a standard test procedure, 
in spite of the fact that it was originally a temporary 
measure to expedite the development of jet power in 
free-body models. Completely satisfactory runs can 
be made up to getaway, and in many cases short 
straightaway flights may be accomplished in complete 
safety to the model. Inasmuch as the speed and ac- 
celeration of the towing launch is limited, this free 
launch procedure allows the survey studies of small 
models to be carried out to high planing speeds and ac- 
tual getaway. When. the larger jet-propelled models 
are constructed, it is likely that free launch techniques 
will be used with, possibly, the addition of aileron con- 
trol and that the weight saved will be put into internal 


instrumentation. 


EXPERIMENTAL RESEARCH 


It has been pointed out in this paper that the intro- 
duction of the dynamically similar model and its asso- 
ciated research techniques has been responsible for 
breaking the period of passive development in the wa- 
ter-based aircraft field. This has been accomplished 
by supplementing design criteria established by the 
rule of thumb and gradual process of experience by full- 
scale trial and error, with fundamental information ob- 
tained from a direct experimental approach where all 
of the parameters affecting hydrodynamic performance 
and efficiency appear as configuration or test variables 
subject to evaluation and While 
many new design criteria have been established since 
introduction of the dynamically similar research tech- 
nique, it is not within the scope of this paper to present 
Rather, it is believed more de- 


accurate study. 


each study in detail. 
sirable to discuss at some length the several broad basic 
concepts that contributed the most toward the radical 


resurgence of water-based aircraft. 


Hull Loading 


The load that must be supported by a flying boat 
hull or seaplane float is probably the most fundamental 
variable in the design of a water-based aircraft, since it 
fixes the basic overall size of the hull. This corre- 
sponds to the selection of the basic wing area for a new 
airplane design. As you recall, the wing area required 
is fixed by the maximum lift (or load) coefficient that 
In ef- 


lect, this determination is strictly a low-speed, maxi 


can be attained and the landing speed desired. 


mum load-carrying criterion and is largely independent 
of the high-speed characteristics of the aircraft in ques 


tion. As a matter of fact, as much research effort and 
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15. The radio gear and power-plant installation for self- 
propelled radio-controlled models of jet-powered water-based 
aircraft. 


Fic. 


expense has been devoted to developing the low-speed 
load-carrying ability of wings as any other phase of 
aerodynamic investigation. The analogy between 
basic hull size and wing area selection is interesting, 
inasmuch as both are a measure of the maximum load- 
carrying ability at the lowest speed where sustention is 
by dynamic forces. In the case of the hull, this speed 
occurs in the transition region between displacement 
and planing and is known as the hump. Just as aero- 
dynamic stability and control reach a critical value near 
the stall and must be considered in the final selection 
of aerodynamic configuration, so do hydrodynamic sta- 
bility and control reach their most critical values in the 
region of the hump. In addition, the water-based situ- 
ation is further complicated by the generation of gravi- 
tational wave systems that reach their peak intensity 
in this region. The bow wave and spray produced often 
limits the loading to a value much lower than theoret- 
ically possible because of physical inundation of power 
plants and cockpit enclosures. It can be seen that the 
major problem facing the seaplane designer in his ef- 
fort to make radical improvements in seaplane efh- 
ciency was to increase hull loading materially while at 
the same time to suppress the spray formation and 
maintain adequate stability. As mentioned previ- 
ously in our discussion of Froude Number, these critical 
low-speed functions lent themselves admirably to in- 
vestigation by dynamically similar models. 

In hydrodynamic analysis, the hull loading is ex 
pressed in terms of the nondimensional load coefficient 
C,, which is based on the hull beam as the characteris 
tic dimension and is equal to A/w6*, where A is the load 
on the hull in pounds, w is the density of water in 
pounds per cubic foot, and 6 is the maximum hull beam 
in feet. As the load on a hull continually varies with 
speed, due to the varying proportion of load carried by 
the wing, it is customary to compare hulls at their 
static load coefficient Cs, = Ao/wb*, where Ap is the total 
gross weight of the airplane. 

Up to and including the seaplanes in general use dur 
ing World War IT, the ratio of the length of the hull to 
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FIG. 

load-resistance ratio at hump speed. 
length-beam ratio; 25° deadrise.) 








The 21-ft. span powered dynamic research model 
Note pro- 
cedure employed for determining influence of power on spray. 


Fic. 17. 
from which the X P5Y-1 high L/d hull was developed 


the beam was usually between 5.0 and 6.0, and it was 
well established that for such proportions it was valid 
to use Cy, as a basis for loading comparison, regardless 
of size. In other words, any two seaplanes having the 
same general proportions of length and beam generally 
had the same hydrodynamic characteristics if the values 
of Cs, were equal. Over a considerable period of 
years, experience had indicated a range of load coef- 
ficients to which could be assigned certain general char- 
acteristics of operation. For example, it was estab- 
lished that, for exceptionally good seaworthiness, short 
take-off time (low resistance), and clean running, a beam 
should be selected which gave a static load coefficient of 
not over Ca, = 0.65. Likewise, it was generally ac- 
cepted that the upper limit of hull loading was around 
Cs, = 1.0. 

To select the beam for a contemplated design, the de- 
signer had a choice, therefore, between these limits of 
load coefficient. With the ratio of 
length to beam more or less constant, it was established 


the condition of 


by experience that, increasing the load coefficient, (a) 
increased the trim and resistance at hump speed, (b) 
increased the height and intensity of the spray, (c) 
decreased the range of hydrodynamic stability during 


RONAUT 


16. Influence of length-beam ratio on the hull load and 
(Hydrodynamic effect of 
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take-off, (d) had little effect on landing stability, (¢ 
decreased low-speed maneuverability and control, ang 
(f) reduced air drag by reducing the size of the hull 
relation to the load. Item (b) or (e) usually deter. 
mined the maximum load that was practical for the de 
sign in question. 


Length-Beam Ratio 


As discussed in the introduction to this paper, th 
first promising approach to the problem of loading was 
revealed in the reports covering the experiments cop 
ducted by the Germans on high length-beam ratip 
hulls.” * In aerodynamics, one of the fundamental 
measures of efficiency of a wing is the LD or ratio oj 
lift to drag. In hydrodynamics, this applies to a hull 
as well and is the ratio of load supported to the resist. 
ance, or A/R. A cross-plot of the German data, where 
the A/R at the hump, for various values of length-beay 
ratio, was plotted against hull load coefficient, showed 
the remarkable trend reproduced in Fig. 16. On this 
curve has been plotted the maximum spray limit con. 
sidered acceptable by the Germans, and it is interesting 
to note the great area available for improvement in 
loading and A/R over the upper limit of conventional 
seaplane practice represented by the Lb = 6.0 curve. 
It appeared from study of these data that increasing 
length-beam ratio materially reduced the energy going 
into the parasitic gravitational wave-making system 
This, in turn, reduced spray and resistance, both fac- 
tors being conducive to higher allowable loadings. 

As these data dealt primarily with the effect oj 
length-beam ratio on resistance, with the inference that 
spray reduction would also be realized, it was apparent 
that many more data, particularly regarding quantita- 
tive results of spray and stability, were required before 
correlated design parameters could be established. A 


broad research program consisting of a family of 21 dy- [ 


namically similar model configurations, covering length- 
beam ratios of 6, 8, and 10, was initiated simultaneously 
by Convair, using free-body radio control, and Stevens 
Institute, employing the Stevens’ method. The radio- 
controlled model with an L,6 of 10.0 hull installed 1s 
shown in Fig. 17. These tests confirmed the fact that 
length-beam ratio has a powerful influence on the al- 
lowable magnitude of Cs, which may be used and stil 
maintain satisfactory spray. It was found also that, 
as the Ld ratio increases, the limiting value of Cs, in 
creases in direct proportion to the function L*), result- 
ing in a smaller beam and, hence, overall size of hull for 
a given load. This was a fundamental finding of im- 
mense importance, for it established the fact that hulls 
of varying length-beam ratio will have equivalent re 
sistance and spray characteristics if a load coefficient 
similar to Cy, but based on L%) instead of 5%, is held 
constant. Dr. Davidson, of Stevens Institute of Tech 
nology, first presented this relationship® and calls this 


coefficient A», which is equal to A/wL*b. Figs. 1S(at 
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18a. Influence of length-beam ratio on loading, resistance, stability, and spray while maintaining constant beam 
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Fic. 18b. Influence of length-beam ratio on loading, resistance, stability, and spray while maintaining constant plan-form area 
(Lb = constant). 
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NACA LENGTH BEAM RATIO SERIES AND AERODYNAMIC RESULTS 
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Influence of length-beam ratio on aerodynamic charac- 
teristics from the N.A.C.A. high L/d hull series. 


Fic. 19. 


summarize the results of this fundamental finding which 
first established a definite design relationship between 
the- influence of length-beam ratio and physical hull 
size. 

Fig. 18a shows the academic relationship of load, re- 
sistance, spray, and stability versus length-beam ratio 
where the beam is held constant and the length is in- 
creased by a proportional increase in station spacing 
along the keel. Here, it will be noted that marked re- 
ductions in spray and resistance are realized as the 
length increases, indicating that the higher length hulls 
can carry much more load; or, for constant loading, the 
size can be reduced while maintaining constant resist- 
ance and spray. In an effort to determine how much 
the load can be increased or the size reduced with in- 
creasing length-beam ratio, a comparison is made in 
Fig. 18b where the hull plan-form area is maintained 
constant as the L/b varies. Here again, the length is 
varied by proportional increase in hull spacing, while 
the beam is decreased at the same time to maintain 
plan-form area constant. It will be noted that the re- 
lationship is still conservative, indicating that further 
increase in load, or less size, can be tolerated for equiva- 
lent resistance and spray. 

In Fig. 1Sc we find the relationship previously men- 
tioned where the product of Lb is maintained constant. 
Resistance and spray are observed to collapse substan- 
tially to a common curve, regardless of length-beam ra- 
tio, when compared at a constant value of the Ke load 
coefficient proposed by Davidson. By studying the 
relations presented in Figs. 18a, 1Sb, and 18c, it can be 
noted that the influence of the planing bottom length 
has been progressively weighed against the effect of the 
beam. It is also apparent that, in the case of constant 
plan form, increasing length-beam ratio would main- 
tain a hull of constant volume and decreasing frontal 
area, whereas, using the L*) relationship, which has now 
been completely substantiated by all research facilities, 
a progressive decrease in hull volume, in addition to a 
rapid decrease in beam required, is made _ possible. 
This latter characteristic is of considerable importance 
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from an aerodynamic standpoint, as indicated jp Fig 
19, where it is noted that there is a 20 per cent reductjoy 
in hull drag in going from a length-beam ratio of 6,0 
12.0. 

To assist the designer in the selection of a beam for, 
new design and to establish clearly the interrelationship 
of loading, seaworthiness, and size, the selection chart of 
Fig. 20 has been prepared. This chart gives the rela. 
tionship of static load coefficient (C,,) to length-beay 
ratio for equivalent hulls based on the L*b or Ky stand. 
ard of loading. The introduction of this relationship 
between length and beam is of extreme importance ty 
the designer, since a knowledge of its effects on the hy. 
drodynamic characteristics of a hull permits greater 
choice in size and shape of hull to arrive at optimum 4. 
mensions to satisfy a particular design specification, 
Also, it permits data obtained on any one hull of a fam. 
ily to be expanded to any other hull of which it is, 
parent. It will be noted in Fig. 20 that certain values 
of Kz have been assigned their characteristic influence 
on spray height; this characteristic, we have noted, js 
the fundamental loading design criterion. In the ini- 
tial stages of a design, a value of Kz = 0.018 should be 
assumed and a value of Ky = 0.022 should never be ex- 
ceeded, since the design loading of a new seaplane 
invariably increases somewhat as the design pro- 
gresses. 

It should be pointed out that, in addition to obtaining 
a much smaller and cleaner hull aerodynamically, in. 
creasing L/b ratio at constant A» will materially reduce 
the hull weight because of reduction in hull volume and 
beam. Structural studies” * indicate that this weight 
trend continues up to a L/d ratio of approximately 15.( 
or slightly higher. With further increase in L/d ratio, 
the hull bending moments become increasingly higher 
because of the added length of hull, and this trend be- 
gins to reverse. Just where this reversal occurs and 
how rapidly it will build up depend upon many func 
tions related to the specific design. It is interesting to 
note that increasing length-beam ratio from 6.0 to 12.0 
has resulted in reductions as high as 8 per cent in struc. 
tural weight of the hull for equal performance.* 


Forebody Proportion 


Throughout this discussion of the parameters affect 
ing hull loading, the conventional standard of using the 
overall length of the planing bottom, which includes 
both forebody and afterbody, is employed to deter- 
mine the L/b ratio of the hull. It is important to note, 
however, that the forebody is the principal load-carry- 
ing portion of the hull and that in many respects it 
bears the most direct analogy to the wing of an air 
plane, being responsive to unit loading, proportion of 
aspect ratio, pressure distribution, and section. Being 
a surface that derives its lift from operation on the sur- 
face of a medium, it depends upon those dynamic 
forces that normally act on the lower surface of a wing 





at a 
tion 
spit 
fact 
tim 
nom 
wat 
erg) 
bule 
con 
flue 
cha 
whe 
reas 
bod 
ove 
7 
the 
mor 
Pat 
loa 
bea 
ing 


+ 


+ 








— << —_ 


= 


aT 


d in Fig 
reduction 
of 6.0 to 


‘a for g 
tionship 
| Chart of 
the rela. 
th-beam 
2 Stand. 
tionship 
tance to 
| the hy- 

greater 
mum dj. 
fication, 
f a fam. 
1 it isa 
1 values 
ifluence 
ted, js 
the ini- 
ould be 
r be ex. 
eaplane 
Nn pro- 


taining 
lly, in- 
reduce 
ne and 
weight 
ly 15.0 
) ratio, 
higher 
nd be- 
rs and 
r func- 
‘ing to 
o 12.0 


struc- 


iffect- 
ig the 
“ludes 
leter- 
note, 
arry- 
cts it 
1 ait- 
on of 
3eing 
» sur- 
amic 


wing 








HIGH-SPEED 


at a finite angle of attack rather than the fluid circula- 
tion resulting from total submergence in a fluid. In 
spite of the loss of lift due to absence of circulation, the 
fact that it is operating in a fluid of approximately 800 
times the density of air results in extremely great lift for 
nominal area. As a consequence of the high density of 
r and resulting high unit pressures, release in en- 


wate 
the boundaries of the surface creates violent tur- 


ergy at 
bulence in the form of waves and spray; this must be 
controlled. As spray is the principal characteristic in- 
fuencing the load that may be carried and since this 
characteristic is fundamentally a low-speed function, 
where the forebody wetted length is a maximum, it is 
reasonable to assume that it is the length of the fore- 
body that furnishes the predominate influence in the 
overall length-beam ratio of the hull. 

This assumption was investigated by Parkinson,’ of 
the N.A.C.A., by correlating all existing full-scale and 
model information with respect to forebody fineness. 
Parkinson’s analysis indicates that, when the static 
load coefficient, C,4,, is plotted against forebody length- 
beam ratio, L;/b, on a logarithmic scale, seaplanes hav- 
ing similar spray characteristics lie along straight lines 
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having equations of the form Ca, = k(L,/b)*, where k 
has the following values: 

0.0975 
0.0825 
0.0675 
0.0525 


Lxcessive spray.......... ‘ , 
Heavy, but acceptable for overload... 
Satisfactory for normal operation..... 
Extremely light spray..... 

It will be noted that Parkinson's k is of the same or- 
der as Davidson's Ke, shown in Fig. 20, indicating that 
the factor L*b is the correct fundamental characteristic 
for comparing equivalent hulls. 

Most successful seaplane designs have incorporated a 
forebody length that varies from 55 to 60 per cent of the 
total planing bottom. Variations greater than this 
range have not proved entirely satisfactory, particu- 
larly from a stability standpoint. Research has shown 
that, where the above proportions of forebody and aft- 
erbody are maintained, either Davidson’s A, or Park- 
inson’s k will give substantially the same design load 
criterion. There has been a strong trend in modern 
seaplane research to separate the forebody and its func- 
tions from the afterbody and to treat these basic planing 
surfaces separately. This line of attack has resulted in 
the development of the planing tail and long afterbody 
form of hull. These types have shown some promise 


HULL BEAM SELECTION CHART 
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Hull beam selection chart giving the relations between length-beam ratio and load coefficient for constant values of Ke 
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with regard to certain specific characteristics; however, 
evidence has shown that there are sufficient numbers of 
interdependent parameters to preclude the listing of 
completely independent criteria. While rigorous pro- 
cedures for such a segregation may be developed in 
time in the laboratories, it is desirable for the designer 
to think of the planing bottom as a coordinated unit and 
to use the overall A» as a load criterion. He should 
keep in mind that, when comparing the load capacity of 
hulls of abnormal afterbody proportion, the additional 
length of afterbody over and above standard propor- 
tions contributes nothing to the characteristics that 
have been attributed to length-beam ratio in this paper. 
This is justified by the fact that it is the forebody that 
carries the load and generates the wave system, while 
the afterbody is merely a trimming device and the 
longer afterbodies are invariably set at a higher after- 
body keel angle so that their effectiveness is still just the 
equivalent of the standard shorter afterbody. Serious 
inconsistencies will result if the total length of these ab- 
normal hulls is used to determine the hull’s L/b ratio for 
comparison with standard forms. 


The selection chart of Fig. 20 has been set up for the 
standard hull proportion. For hulls of abnormal after- 
body length, an equivalent length-beam ratio, L’/d, 
must be employed in all cases where length-beam ratio 
is used as a standard of comparison. The effective 
planing length, L’, is the length that would result when 
the forekody of the hull in question is assumed to be 58 
per cent of the total length, as in the expression L, = 
0.58L’, where L, is the actual forebody dimension and 
L’ is the total effective length for comparative pur- 
poses. 

All of the parameters discussed thus far in this sec- 
tion have dealt with the fundamental problem of in- 
creasing the load-carrying capacity of a hull. Earlier 
in our discussion of hull loading, we found that, for 
years, established practice had dictated a design static 
load coefficient of Cs, = 0.65 with an absolute maxi- 
mum of Ca, = 1.0, above which any further loading 
would increase resistance, spray, instability, or lack of 
maneuverability to unsatisfactory values. From our 
study of Fig. 20 we see that, through the results of dy- 
namic model research, the powerful influence of length- 
beam ratio as a basic design parameter has made avail- 
able to the designer as much as 300 per cent increase in 
this critical value with no loss of hull efficiency or sea- 
worthiness. 


Hydrodynamic Stability 


Established design criteria, prior to the introduction 
of high length-beam ratio as a fundamental parameter, 
indicated that increasing the static load coefficient be- 
yond a Cs, = 1.0 would increase trim, resistance, 
spray, and instability to unacceptable values. We 
have seen from Figs. 18 and 20 that, through the intro- 
duction of length-beam ratio as a basic design parame- 
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ter, the major problem of greatly increased hull Jogq. 
ings was impressively solved with an actual overall jm. 
provement in trim, resistance, and spray. Stability 
however, has been another story and has been purposely 
overlooked so far in this discussion because of the spe. 
cialized research programs required for its solution, In 
returning to Figs. 18, it will be noted that, while resjg. 
ance, spray, and trim could be controlled by appropri. 
ate length-beam ratio for a given increase in loading, the 
lower trim limit of stability increased rapidly, with any 
increase in loading, to values that would not permit a 
stable take-off to be made. During the early phases of 
length-beam ratio research, this persistent deterioration 
of stability promised to nullify all other gains in this de. 
velopment. Conventional hydrodynamic stability had 
been none too good, and it was desirable to improve this 
characteristic; further with 
could not be tolerated. 

To review for a moment, it will be recalled that, dur. 
ing the take-off of a seaplane, after it has passed the 
hump speed, the planing trim must remain between two 


compromise stability 


limiting curves of trim or it will encounter instability 
known as “‘porpoising.”’ 
tion in pitch and rise, and, if allowed to persist, it ‘may 


Porpoising is a cyclic oscilla- 


become divergent and result in destruction of the sea- 
plane. These limiting curves are shown in Figs. 18 and 
consist of an upper and lower limit of trim, versus speed, 
between which the aircraft will be stable. It has been 
established that the upper trim limit is dependent upon, 
and is a function of, the afterbody form and _ position. 
The lower limit depends upon forebody lines. To sat- 
isfy the stability requirements for take-off, the designer 
must provide an adequate range of stable trims to ac- 
commodate the normal trim range of the hull when op- 
erating anywhere within the design range of center-of- 
gravity loading, flap, and elevator position. If the trim- 
ming moment available between the limiting curves is 
not sufficient to accommodate all normal operation con- 
ditions, the aircraft must be placarded for limited oper- 
ation or modified to broaden the limits. It is obvious 
that the critical condition for upper limit is with maxi- 
mum aft center of gravity, flaps up, power off, and up- 
or, in other words, conditions that produce a 
Conversely, 


elevator 
bow-up moment and, hence, high trim. 
lower limit is critical for conditions that produce bow- 
down moments. Inasmuch as low trim is desirable 
from a resistance and spray standpoint and the take-off 
conditions of high power and flap deflection produce 
strong bow-down moment, the lower limit is of major 
importance during the heavily loaded take-off run. As 
we have noted, the lower limit of trim increases steadily 
to higher trims with increase in loading, which is di 
rectly opposed to the desired trend. It was found in 
the past that only moderate increases in loading raised 
the lower trim limit to such an extent that it was found 
necessary to reduce flap and limit the forward range of 
center-of-gravity position in order to take off. This, in 
turn, increased the resistance and left little up-elevator 
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control available, which further decreased the load- 
carrying ability of the aircraft. 

With the large increases in loading made available 
through the employment of high length-beam ratio, this 
trend approached intolerable proportions, which indi- 
cated that a stable take-off would be impossible. Inas- 
much as the forebody of the hull determines the posi- 
tion and form of the lower trim limit, an extensive re- 
search program was initiated in an effort to isolate the 
forebody design parameters associated with lower limit 
stability. All there was to go on at the start of these 
studies was the knowledge that for satisfactory planing 
stability the forebody contacted by the water during 
planing must be free of all longitudinal curvature. At 
speeds just beyond the hump, where lower limit por- 
poising is most likely to be encountered, this wetted 
length had been found to extend approximately one and 
one-half beams forward_of the step. As a consequence, 
it had been rigid design practice to maintain an area of 
constant deadrise for approximately one and one-half 
beams forward of the step. This region was commonly 
called the “‘forebody flat.”’ 

It will be recalled that, when the first families of 
hulls of increased length were investigated, the hulls 
were lengthened by proportionate increase in spacing 
between the transverse stations of a good low-length 
parent hull. This was an obvious method and was be- 
lieved to be sound inasmuch as the flat was rigidly 
maintained, and actually lengthened, by this procedure. 
It was noted, however, that the static trim of these long 
hulls was increasing to extremely high values, indicating 
that the buoyant power of the forebody was increasing 
at a greater rate than the afterbody. In an effort to re- 
duce this power by reducing the volume through in- 
creased deadrise forward of the step, it was found that 
the lower limit of stability also improved with reduced 
forebody effectiveness. In attempting to bring the 
lower trim limit down by cutting further back into the 
flat with increased deadrise, curvature was introduced 
and improvement ceased. 

With this new trend to work on, it gradually became 
apparent that, while longitudinal curvature must be 
eliminated from the planing bottom ahead of the step, 
this need not be accomplished by maintaining a con- 
stant transverse section but could be a warped surface 
of varying deadrise so long as the longitudinal elements 
of this surface were straight lines for one and one-half 
beams forward of the step. Upon establishing the basic 
fact that the trim at which lower limit instability is en- 
countered is dependent upon the linear rate of change of 
deadrise, or forebody warping, in the planing area, it 
became a relatively simple procedure to set up design 
criteria that would provide adequate stability for the 
various combinations of length-beam ratio and loadings 
previously established. Fig. 21 gives the linear varia- 
tion of deadrise forward of the step in degrees per beam 
of forebody length required to give adequate lower limit 


stability for various forebody lengths. The degree of 
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Fic. 21. Linear variation of deadrise forward of the step in 
degrees per beam of forebody length to maintain lower limit 
stability with increase in length-beam ratio. 


warping given in Fig. 21 extends for only one and one- 
half beam-lengths forward of the step, at which point 
the buttock lines may then curve as gradually as pos- 
sible to the deadrise required at the bow for the sup- 
pression of spray. 

It will be noted in Fig. 21 that for a forebody L,/b ra- 
tio of 3.5, which corresponds to an overall L/6 ratio of 
approximately 6.0 for a hull of conventional propor- 
tions, the forebody flat has a zero rate of warp in de- 
grees per beam, resulting in a flat area of constant dead- 
rise which conforms to the old rule-of-thumb concept. 
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However, for all L,/b ratios greater than this value, 
the gradual increase in warping, to maintain equivalent 
lower limit trim, results in no region of constant fore- 
body deadrise. These criteria have been incorporated 
into the plot of Fig. 22 which shows the forebody dead- 
rise variation for a typical family of hulls of varying 
length-beam ratio and standard proportion of forebody 
to afterbody. The dashed curve, labeled ‘‘Station G”’ 
in Fig. 22, represents the one and one-half beams of fore- 
body flat or linear warp that must be free of longitu- 
dinal curvature. It is apparent that, because of the 
long forebody flat required for the lower L/d ratios, it is 
difficult to warp the bottom to the high sharp lines re- 
quired at the bow without utilizing rapid buttock 
curvature. As the L/b ratio increases, it becomes 
increasingly easy to obtain the deadrise angles of the 
order of 70° to 80° at the bow which are desired from a 
seaworthiness and low buoyancy standpoint. 
It is this easing of the forebody lines with increasing 
length which results in the major increases in per- 
formance attributed to high L/6 ratio hulls, and the 
proper use of linear warp will assure this performance 


bow 


without loss of stability. 

It was stated earlier in this paper that no attempt 
would be made to cover the entire scope of research 
investigations since the advent of dynamic model 
testing. These detailed findings are adequately re- 
ported in applicable scientific literature and do not 


warrant repetition here. However, the discussions 
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Fic. 2¢ Local spray suppression devices including the highly 
effective spray dam utilized on high-speed jet-propelled water- 
based aircraft. 
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of increased hull loading contained herein are believed 
to be of such fundamental importance that any othe 
studies may be classified as refinement and extensioy 
to this basic concept of modern hull design. Througt 
the application of the broad design criteria summarize 
in Figs. 20 and 21, we find the water-based aircraft 
potential suddenly Whit 
detailed exploration and refinement will continue un 
abated, the foundation for a revitalized program ¢ 
development is wel 


expanded many times. 


modern high-speed seaplane 


established. 


Spray Control 

During the course of the research investigations tha 
have been discussed thus far and prior to the realiza 
tion of the full impact that the high length-beam ratip 
studies would have on the course of seaplane design, 
many other approaches to the problems confronting 
water-based aircraft development were initiated. [pn 
asmuch as excessive quantities of drenching spray 
was the principal factor restricting the loading 
conventional seaplanes, the major effort was directed 
toward the development of improved spray suppression 
devices. Even though these studies proved to be oj 
little value at the time and were eventually dropped 
in favor of full exploitation of the much more promising 
influence of length-beam ratio, one development during 
this period was found much later to have a most pro 
found influence on the future course of transonic and 
supersonic water-based aircraft. For this reason, it is 
appropriate to discuss briefly localized spray suppres 
sion. 

A conventional “‘V”’ bottom planing surface, such a 
shown in Fig. 23(A), is not satisfactory for a seaplan 
hull due to the great height to which the spray climbs 
upon leaving the chine; therefore, most conventional 
hulls have incorporated some form of transverse curva 
ture in the bottom called chine flare, as indicated in Figs 
23(B) and 23(C). The most effective control through 
the use of flare is obtained by a generous radius that 
terminates with a horizontal or slightly downwar 
direction at the chine, as shown in Fig. 23(B). If: 
radius that is too small or abrupt is used, only the 
layer of water adjacent to the hull is acted upon, am 
the inertia of the unaffected mass carries it by the loci 
curvature with little change in direction. Likewis 
if exaggerated flare is employed, as in Fig. 23(C), the ad 
jacent layer will actually rebound from the unaffecte 
mass and rise higher than if no flare at all had beet 
used. 

During the course of isolating the effect of downwar 
deflection of the chine, a simple expedient was et 
ployed which consisted of attaching a metal strip to 
hull without flare, such as shown in Fig. 23(A), and thet 
progressively deflecting this strip and recording th 
spray height. As one would expect, there was no im 
provement beyond a few degrees of deflection at whit! 
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t the reflected spray became progressively worse, 
As in the course of 


poitl 


as in the condition of Fig. 25(C). 
manly research investigations, the study was continued 
“ad absurdum'’ for the record and upon approaching 
4 nearly vertical position, as shown in Fig. 23(D), an 
that marked 

It was con 


interesting phenomenon occurred in 
reductions in spray height were noted." 
cluded from these tests that the increased effectiveness 
is due to the sharp intersection of the strip with the 
bottom, which in effect acts as a dam. This powerful 
effect is lost where the lowered chine is faired by using 
reflex flare or a fillet of any nature. Because of ob- 
vious structural limitations, this highly effective ver- 
tical spray strip, shown in Fig. 23(D), was not con- 
sidered feasible for application to conventional hulls 
and for many years remained principally of academic 
interest only. 

With the advent of jet propulsion, which has en- 
couraged the development of extremely low, blended- 
hull type configurations for seaplanes, the vertical 
spray strip (or spray dam as it is now commonly called ) 
has shown great promise in being able to control 
eflectively all spray without the necessity for any 
flare, any sharp chine, or, in fact, any discontinuity 
of any sort on the basic hull bottom. This form of 
spray dam has been extensively developed by the Con- 
solidated Vultee Hydrodynamics Laboratory, and in a 
broad sense this device, as illustrated in Fig. 23(E), has 
made the high-speed jet-propelled water-based aircraft 
feasible. Because of the inherent unbroken form of 
the aerodynamically smooth blended-hull configura- 
tion, the spray dam feature may be maintained struc- 
turally by placing the structural depth to the dam on 
the outer face where it can be adequately tied into the 
continuing hull structure. By rotating the spray 
dam wedge about the upper inboard edge, it may be 
retracted flush in flight, leaving an unbroken blended 
wing and hull contour. The spray dam obtains its 
effectiveness by violently agitating the main spray 
blister and thoroughly mixing it with air. This aerated 
mass is deflected downward with great force in a high- 
velocity jet. Because of the high content of entrained 
air, this water mass penetrates the free water surface 
with little or no rebound or reflection, and the high 
velocity curtain generated effectively retains the mass 
of water not directly contacted by the dam. 

An impressive example of the spray dam in operation 
is shown in the photograph of Fig. 24, which was 
taken at the speed for maximum spray. Here it will 
be noted that the bow wave is completely suppressed, 
and at no point does the spray reach a height greater 
than the dam itself. At a comparable loading, the 
bow wave on a conventional hull, utilizing a sharp 
chine and flare, would be going completely over the 
In combination with the 
the 


the spray dam principle on a com- 


wing in the example shown. 
powerful high 
application of 
pletely smooth aerodynamic form of high critical Mach 


influence of length-beam ratio, 
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The forebody of a supersonic water-based aircraft 


24 
showing the action of the retractable spray dams at the speed for 
maximum spray and the relationship of the jet air intakes 


Number has made possible jet-propelled supersonic 
water-based aircraft incorporating hydrodynamic sea- 
worthiness and stability comparable to the best pro- 
peller-driven seaplanes of conventional form. Security 
prevents a detailed analysis of the parameters and 
factors uniquely associated with the transonic and 
supersonic water-based problem. However, it suf- 
fices to say that recent hydrodynamic research, based 
upon the fundamental principles discussed in_ this 
paper, has placed the water-based aircraft high on the 


priority list of National Defense once again. 


CONCLUSIONS 


Extensive refinement of the dynamically similar 
model and its associated research techniques during the 
years since 1938, when this new concept was first in- 
troduced in this country, has established the science 
of seaplane hydrodynamics on a par with contemporary 
fields of aeronautical research and development. 
Extensive experience has indicated that the funda- 
mental laws of dimensional similitude can be matched 
or compensated to the extent that accurate deter- 
mination of the complex integrals of unrestrained 
motion, in discontinuous media, are possible with an 
economical and rapid experimental procedure. 

One of the greatest contributions of the dynamic 
research program to the advancement of water-based 
aircraft design has the establishment of hull 
length-beam ratio as a fundamental design parameter. 
Through the application of this basic criterion and its 


dependent variables, the hull loading of seaplanes has 


been 
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been increased over 300 per cent with an accompanying 
improvement in performance, seaworthiness, and 
stability. Design charts have been prepared which 
provide the designer considerable latitude in selection 
of optimum characteristics with a high measure of 
confidence. 

Through a new blended wing-hull design approach, 
coupled with the development of an extremely effective 
spray suppression device known as a “spray dam,”’ it has 
been possible to water-base idealized high Mach 
Number aerodynamic forms without compromising 
hydrodynamic performance, seaworthiness, or stability. 
By utilizing jet power, it has been demonstrated that 
water-based aircraft can be developed into efficient 
and dependable transonic and supersonic aircraft. 

It appears reasonably certain that established re- 
search techniques, such as discussed in this paper, 
can keep the fundamental design of water-based air- 
craft abreast of current aerodynamic development. 
It is just as reasonable to assume that, without this 
unique development, water-based aircraft today would 
be relegated to nothing more than auxiliary craft, 
used in those few instances where operation on water 
became mandatory for a particular specialized mis- 
sion. 

Many of the principles and new design criteria dis- 
cussed in this paper appear for the first time in the 
Convair XP5Y-1. This high-speed seaplane, con- 
structed for the U.S. Navy, not only embodies excep- 
tional hydrodynamic performance but also holds the 
distinction of being the first multiengined aircraft in 
the United States to be powered with propeller turbines. 
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The technical mating of these two outstanding «. 
velopments is indicative of the new role of water-bage 
aircraft. 
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Application of Frequency-Response 
Analysis to Aircraft-Autopilot Stability 


CHARLES L. SEACORD, JR.* 
Georgia Institute of Technology 


SUMMARY 


A description of the frequency-response method for analyzing 
an aircraft-autopilot combination is given, and examples are used 
to demonstrate various aspects of the method and the criteria 
for stability. In particular, one numerical example illustrates 
the method of calculating the aircraft response curves. The 
procedure for determining the autopilot response curves is out- 
lined. The use of Nyquist diagrams for the analysis of fre- 
quency-response data is‘discussed, and Nyquist diagrams for 
the example are correlated with the separate frequency-response 
curves. The complete longitudinal and lateral equations of 
motion for an aircraft are presented, and, in addition, certain 
simplified equations for special cases are included. Sources of 
information concerning the aerodynamic derivatives are given 


in tabular form. 


INTRODUCTION 


oo RECENT PAPERS, among them references 1 
and 2, have been written concerning the use of a 
combination graphical-analytical method for deter- 
mining the stability of a given autopilot-aircraft com- 
bination. This frequency-response method, as it is 
called, is based upon principles much used in electrical 
engineering, but only during the last 3 years have these 
principles widely been applied to aircraft stability prob- 
lems. The method makes it possible to use measured 
autopilot characteristics, along with calculated aircraft 
characteristics, to determine the stability of the com- 
bination. It is the possibility of using actual test data 
on autopilots, thus eliminating the use of complicated 
and inaccurate equations for representing the autopilot, 
which constitutes one advantage of the frequency- 
response method over the direct analytical solution 
of the equations of motion. Another advantage is 
that various hypothetical autopilot systems may be 
studied by its use with less labor than with the com- 
plete analytical solution. 

The ever increasing use of autopilots makes it de- 
sirable for the aeronautical engineer to understand the 
methods used in autopilot work and the limitations of 
such methods. The object of this paper is to present, 
in an elementary manner and with as little specialized 
mathematics as possible, the procedure to be followed 
in the solution of a problem in autopilot analysis. 

The basic equations of motion of an aircraft are pre- 
sented, along with some simplified forms, and an ex- 

Received August 15, 1949. 
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ample of the computation of the response curves is given 
in detail, including references for the aerodynamic 
derivatives used. The method for obtaining auto- 
pilot response curves is given, and the use of both the 
aircraft and autopilot response curves in determining 
the stability of the combination is illustrated. The pro- 
cedure for converting the frequency-response data 
into Nyquist diagrams and the criteria for using the 
Nyquist diagrams are also given. 


In Appendix A, an idealized case is used to illustrate 
the similarity of the results obtained through the use of 
frequency-response curves, Nyquist diagrams, and the 
simultaneous solution of the aircraft and autopilot 


equations of motion. 


SYMBOLS 


R. = amplitude of aircraft oscillation produced 
by unit amplitude oscillation of control 
surface 

R, = amplitude of aircraft oscillation required 
to produce unit control surface deflec- 
tion by autopilot 

€ = phase angle in degrees between 6 and y, 
positive when y is ahead of 6;  sub- 
scripts are a, aircraft; p, autopilot 

ri) = angle of roll, rad. 

y = angle of yaw, rad. 

B = angle of side slip, 8 = tan™! (sideslip 
velocity / V), rad. 

F = control deflection, rad.; subscripts refer 
to a, aileron; e, elevator; r, rudder 

y = angle between flight path of center of 
gravity and horizontal, rad. Nore: 
Subscript 0 indicates steady-state value; 
symbol with no subscript indicates incre- 
ment from steady state 

bo = amplitude of sinusoidal forcing control 
motion, rad. 

r = yawing angular velocity dy /dt, rad. per sec. 

p “= rolling angular velocity d@/dt, rad. per 
sec 

t = time, sec 

w = angular frequency of oscillation, rad. per 
sec. 

We = angular frequency at which eg = e,—that 
is, at which hunting oscillation may 
occur 

m = mass of aircraft, slugs 

ky = radius of gyration about longitudinal 
principal axis, ft. 

kz = radius of gyration about vertical principal 


axis, ft 
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FREQUENCY, , rad./sec. 


Fic. 1. Frequency response curves in roll for a typical aircraft 
and autopilot. Curves indicate stability. 

n = angle between principal longitudinal axis 
and flight path, positive when principal 
axis is inclined above flight path, deg. 

b = span ot aircraft, ft. 

S = reference wing area of aircraft, sq_ft. 

V = velocity along flight path, ft. per sec. 

q = dynamic pressure, lbs. per sq.ft. 

p = air density, slugs per cu.ft. 

D = differential operator d/dt 

i = imaginary unit, i = ~/—1 

Ce = yawing moment coefficient, yawing mo- 
ment /qgsb 

C1 = rolling moment coefficient, rolling mo- 
ment/qsb 

Cy = sideforce coefficient, sideforce /qs 

Cug Cig Crg = 2C,/B, OC)/dB, Cy /2B 

Cnip Srtes &S = 0C,/O0(rb/2V), 0C,/d(rb/2V), OCy/d(rb+ 
2V) 

Cups Crp, Cry = 0C,,/0(pb/2V), 0C:/d(pb/2V), OCy-/d(pb + 
2V) 

Cu Cie Crs = OCn/06, 9C1/d8, OCy/d5 

Cans Cty Crps = WCn/2(D5b/2V), OC1/2(D6b/2V), dCi 


0(Ds6b/2V) 


K = control gearing ratio of autopilot; at zero 
frequency, 6 = Ky, K¢@ 
P4(tw) = vector ratio of output displacement motion 


to input control motion for aircraft, as a 
function of angular frequency 


P (tw) = vector ratio of output control motion to 
input displacement motion for autopilot 
KG (tw) = P,(iw):P»(iw), see reference 2 
FREQUENCY-RESPONSE METHOD 


General Description 


The frequency-response method of autopilot-air- 
craft analysis involves the use of four curves as shown 
in Fig. 1. These four consist of two pairs. One pair of 
curves describes what a forcing oscillation of the con- 
trol surface (input) does to the motion of the aircraft 
(output). The other pair of curves describes how the 


autopilot moves the control surface (output) in response 
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to the motion of the airplane (input). 
tities plotted to describe each response are the amplitude 
ratio between the airplane motion and the control syr. 
face motion, R, and the phase angle, e, between the air. 
plane motion and the control surface motion. The 


The two quan. 


relationship of the body motion both in amplitude and 
phase (time relationship) to the forcing control de. 
flection depends solely upon airplane parameters. The 
relationship of the enforced control motion to the body 
motion both in amplitude and phase depends solely 
upon autopilot parameters. These response curves 
are usually determined for a range of frequencies from 
zero to about three times the natural frequency of the 
airplane. Both sets may be determined either ex. 
perimentally or by use of equations of motion. Gen. 
erally, the response of the airplane is calculated and the 
response of the autopilot is measured. 

The response of the aircraft is obtained by calculation 
of steady-state solution of the equations of motion.' 
The procedure involves the replacement of the differ- 
ential operator D by a sinusoidal frequency variable 
‘“tw.’’ This allows the response of the airplane per 
unit disturbance to be plotted as a function of the 
frequency of the periodic disturbance. 

The response data for the autopilot are obtained by 
oscillating the autopilot at various frequencies and 
amplitudes and measuring the control or servo motion 
produced. In those cases where the control motion is 
not a sine wave, it is necessary to use some degree of 
approximation for an equivalent sine wave. The pro- 
cedure for determining the equivalent sine wave is 
given in a separate section. 

The two sets of curves on Fig. 1 were calculated and 
measured for a typical air frame and autopilot. They 
are used to determine stability as follows: The inter- 
section of the e curves is determined. This intersec- 
tion, or critical frequency, determines the frequency at 
which a neutrally stable oscillation may exist. The 
values of the R curves at the critical frequency w, de- 
termine whether the oscillation can exist. That is, if, 
at the critical frequency, the R values are equal—e., 

then there will be a neu- 
If, however, R, is greater 


the curves intersect at w, 
trally stable oscillation. 

than R,, as in this example, the oscillation will be 
damped, and, if R, is greater than R,, instability will be 
indicated. The existence of a stable or an unstable os- 
cillation is the only definite conclusion that can be drawn 
directly from the two sets of curves. By definition,* * 
the response curves are valid only for the condition 
of zero damping of the oscillation, and, therefore, if the 
R curves do not intersect at the same frequency as the 
€ curves, the « intersection no longer determines the 
frequency of oscillation. This can be stated another 
way—if the airplane-autopilot combination is either 
stable or unstable, the basic response curves will indicate 
the fact but will not define the frequency or the damp- 


ing of the oscillation. 
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Computation of Airplane Response Curves 


As mentioned in the section above, the response 
curves for the airplane may be obtained by utilizing the 
equations of motion of the airplane to calculate the air- 
plane motions in response to a sinusoidally varying 
control motion. The motion of an airplane is generally 
represented by a set of six differential equations.” 6 Be- 
cause of symmetry, there will be a negligible amount of 
coupling between the lateral (rolling, yawing, and side- 


Yawing moment equation: 


m 


i 
© th cost ¢ + bx! tin’) D* — C.. — ply + | — (ke? — ky? 
Y 21 gSb 


gSb 


Rolling moment equation: 
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slipping) motions and the longitudinal (pitching, ver- 
tical, and forward) motions, provided the oscillations 
Therefore, the equations may be divided 
In some special cases where 


are small. 
in two sets of three each. 
simplifying assumptions may be made, it is possible 
to eliminate one or two of the equations. Examples 
of these simplified systems and the complete longitu- 
dinal equations are given in Appendix B. The gen- 
eral equations for lateral motion are given below: 


b 
)cos n sin n D? — C,, OV D| o- 


Ong 2V 


b 
B = Crs, ba oa (cu, a Cups, D)s. (1) 


b m ? : . " b 
= (kx? cos? n + kz? sin? »)D? — C, o7 D| o+ Fe (kz? — kx?) cos n sin » D? — C,, D 5 y- 
, . b ; , 
Ci, 8 => Cis, + Cis, OV D 0a + Ci; 6, (2) 


Sideforce equation: 


VD W l W cos yo l 
k — -<sin Yo — Cr, D |y-| ——" + Cry sD |e + 


gS q: 21 qs 


The equations given above are written with the usual 
nondimensional coefficients, but the equations them- 
selves are dimensional because the operator D implies 
differentiation with respect to time. The dimensional 
form was used because the resulting curves must be 
compared with the autopilot test curves that are ob- 
tained as functions of time or frequency. Examples 
of completely nondimensional equations may be found 
in references 7 and 8. 

It will be noticed that there are two classes of terms in 
the equation—namely, mass and aerodynamic terms. 
The mass parameters may be obtained from a summa- 
tion of the contribution of each part of the aircraft or 
by experiment. The experimental techniques are 
described in reference 9. 

Aerodynamic Parameters.—The aerodynamic deriva- 
tives are divided into the static and the rotary deriva- 
tives. The static derivatives are those arising as a re- 
sult of some angular displacement of the airplane in 
respect to the relative wind, while the rotary deriva- 
tives are those due to an angular velocity of the air- 
plane referred to the air. 

If extensive flight tests are excepted, the most accu- 
rate method of determining all of the static and certain 
of the rotary derivatives is by means of wind-tunnel 


mVD ; aes 
| gS ice Cc vol B - Cr, 9a + 


b 
Ds, (3) 


(cy, + Cyps, 2] 


tests. A tabulation of both lateral and longitudinal 
derivatives and methods of obtaining them are given in 
Table 1. 

The lateral derivatives to be used in the example are 
given below, together with a brief note as to their 


source: 

Cn, = ().2292, from wind-tunnel force tests 

Ci, = —(.1375, from wind-tunnel force tests 

Cy, = —90.9168, from wind-tunnel force tests 

Cy, = —0.699, estimated by Cy, = Cirying + 2 X 
(1/b)? Cr, (see reference 10) 

ce = ().1719, estimated from reference 11 

Cy» = —(0.0040, estimated by CG, = -(C - 
Cy,)/8 from reference 5 

C, = —0.6429, estimated from reference 11 

Cn, = — 0.0146, from wind-tunnel force tests 

Cis, = —(0.0573, from wind-tunnel force tests 

Cips, = (, estimated from reference 64 

Cy, = 0, from wind-tunnel force tests 

Cy, = (0), estimated from reference 76 


Cy, = 0,estimated from reference 76 
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Name 


Yawing moment coefficient due to 
sideslip 


Yawing moment coefficient due to 
lag of sidewash or sidewise accel- 
eration 

Yawing moment coefficient due to 
yawing angular velocity (damp- 
ing in yaw) 

Yawing moment coefficient due to 
rolling angular velocity 

Yawing moment coefficient due to 
aileron deflection 

Yawing moment coefficient due to 
rudder deflection 

Yawing moment coefficient due to 
rudder angular velocity 

Rolling moment coefficient due to 
sideslip 


Rolling moment coefficient due to 
yawing angular velocity 

Rolling moment coefficient due to 
rolling angular velocity (damping 
in roll) 

Rolling moment coefficient due to 
aileron deflection 


Rolling moment coefficient due to 
aileron angular velocity 

Rolling moment coefficient due to 
rudder deflection 

Rolling moment coefficient due to 
rudder angular velocity 

Sideforce coefficient due to angle of 
sideslip 

Sideforce coefficient due to lag of 
sidewash or sidewise acceleration 

Sideforce coefficient due to yawing 
velocity 

Sideforce coefficient due to rolling 
velocity 

Sideforce coefficient due to rudder 
deflection 

Sideforce coefficient due to rudder 
angular velocity 

Lift coefficient due to angle of at- 
tack 

Lift coefficient due to lag of down- 
wash or vertical acceleration 

Lift coefficient due to change in ve- 
locity (Mach effect ) 

Lift coefficient due to pitching ve- 
locity 

Lift coefficient due to elevator de- 
flection 

Lift coefficient due to elevator angu- 
lar velocity 

Drag coefficient due to angle of at- 
tack 

Drag coefficient due to lag of down- 
wash or vertical acceleration 

Drag coefficient due to change in 
forward velocity 


Drag coefficient due to pitching ve 
locity 


major factor; B—minor; C 


usually 
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TABLE 1 
i a ; Relative* - References 
Test Technique Importance Subsonic Supersonic 
Force tests, slope of C, vs. y. A 5, 11, 16, 17, 18, 27, 32 3 
Note: In wind tunnel 8 = 19, 20, 22, 24, 25, 
—y 26, 28, 29, 30, 34, 
35, 36, 37 
See Cmp, B 10, 64 27, 32. 33 
Oscillation tests or curved flow A 5, 10, 11, 34, 40, 41, 32, 33. 97 
test in wind tunnel 42, 43 ‘ 
Rotating model or rotating air- A 5, 10, 11, 34, 42 27, 32, 33 
stream tests in wind tunnel 
Force tests, slope of C,, vs. 5a B 18, 71 
Force tests, slope of C, vs. 6, A 40, Fi 73, 80, 81. 89 
No measurements reported B 64 
Force tests, slope of Ci vs. 8 A 5, Li, 6, 17, 21. 27.31.90 e 
22, 23, 24, 25, 26 
28, 29, 30, 35, 37, 
39, 34 
Oscillation tests or curved flow A 5, 10, 11, 34, 42,43, 32, 33, 27 
tests in wind tunnel 76 
Rotating model or rotating air- A 5, 10, 11, 29, 30,34, 27,32.33. 46, 
stream tests in wind tunnel 39, 42, 44, 45, 76 47 
Force tests, slope of C; vs. a A 64, 71, 84, 85, 86,87 73, 75, 79, 
80, 81, 82, 
83 
No measurements reported ' B 64 
Force tests, slope of C; vs. 6, B 
No measurements reported C 64 
Force tests, slope of Cy vs. 8 A 5, 17, 22,24, 29.35, 27, 32 33 
36, 37, 38 
No measurements reported A 64 27, 32, 33 
Curved flow force tests B 34, 77 27 ,32,33.2 
Force tests in rolling flow or on B 34, 76, 77 23, 33, 72 
rolling balance 
Force tests, slope of Cy vs. 5; A 18, 21, 25, 26, 35, 70 79, 80, 81, 82 
No measurements reported B 
Force tests, slope of Cz, vs. a A 35, 39, 30, 53, 67 27, 32, 33, 5 
55, 68, 57, 78 
No measurements reported B 64 
Cross-plot of Cy, curves at various C Subsonic 32, 33 
Mach Numbers B Supersonic 
Curved flow force tests B 34, 60, 74 32, 33, 72 
Force tests, slope of Cz vs. 6, A 64, 71, 84, 85, 86,87 73, 75, 79 
80, 81, 82 
No measurements reported B 64 
Force tests A 29, 39, 59, 69, 78 27, 32, 33 
No measurements reported Cc 
Force tests, cross-plot of Cp vs. J C WU <0.7 
or velocity B M>0.7 
B Supersonic 
No measurements reported ® 


negligible 


(Continued on facing page 
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( Continued ) 
Sym- Relative* - References ~ 
‘tel Name Test Technique Importance Subsonic Supersonic 
Cys, Drag coefficient due to elevator de- Force tests, slope of Cp vs. 6, B 
: flection 
rm Drag coefficient due to elevator No measurements reported c 
%% angular deflection 
c. Pitching moment coefficient due to Force tests, slope of C,, vs. a A 8, 29, 30, 34, 35, 39, 27, 31, 32, 
4 angle of attack (static longitudi- 48, 49, 50, 51, 52, 33, 56, 72 
nal stability factor ) 53, 54, 57, 58, 59 
C. Pitching moment coefficient due to Oscillation tests and curved flow A 8, 62, 63, 74 31, 32, 33, 72 
Da “lag of downwash or vertical ac- tests 
celeration 
de Pitching moment coefficient due to Cross-plot of C,, curves at various C M < 0.7 32, 33, 72 
change in velocity (Mach effect) Mach Numbers B_ Supersonic 
B M>0.7 
é. Pitching moment coefficient due to Oscillation tests or curved flow A 8, 10, 34, 49, 60,61, 27, 31, 32, 
’ pitching angular velocity (damp- tests in wind tunnel 62, 63, 64, 65, 74 33, 72 
ing in pitch) 
Cus Pitching moment coefficient due to Force tests, slope of C,, vs. 5. A 65, 54, 64, 70 73, 75 
' elevator deflection 
Cn. Pitching moment coefficient due to No measurements reported B 64 
>, — elevator angular yelocity 
T.' Thrust coefficient Flight tests, force tests B 50 
* A—major factor; B—minor; C—usually negligible. 
Cy, = Cn = Cis, = Cups, = CY, = 0, no rudder (1) (0.01231D? + 0.006859D)y + 
r rT NNOF TN) 990n<« > 
on subject airplane (0.00003935D)¢ — 0.22938 = — 0.011466, 
os = 60.5 slugs (2) (0.000591D? + 0.006308D)¢ — (4) 
k 4.23 {t.. calculated (0.0001686D)y + 0.13758 = —0.057306, 
Z _ ae +y Le c 
- : 38 — 0.255 
ky = 0.927ft., calculated (3) 4.638Dy — 0.25506 + b 
7 = (4.638D + 0.9168)8 = 0 
0 — 
n = 0°, calculated a ‘on ; — 
b 11.5ft To utilize the frequency-response technique, it is 
= .O ft. - : 
264 Ib sa ft necessary to know the response of the airplane to a 
= 90 S. per sq.It. ° ° ° ° 
q niined _—s sinusoidal control motion of constant amplitude and 
J = 586 ft. per sec. : : bs 
5 21 sa ft at various frequencies. The response can be computed 
‘ = £1 S8q.It. . ° 
by replacing the operator, D, by a frequency variable 
Calculations.—The values of the derivatives are sub- (iw). This process, the derivation of which may be 


stituted in the equations, giving the following results: 


Eqs. (1), (2), and (3) are a system of simultaneous equations in three unknowns, ¢, 6, and y. 


steps given below. 


in @ is desired, the equations are solved for ¢ by the method of determinants. 


yields: 


This can be evaluated and simplified to an expression of the 


op 10-*( —327,300D* — 247,900D? — 6,858,000D) 


0.01231D? + 0.006859 


—0.0001686D 
4.638 


—0.001686D 
4.638D 


0.01231D? + 0.006859D 


—0.01146 
— 0.05730 
0 


0.00003935D 
0.0005914D? 
— 0.2550 


At this point the substitution D = iw is made, giving 


p 


+ 0.006308D 


meee 


0.1375 


4.638D + 0.91686, 


0.1375 


4.638D + 0.9168 


5,  10-%(3,378D® + 38,580D* + 90,440D* + 720,300D? + 23,070D 


(327,300w? — 6,858,000) + t(—247,900w) 


Zt, 


5, (3,378 — 90,440w? + 23,070) + 1(38,580w* + 720,300w) 


» ratio of ¢ to 6, as an expression in D 


found in references 3 and 4, is described in detail in the 


If the response 
In the numerical example, this 


(5) 


(6) 
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The expression in iw can be simplified and rationalized to the form 


o/b. = A + 1B 


where 


and 


1,179,000w*® — 4,779,000w* + 493,400,000w 


110,600w* — 4,320,000w* + 44,930,000w* — 1,582,000 
1,141w* + 87,740 w® — 4,724,000w* + 51,460,000w? + 52,310 


(9 


denominator of Eq. (8) 


The quantities with which it is convenient to work are €,, the phase angle between ¢ and 6,, and R,, the absolute 
value of the ratio dmar./(59)¢. These can be written ase, = tan~!(B/A) and R, = VA? + B’. By substituting 
the values of A and B into the expressions for e, and R,, the formulas for R, and ¢, as functions of w (the frequency 


of the forcing control) are found. 


curves varies considerably, depending on the particular response being studied. 


It is the graphs of these formulas that are plotted on Fig. 1. 


The form of the 
As an example of the variation 


that may occur, the yaw response curves for the same aircraft used in the preceding example are shown in Fig. 2. 


Determination of Autopilot Response Curves from 
Oscillating Table Tests 
It was mentioned that the frequency-response anal- 
ysis requires the use of two sets of curves. In the 
preceding section, the method of calculating the air- 
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Oscillating 
Table ee 





Typical test setup for measuring frequency response of 
an autopilot 
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plane response curves is given; in this section, the 
method is given for transforming autopilot test data 
into frequency response curves. 

The autopilot to be tested is placed upon the oscillat- 
ing table, as shown in Fig. 3, and position recorders are 
set up to measure the position of the table and the posi- 
tion. of the control surface. A block diagram of the 
setup is shown in Fig. 4. There may arise cases in 
which it is impossible to place the autopilot on the 
table and, at the same time, have it actuate the controls 
in the aircraft itself. In such cases it is necessary to 
measure the output of the servo without going through 
the control linkage or to build a mock-up of the control 


i] 


linkage. It can be seen in Fig. 3 that the control 
surfaces are spring-loaded in order to simulate aero- 
dynamic hinge moments because the response of a servo 
sometimes varies considerably with the load on the 
servo. The traces of table position and control posi- 
tion, when recorded on film along with suitable timing 
marks, give a record much as shown in Fig. 5. It 
should be mentioned that the use of an oscillating 
table as a source of disturbance is not universal; if the 
sensing element of the autopilot does not respond to 
angular position, velocity, or acceleration or if the 
signal from the pickoff can be supplied by simpler 
means, then some means of disturbance other than the 
oscillating table is used. 


Once the records are obtained, there are three meth- 
ods of approximating the equivalent sine wave (as de- 
fined in reference 1) which may be used. These are 


listed below: 


(1) If the control motion is a regular symmetrical 
wave (such as a square wave or sawtooth wave), it is 
possible to use the expression for the first harmonic of 
the Fourier series representation and from the dimen- 
sions of the control wave obtain the amplitude and 
phase of the harmonic. 
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CONTROL LINKAGES 
AUTOPILOT LOOP 
FEEDBACK LINK 
ROLL TIMER | CONTROL | 
L—| POSITION POSITION 
PICKUP ines PICKUP 

















CIRCUIT 
—{ recoroer |— a 


Block diagram of setup for measuring frequency response of an autopilot 





TABLE A 


. 


(1) Measure and compute the following quantities: 


Units 


Symbol Description 
to The period of the table oscillation inches 
F Speed of recorder film inches per second 
F, Calibration of control record degrees per inch 
7 Period of table oscillation second 
wo = 2xr/T Angular frequency radians per second 
A, Area under 6 curve for one-half cycle (by graphical 

integration ) square inch 

B, Area under curve of 6:(cos wt) by use of following table 


and graphical integration 


to t = (360/T) to 
in. deg. 


cos wt 


square inch per second 
6 6 cos wt 
in. in. 


t = () at start of table cycle 


(II) The quantities desired are Rp and €, where: 





Rp = 1/V A? + B? and ep = tan! (B/A) 
Aand B are dimensionless amplitude factors of the in- and out-of-phase components of the control motion, 
ie w: A; F, 
‘ 2F |amplitude of table oscillation (deg. )] 
B= By Fi 


(2) In cases where it appears that an exact answer 
is not needed (i.e., where a large stability margin exists), 
a sine wave may be faired through the control wave and 
the phase and amplitude read directly. 

(3) If method (1) is not applicable or if the results 
using method (2) indicate the need for more accurate 
data, a graphic analysis may be made as outlined in 
Table A. 

In those places where reference is made to the ‘‘area 
under a curve,” the area is that between the curve and 
the zero axis of the curve, except in those cases where 
the curve is not symmetrical about the zero axis. This 
means that, if the control curve is symmetrical, the in- 
tegration can be made without any manipulation, but, 
something must be 
The procedure 


if the curve is not symmetrical, 
done to obtain an axis of symmetry. 
used is to find a line parallel to the zero axis about which 
the curve will be symmetrical in area. This may be 
done by integrating the area between the record refer- 
ence line and the control motion and then dividing by 


«F |amplitude of table oscillation (deg. )] 
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Sample of record obtained from frequency-response test 
of an autopilot. 
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the length of one cycle. The result will be the dis- 
tance from the reference line to the 
It has been found that the accuracy and speed of all 
integrations and measurements are increased by en- 
larging the records before the computation is started. 

The above process must be carried out for records 
taken at a number of frequencies, and the values of 
R, and ¢, must be obtained plotted against w. A set 
of curves obtained from tests of a particular electric 
l and 2. 


“equal area line.” 


servo autopilot are shown in Figs. 
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Analysis of Frequency Response Data 


Plots of R and « Versus w.—In the previous sections 
there have been described the methods for obtaining 
frequency-response curves for both the aircraft and 
autopilot. As plotted on Figs. 1 and 2, the use of these 
curves, in determining the stability of the combined 
system will now be illustrated. On Fig. 1 it can be 
seen that at w, the value of R, is about three times that 
of R,, thus amply satisfying the condition for sta- 
bility. In fact, the ratio of R, to R, is such that the 
motion following a disturbance will probably be well 
damped. 

Another way of estimating the amount of damping is 
to consider the phase margin, or €, €p, at the point 
where R, = R,. A value of plus 30° (when using the 
present convention), based chiefly on experience, is 
given in reference 12 as the amount of phase margin 
which will provide satisfactory damping. The con- 
cept of phase margin can be grasped more clearly from 
the Nyquist plots, which will be presented later. 

Fig. 2 illustrates a case of indicated instability. At 
w,, the value of R, is approximately one-half that of R,, 
and the phase margin at the point where R, = R, is 
negative. If the airplane were flown in this condition, 
it would exhibit an oscillation of increasing amplitude. 
The oscillation would increase without limit, unless 
some nonlinearity in either the aircraft or autopilot 
acted to change the responses and thus limit the ampli- 
tude to some constant larger value (i.e., steady ‘‘hunt- 
ing’). 

It would be possible to make the system stable if 
values of any of the four curves could be changed suffi- 
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ciently. It so happens in the present case that the 
shape of the aircraft curves in the critical region ig jp. 
fluenced chiefly by such basic parameters as yawing 
inertia and wing loading; thus, there is no simple modi. 
fication that might be made to this particular aircraf; 
to obtain stability. It is possible, however, to change 
the autopilot curves. The means of accomplishing 
such changes will be discussed in a later section. 

Nyquist Plots: The data obtained from the above. 
mentioned frequency-response measurements and cal. 
culations may be combined into one curve and plotted 
on polar coordinates. This is one form of the type of 
plot known as a Nyquist diagram, and it yields basically 
the same type of information as can be obtained from 
the separate frequency-response curves but is, in some 
cases, simpler to apply and interpret. 

The Nyquist diagrams for the two cases of stability 
mentioned above are plotted in Figs. 6 and 7. The 
plots are made by plotting the ratio of the aircraft and 
autopilot responses expressed in polar coordinates as a 
vector. This can be done because the response curves 
represent vectors of amplitude R at a phase angle ¢ and 
the ratio becomes a vector with an amplitude of 
R,/R, at a phase angle of (eg — €,). The magnitude of 
(y/6)a/(/6), is plotted at the angle (e, — ¢,). Thus, 
both responses are combined into one set of coordinates 
and plotted as a single curve. 

It will be noticed that, for the case of instability, the 
curve passes outside the point 1, 0° and that the fre- 
quency at which the curve crosses the 0° ray is the same 
as the critical frequency found by the use of the four 
separate frequency-response curves. The reason for 
the use of the (+1) critical point rather than the more 
usual (—1) is explained in Appendix A. A rigorous 
derivation and proof of this criterion—that is, neutral 
stability indicated by the curves passing through 1, 0 
involves advanced mathematics and is beyond the scope 
of this paper. It is possible, however, to show briefly 
the correspondence between the coincident intersec- 
tions of the two pairs of frequency-response curves and 
the point 1, 0° on the Nyquist polar plot. The quan- 
tity plotted on the polar diagram is known as the trans- 
fer function and, for systems similar to those under 
(W/5)a'(6/p),. This 
uct, which is written in reference 2 as AG(iw), can be 
The significance of the 


discussion, is equal to prod- 
expressed as (R,/R,)\(€a — €). 
value of the product of the vector amplitudes, which is 
written R,/R,, becomes apparent when the condition for 
neutral stability as stated for the separate response 
curves is considered. If R, = R,, the term R,/R, be- 
comes equal to 1; if R, > R,, then R,/R, > 1. Thus, 
the phase and amplitude characteristics that cause the 
e and R intersections to occur at the same frequency ina 
plot of the separate response curves will also cause the 
polar plot of R,/R, versus (€, — €,) to pass through the 
point 1, 0 on the Nyquist plot. Similarly, if R, > 
R, at € €, = 0, the Nyquist plot will pass outside 
the point 1, 0 (see Fig. 7) and will indicate instability, 
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and if Ra < Ry, at €¢ — & = 0, the Nyquist plot will 
pass inside 1, 0 and will indicate stability (see Fig. 6). 

The quantity “phase margin,’ which was mentioned 
eviously, is shown graphically on Fig. 6 for the stable 
case. The phase margin is the angle between the 0° 
ay and a line drawn from the origin to the intersection 
of the transfer function plot with a circle of unit ampli- 


pr 


tude. 
The notation used in reference 2 and that of this 


paper are related as follows: 


Herein Reference 2 
amplitude of (¥/5)a = Ra = Pa(tw) 
amplitude of (¥/6)p = Rp = 1/P» (tw) 


The Nyquist method and its derivation are described 
completely in references 2, 12, and 13, and a number of 
uses for such diagrams are given. The chief advantage 
of the Nyquist plot is that it allows the combination of 
an indefinite number of elements into one curve. 
Handling a system with more than two elements would 
be somewhat awkward if the separate frequency re- 
sponse curves were used, but a Nyquist plot for any 
number of elements may be made merely by multiply- 
ing together the separate ratios of output to input and 
plotting the resultant quantity at the sum of the indi- 
vidual phase angles. This may be stated in the terms 
of reference 2 by saying that the product of the trans- 
fer functions of the various elements will equal the 
transfer function of the system—that is, AG(tw) system 
= P,(iw):P,(iw). This makes it possible to break up 
a complex autopilot system into its component parts 
and to isolate those parts that need adjustment. The 
method may also be used to predict the performance of 
an autopilot system composed of several existing com- 
ponents on which data are already available. 


DISCUSSION 


Applicability 


A frequency-response analysis that, like the example 
given, involves only the two transfer functions or sets 
of response curves has definite limits of applicability. 
These limits are usually set by the mechanical arrange- 
ment being studied. The method is valid, as pre- 
sented, for those cases where the autopilot is sensitive 
to one variable or where the primary motion does not, 
through aerodynamic coupling, cause other secondary 
motions that feed back into the autopilot. 


An example of the first case would be the roll stabili- 
zation system, which was analyzed above, if the auto- 
pilot was sensitive to roll only. If the autopilot sys- 
tem is modified in such a way as to make it sensitive to 
yaw, but acting through the roll control, the analysis 
grows more complex. For such a system, the method 
of analysis must be modified in such a way that both 
simultaneously. A_ block 
diagram of such a system is given in Fig. 8, and the 


responses are considered 
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Ey = Yn Wout 
E® = Ey- Pout 


Block diagram of an aircraft controlled through ailerons 
by an autopilot sensitive to both roll and yaw. 
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procedure for such an analysis is described in reference 
14, page 140. 

An example of the second case is that of an airplane 
controlled in pitch. Here, the longitudinal motion does 
not couple with any other motion, and, thus, even 
though the autopilot is sensitive to quantities other 
than pitch, the motions of the airplane will not feed 
back into the autopilot except as variations in pitch. 

A third case in which the simple analysis can be used 
is that of an aircraft with two independent autopilots 
operating separate controls when the aerodynamic cross- 
coupling of the controls is not large. The procedure 
here is to analyze the stability of one motion, say roll, 
and adjust the system in such a way that stability is 
assured. The other motion, yaw, can then be analyzed, 
assuming that the aircraft is stabilized in roll, thus 
eliminating roll as a variable. 

The critical control gearings thus calculated are only 
approximate. If it is practicable to use control gear- 
ings sufficiently different from critical in the stable di- 
rection to imply adequate damping, it will be unneces- 
sary to refine the calculation. If there is any doubt as 
to the margin of stability, the more rigorous method 
of reference 15, which considers the aerodynamic cross- 
coupling of the controls, should be used. 


Nonlinearities 


The fact that most present-day autopilots are more 
or less nonlinear requires special consideration when 
an attempt is made to apply the frequency response 
analysis to a specific design problem. Linear implies 
that, if the response of the autopilot to a sinusoidal dis- 
turbance of unit amplitude is known, the response to 
a disturbance of any amplitude may be obtained by 
multiplying the unit response by a suitable constant. 
Nonlinear implies that a change in the amplitude of the 
disturbance will.cause a change in the shape of the 
response curves, as well as in the magnitude of the 
response. Many elements in the system, such as the 
sensing element and the servos, may be inherently non- 
linear, and, in addition, the presence of limiting stops 
or dead spots in the system constitute nonlinearities. 

These factors may affect either the amplitude or the 
phase of the response, and the effects may show up at 
almost any oscillation amplitude. Thus, it is necessary 
when analyzing a new or unfamiliar system to measure 
the response of the autopilot over a range of amplitudes 
from one below the desired maximum hunting in flight 
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tion for an air-powered autopilot. 


up to an amplitude above the largest disturbance likely 
to be encountered. 

There may also exist nonlinearities in the aerody 
namic derivativesused, particularly those concerned with 
control surfaces, in computing the response curves for 
These also show up as a variation in 
If such variations do 


the aircraft. 
characteristics with amplitude. 
exist, they may be accounted for by making separate 
response calculations for a range of the derivative in 
question. 

After plotting autopilot frequency-response curves 
for various amplitudes of oscillation, 1t will usually be 
seen that the results, when used with the airplane 
curves, indicate a change from stable to unstable at 
some amplitude. It is possible for the unstable régime 
to be either above or below the critical amplitude, ' 
but, fortunately, it usually occurs in the lower range, 
thus ensuring that the airplane, if disturbed, will damp 
or build up to the critical amplitude and then hunt 
at that amplitude. If this hunting amplitude is higher 
acteristics of either the airplane or the autopilot in such 
a manner that the response curves will indicate stability 
at the desired amplitude. 

An example of the variation in autopilot response 
with amplitude of disturbance is shown in Figs. 9 and 
10. In Fig. 8 the curves are for a rate-displacement sys 
tem with air pickoffs and servos (similar to the German 
V-1 autopilot). The rise in R, for the greater ampli 
tudes of oscillation is caused by stops that limit the 
maximum travel of the controls; the change in the 
phase angle at low amplitudes is probably due to dead 
spot in the rate gyro and, to some extent, in the servos. 
In Fig. 10 the curves are for a rate-displacement sys 
tem with autosyn pickofls and an electric servo. It 
will be seen here that an increase in R, occurs at the 
higher frequency when the amplitude is changed from 
1° to 3°. This was not due to stops, which were not 
used, but to the low maximum speed of the servo. The 
low speed made it impossible for the servo to attain full 
deflection in the limited time available at high fre- 
The flat phase angle curve for the 1° ampli- 
One, the rate 


quencies. 
tude is due to possibly two factors. 
gyro has limiting stops. This makes it impossible for 
the signal from the rate gyro to increase linearly with 
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frequency. Thus, the 3° amplitude shows the infly 
ence of a large rate signal at low frequencies but not at 
high frequencies. The second factor is the maximyy 
speed of the servo motor. This acts in much the sang 
manner as the rate gyro stops and prevents the fy 
effect of the rate signal from being felt at high fy 
Thus, the | amplitud 


curve shows little rate effect at low frequencies by 


quencies and amplitudes. 


does not drop off at the higher frequencies 


The above examples are not intended to indicat 
general trends, although they are fairly representatiy 
of the difference between air and electric servo systems 
but are included to give some idea of the types of non 
linearities that may show up when actual autopilots an 


tested. 


Control Gearing 


In those cases where the response of the autopilot is 
determined at the servo, it is possible to use the method 
outlined in reference 1 and find the critical control 
gearing. In the notation used in the present report, 
this process consists of finding the intersection of 
e, and e, and then determining the factor A that will 
make R, K(R,). 
then be K° 6 for every 1° yp. 
mum ratio of 6 to Y which can be used without making 


The critical control gearing will 
This will be the maxi 
the system unstable. If R, is greater than R,, the sys 
tem is indicated to be stable, but there are other con 
siderations than stability alone. One is the amount ol 
stability or phase margin. This is usually increased by 
using a control gearing that is lower than critical. An 
other and counter consideration is that of providing ade 
quate control. The control gearing must be _ high 
enough to ensure that any possible out of trim condition 
can be overcome without excessive airplane deviation 
and that sufficient control is available to meet maneu 
vering requirements. In some cases these conditions 


cannot be met with the autopilot as originally con 
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sidered, and it will become necessary to make altera 
tions to the autopilot system which will shift the phase 
response in such a way that the intersection of ¢, and 
¢, will occur at a point that will allow the conditions on 


R, and R, to be met. 


Modifications to Improve Stability 


Either the aircraft or autopilot may be changed 
when it is necessary to improve the stability of a sys 
tem, but the basic aircraft design is usually dictated by 
performance requirements; therefore, the autopilot 
is the most convenient item to modify. An example 
of a possible modification is suggested by the curves 
on Fig. 2. It is obvious that, if the intersection of the 
ecurves could be made to occur at a higher frequency, 
the value of R, at the intersection would be low enough 
to make the system stable. The intersection can be 
moved by causing the 6 curve to have negative values in 
the region from w = 0 to |. This change may be ac 
complished by adding first derivative lead to the auto 
pilot! Common ways of obtaining lead are by use of 
arate gyroscope or by use of a rate-compensated elec 
tronic circuit. The effect of a rate gyroscope may be 
seen by comparing the e, curve of Fig. 2 with the 3 
amplitude «, curve of Fig. 10, the chief difference be 
tween the two autopilots being that the data on Fig. 10 
are for a system having a rate gyroscope. It should be 
noted that any actual rate gyro will have a resonant 
frequency beyond which the rate signal will drop off and 
that any compensating circuit will also vary in its effect 
as the frequency is increased. Use of higher deriva 
tives (e.g., acceleration signal provides second deriva 
tive lead) and of integral signals can also be made to 
extend the range of the corrections, but in each case 
the modifying element or adjustment must be designed 
to meet the particular requirements of the system at 


hand 


Autopilot Cross Coupling 


It happens that some autopilots, when subjected to a 
disturbance about one axis, will respond with control 
movements about two axes. That is, a disturbance in 
pitch may cause an aileron movement as well as the 
desired elevator movement. The aileron 
motion will, of course, cause an equally unwanted roll 


unwanted 
ing motion. Thus, it is advisable to make tests to de 
termine the existence and amount of cross-coupling. 
The amount that can be tolerated depends upon the use 
to which the aircraft is to be put and upon the general 
response characteristics. If the motion caused by cross 
coupling does not affect the performance of the air 
craft or if the period of the natural oscillation is appre 
ciably different from that caused by coupling, it is 
possible to accept the cross-coupling. Generally, how- 
ever, it is possible and advisable to eliminate such ef- 
fects entirely, inasmuch as they are due to mechanical 
imperfections rather than inherent design features. 


AUTOPILOT 
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CONCLUDING REMARKS 


The frequency response method, either by use of 
separate response curves or by use of Nyquist diagrams, 
offers a convenient means of analyzing autopilot air 
craft stability problems. The method may be used to 
determine the characteristics of an existing system or 
to design a new system and is especially useful in han 
dling certain types of nonlinear systems. The results of 
a frequency-response analysis can be used to indicate 
what sort of changes should be made to either the air 
plane or autopilot in order to obtain satisfactory sta 
bility characteristics. 

There are, however, limitations to the accuracy of 
such an analysis. These limitations are chiefly due to 
the inaccuracies in the aerodynamic data for the air 
plane, the assumptions that are made in setting up 
the airplane equations, and the existence of nonlineari 
ties in the autopilot. These errors can be minimized 
to some degree by making the airplane calculations for a 
sufficient range of aerodynamic derivatives and as 


sumptions and by testing carefully the autopilot 
over a sufficiently large range of frequencies and ampli 


tudes. 


Appendix A 


SHOWING CORRESPONDENCE BETWEEN 
AND MATHEMATICAL METHODS 


EXAMPLI 
FREQUENCY- RESPONSE 
OF STABILITY ANALYSIS 


To illustrate the validity and limitations of the fre 
quency-response method, an analysis will be made of a 
one-degree-of-freedom system. The procedure is as 
follows: The frequency-response curves for the air 
plane and autopilot will be obtained, and from them 
the stability will be determined as outlined above. To 
check the results thus obtained, the differential equa 
tions of the complete system will be set up and solved. 
The frequency and damping of the system obtained 
from the complete solution will be compared with those 
The 


system chosen is convenient for purposes of analysis and 


obtained from the frequency-response method 


is used because its mathematical solution is simple and 
well known. The new quantities used are defined as in 


Table B. 


TABLE B 


Symbol Quantity Magnitude 


Iz Moment of inertia about Z axis, slug ft 22.5 
squared 
NY Rate change of yawing moment with 105 
yawing velocity, ft.lbs. per rad. per 
SCC 
Vy Rate change of yawing moment with 
yaw angle, ft.lbs. per rad — 1,730 
Ns Rate change of yawing moment with 
control deflection, ft.lbs. per rad — 1,730 
l Control lag factor, se« Variable 
K Control gearing 0.6 
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Fic. 11. Frequency-response curves of an aircraft free in yaw 


only and an autopilot which follows the expression 6 = K(1 — 
ID)jy, K = 0.6,1 = 0.1. The curves indicate neutral stability. 


The equation 

(IzD*? — NjD — N,v = N36 (Al) 
represents the motion of an aircraft restrained in bank 
and sidewise motion but completely free in yaw. The 
equation 


5 = K (1 —ID)y (A2) 


used for the autopilot, represents an inertialess system 
so arranged that the control deflection is proportional 
to the airplane deviation and subject to a viscous lag 
varying with the rate of airplane deviation. 

First, to find the critical or natural frequency, the 
differential equation of the complete airplane-autopilot 
system issolved. Thus, taking 

(IzD? — N,D — N,)¥ = No 
and 
6 = K(1 — lD)y 
and substituting for 6 gives 
(2D? — NyD — N,v = N;K(1 — ID)y 
Combining terms gives 
{1,D? + (—Ny + N,KA)D — N, — N;Kly = 0 (A3) 
As mentioned before, the frequency response curves are 
computed for a condition of neutral stability. There- 
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fore, the values of Ny and N,/A are so chosen that the 
total damping in the system is zero. This condition 
will be met if Ny = N/K. Substituting thus in Eq. 


(A3) yields 
[[zD? + (—N, — N;K)]y = 0 (A4) 


This is an equation for simple harmonic motion with 
natural angular frequency, 


wo, = V(N, + N;K)/—Iz 


It will now be shown that the same value for w, can 
be obtained from the frequency-response method. The 
first step is to compute the response curves for the air- 
plane. This process is completely explained in the body 
of the report and in references 3 and 4 and is only out. 
lined here. 


The equation of motion for the airplane is 
(IzD? = N,D — Ny yp = N56 


Using the method previously described, the following 


expressions are obtained: 
(5) ae Se 
6 a 7 (zw? 4 N,)? + N,2w* 


ja) 


. N;Ny 
1 = 
(zw? + Ny)? + Nya 


which is of the form A + 2B, a complex number. 

The above relation describes the output motion > 
caused by a sinusoidal input of unit 6 as a vector, the 
characteristics of which depend upon the frequency 
of input, w. The properties of the vector (W/6), may 
be described in terms of the ratio of the amplitudes of 
y and 6—that is, Wmar./Smar. and the phase angle be- 
tween y and 6. 
Wmaz./Smaz.s 


The expressions are R, = V A? + B= 


N. 
R, = 5 : = amplitude ratio 
Vv (Izu? + Ny )- + Ny 2qy? 
(A5) 


ande, = tan—!(B/A). 


— Nyw 


Izu" + Ny 


€, = tan! = phase angle of y referred to 4 


(A6) 


These curves are shown in Fig. 11. 
The equation for the autopilot is treated in a similar 


manner. 
6 = K (1 —lD)y (equation of autopilot) 
ep = tan lw (phase angle of y referred to 6) (A?) 
p= (AS) 
KV14+ Po 
These curves are also shown in Fig. 11. Now the 


condition for the possibility of a hunting oscillation is 
that the phase angle curves of aircraft and autopilot 





in 


gi\ 


wh 
fer 


fre 


col 
tio 
bot 
is S 
tion 


sec 
fine 


Thi 


> pl 


€ = 


Thi: 


whi 
Syst 


lat the 
dition 
in Eq. 


(A4) 


1 with 


Wy Can 

The 
he air- 
> body 
y out. 


owing 


Nya? 


ion 
r, the 
1ency 
may 
les of 
e be- 


B? = 


(A5) 


l tod 


(A6) 


nilar 


(A7) 
(AS) 
the 


yn is 
pilot 








AIRCRAFT-AUTOPILOT 


That is, 


intersect. 
€g = € 
To determine analytically the frequency at which the 
intersection will occur, the above expressions for €, 
Eqs. (A6) and (A7), are equated, using the value for 
N,of Ny = IN, K. 


€a = & 
Niw 
tan- Vw = tan™! ( — -) 
—Izw* —- Ny 
gives, as the critical frequency, 
N,I + Ny N, + NK 

a = = (AQ) 

—IIz eae Iz 


which is identical to the result obtained from the dif- 
ferential equation of the complete system. 

The expressions for the amplitude ratios at the critical 
frequency are now compared by substituting Eq. (A9) 
in Eqs. (A8) and (A5). When this is done Eq. (A8) be- 
comes equal to Eq. (A3) and thus satisfies the condi- 
tion for neutral stability, as shown in Fig. 11. Thus, 
both the necessary conditions are met, and the system 
is shown to be neutrally stable. 

Determination of the Frequency of a Damped Oscilla- 
tion.—As mentioned in the body of the report, the inter- 
section of the e curves for neutral stability will not de- 
fine the frequency of the oscillation if the oscillation is 
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that is, if the R curves do not inter- 
The cor- 


stable or unstable 
sect at the same frequency as the e€ curves. 
rect frequency can be determined, however, if the re- 
sponse curves are computed or measured for the actual 
amount of damping present in the system. The an- 
alytical solution is as follows: 

The equations of motion are written in the same 
manner as before, but the condition of neutral stability 
is avoided by specifying that Ny, > 1N;K, thus ensur- 


ing that the motion will be convergent. Thus, solving 


(IzD? — NjD — N,v = N;K(1 — ID)y 
gives, for the damped natural frequency, 


—(— Ny, + N,Kl)? _ (Me + = 


Ond = °° 
4T,? Iz 


To check this result, the frequency-response equations 
are written in the same manner as before, but, when the 
substitution for D is made, a damping term is included. 


0.- it 
5), IgD? — ND — N, 


(instead of D = tw use D = —u + iw) where u is equal 
to the sum of the damping terms (i.e., the coefficient 
of D divided by twice the coefficient of D*) in the com- 
The equation 


Thus, in 


bined airplane and autopilot equation. 
then becomes 


(y/5)a = N3/Uz (—u + iw)? — Ny (—u + Ww) — Ny) 


This is solved to give 


2 N,V (Izu? — Igw*? + Nyu — N,)? + (27,uw + Nyw)? 
par (Izu? — Izw* + Nyu — N,)? + (21zu 0 + Ny)? 


Il 


Ea 
The autopilot equation is treated similarly and yields 


1+ lu 


(3), - 
5), K{(Q + lu)? + Le?) 


II 


Ep 


tan [(27zuw + Nyw)/(Uzu? — Izw? + Nyu — Ny)] 


tlw 


t K[(1 + lu)? + Pw?) 


tan~! [lw/(1 + lu)] 


The curves for these expressions are plotted for two values of / in Figs. 12 and 13. 
As before, the expressions for the phase angles of the airplane and autopilot are equated and solved for w,: 


& Ss 


tan-! = tan 


lw 
1+ lu 


This reduces to 


—(—N, + IKN,)? 


4 ];? 


2Tzuw + Nyw 


Tou? a I zw" + Nyu _ N, 


(% +N; 
lz 


which is the frequency of the damped oscillation, as predicted by the solution of the differential equation of the 


System. 
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The three Nyquist plots that correspond to the curves 


















































Ig80 
re shown in Figs. 11, 12, and 13 are presented in Fig. 14. 
It can be seen that each of the curves crosses the 0° ray 
20 160 ve ~ 
Ra at the critical frequency as defined by the separate 
\ phase angle curves of Figs. 11, 12, and 13. Further, 
i 140 applying the criterion that those curves that pass oyt. 
faa side (0°, 1) indicate instability and those that pass ip. 
oO “ side (0°, 1) indicate stability, it is evident that the re. 
S IS 120 \ sults obtained by both methods are identical. If the 
o . a | Nyquist plots in this paper are compared with those of 
ac i 100 —— the references, it will be noticed that the critical point 
¢ a : \ y-€p in the references is usually at the 0°, —1 point instead of 
< Iog 80 M. / __ | at the 0°, 1 point, as in Figs. 6, 7, and 14. This discrep- 
ne i a \ ce ai ancy arises because of the sign convention used in the 
o a >< “<< p autopilot equations and tests. It is usual in electrical 
x 60 aa — Paarccenenene work to represent the correction to a positive error bya 
a sf negative quantity, but, because of the N.A.C.A. sign 
= 5 40 Fal ; conventions, the correction to a plus or upward pitch er- 
J y ror is represented by a plus (down elevator is positive) 
ff Pratae. quantity. This difference in conventions causes the 
20}; SSS critical point to shift 180°. The Nyquist plots are 
f equally valid for either case, and the method of inter- 
0 0 : pretation remains the same. 
O 4 8 12 16 20 
ANGULAR VELOCITY, w, RAD/SEC Appendix B 
Fic. 12. Frequency-response curves of an aircraft free in yaw 
only and an autopilot which follows the expression 6 = A(l — 
EQUATIONS OF MOTION 


lD)y, K = 0.6,1 = 0.2. The curves indicate instability. 
Both longitudinal and lateral equations are subject to 


the following assumptions: 





: (1) All deviations from steady state are small 















































180 = that is, incremental value < |. 
rEg (2) The products of two or more incremental quan- 
20 160 N R tities may be neglected. 
Q) (3) The combined effect of two or more separate 
140 if | components is equal to the sum of the individual effects. 
a as | | (4) The changes in aerodynamic forces and moments 
a re ie © q | due to a deviation are directly proportional to the devi- 
Z ‘5 sad 20 ~{O —/AtRp 2 a ation—that is, the stability derivatives are constants. 
oe | ? ~ % | (5) The lateral and longitudinal motions are inde- 
° . [00 7 a pendent of each other. 
° uy | | 
& S Equations of Longitudinal Motion 
cl0 Zz 80 | } | 
at | | Equations in Three Degrees of Freedom.—These 
=| 2 60 equations are written in terms of the more commonly 
5 a used parameters of longitudinal motion: w, the ratio ol 
oa change in forward speed to basic forward speed; 6, the 
= 40 angle of inclination of the body axis to the horizontal; 
and a, the angle of attack. The derivatives are writtet 
20 in terms of Cp, C,, and C,, and include all those that 
a have been considered appreciable for conventional ait 
0 O a planes, plus other derivatives that may be important 
0 4 in pilotless aircraft work. It is possible that there are 
other derivatives that should be used in specific cases 





ANGULAR VELOCITY, w, RAD / SEC These can be inserted in the equations as needed 
Fic. 13. Frequency-response curves of an aircraft free in yaw provided that they are expressed in the same units as 


only and an autopilot which follows the expression 6 = A(1 , d : 
the rest of the derivatives. 


Dy, K = 0.6,1 = 0.05. The curves indicate stability. 
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AIRCRAFT-AUTOPILOT 


Drag force equation: 


mn Du = —2(Cpv + Cp,)u — (Cox + Cong D = ) x 
q- ~ 

9 * 
ai (cy, + Coy, ov D) — Cope De — T¢' sin ao X 


W 6 ' 
= "OS Yo —_ ) 
a 5 cos Yo ( a 
Lift force equation: 
mV mV a 
6 — = Dea = (2C, ob Cry) u ot 
gS gs 
" c . Cc 
(a, + Cidae 2V D) er (c, + Crp, oV D) x 
. c i W . 
5. + Cuno = DO + Tc’ cos ao a + — sin yo (0 — a) 
21 gS 
Pitching moment equation: 
a : . 2 = 
D®— = 2C,u + (Ca, + Cr, -D)x 
gSc a a 2} 


a + Cutt + Cong = Dé + (Cin + Cons, ca D) x 


y 
4 TS 
c 
where the quantities are defined as follows: 
m = mass of aircraft, slugs 
q = dynamic pressure [(p/2)V?], Ibs. per 
sq.ft. 
p = air density, slugs per cu.ft. 
5 = wing area, sq-ft. 
( = mean aerodynamic chord, ft. 
D = differential operator d/dt 
V = forward velocity, ft. per sec. 
u = ratio of incremental forward velocity to 
basic forward velocity, AV/V 
Cp = drag coefficient, drag /gS 
Cr = lift coefficient, lift /gS 
Ca = pitching moment coefficient, pitching 
moment /gSC 
Te’ = thrust coefficient, thrust /gS 
a = angle of attack, radians 
9 = angle between aircraft body axis and 
horizontal, rad. 
7 = angle between flight path of center of 
gravity and horizontal, rad. Nore: 
Subscript 0 indicates steady state 
value; symbol with no subscript 
indicates increment from steady state 
q = rate change of 6 with time, d6/dt 
6 = deflection of control surface, rad. 
y = perpendicular distance from thrust line 
to center of gravity, ft. 
Coy CL Cmo = OC p/da, 0C,/Oa, OCm/Oa 
Coa Chiw Do = 0Cp/0(Dac/2V),0C,/0(Dac 2V),0Cm = 
0(Dac 2V) 
at He = 0Cp/du, 0C,_/Ou, OCn/du 
Coy, Cia ca - = O0Cp/06, 0C,/00, OC m/O8 
Cop: Ching Cmpg = WCv/2(ge/2V),0C1/0(gc/2V),OCm/O(ge* 
2V) 
Co;, Cis, Cm = 0Cp/d5, OC /05, OC, /05 
ae = 0C,/0(Dic/2V), dC,,/0(Dic/2V) 


Ds’ ~™ ps 
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Fic. 14. Nyquist diagram for unstable, neutral, and stable air 


craft-autopilot combinations. 


Equations in Two Degrees of Freedom.—The three 
equations given before describe a motion that is usually 
the sum of two well-defined oscillatory motions. One 
of these, the phugoid, is a lightly damped, long-period 
motion in which the angle of attack remains nearly con- 
stant and the forward speed and horizontal inclination 
vary. The other is a short-period, heavily damped 
motion in which the forward speed remains constant 
and @ and a vary. It is chiefly the characteristics of 
this motion which affect the requirements that the 
autopilot must fulfill; therefore, the drag force equa- 
tion and all terms due to a variation in forward speed 
may be dropped, leaving two equations in @ and a. 
These equations are presented below: 


Lift force equation: 
mV mVDe ’ 2. 
- DO — — — (Cr, + Conn DD] X 
gs gs ™ 21 
j W 


a—-({C, : -D)0 Tc’ cos a9 a — — X 
De 2) gs 


sin yo (@ — a) = (c, + ¢,., <p) 6, 


Pitching moment equation: 


D°*6 : a c : r 
(« Mo + Cm ne oV D) an- Cg oY Do — 


qgSc >. 
pers 3 - Cc 
Tc % = (Cn + Cmps, OV p) é. 


c 


Equations of Lateral Motion 


Yawing and Sideslipping.—The complete set of three 
lateral equations was given in the body of the report. 
For certain types of aircraft, these equations may be 
simplified to two- and one-degree-of-freedom systems. 
If the aircraft is symmetrical and identical about the 
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X-—Z and Y-body axis planes (i.e., of cruciform con- 
figuration), it can usually be assumed that there is no 
coupling between the rolling and either the yawing or 
sideslipping motion. This assumption can be proved 
by a wind-tunnel test that shows the derivative C,, 
to be zero. When the rolling terms are omitted, the 
equations become: 


Yawing moment equation: 


kz*D? b 
m Z ¢ ri a Dy — Cn B = Crs, 


gSb _ atid 


a b 
(Cas, + Cus 57, D)s, 
Sideforce equation: 


( W ., . 
— — sin yo — Cy 


gS gs "2 
mVD ‘ , 7 
er oT 


b 
(Cr, aa Crys oV p) 6; 


The assumption upon which the above simplification 
was based is valid in one direction only. That is, it can 
be assumed that the side slipping and yawing does not 
cause rolling, but it cannot be assumed that a rolling 
motion will not cause yawing moments. Therefore, 
the two-equation method is valid only when the roll 
stabilization is good enough to limit the rolling veloci- 
ties to values that will not cause appreciable yawing 
moments. 

Rolling. 
ing and sideslipping allows the rolling motion of a 


Using the same assumptions as for the yaw- 


cruciform aircraft to be expressed by only one equa- 


tion: 
mkx*D? ~ oO i Q b ) 
- qSb — Ci, ope (09) — Cis, + Cis, oye 4 


ba + Cig, 0; 
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Third Symposium on Plasticity 


The Third Symposium on Plasticity will be held on September 8 and 9, 1950, under the auspices 
of the Graduate Division of Applied Mathematics at Brown University, Providence, 


American and foreign experts in the field will present papers covering their most recent experi- 
mental and analytical research on ‘‘Stress-Strain Relations’ and “Problems and Methods in the Analysis 
Institute members are invited to participate in the discussion. 


For further information, please write to Prof. D. C. Drucker, Brown University, Providence 9, R.I. 
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Torsion of Thin-Walled Closed Cylinders 
Beyond the Elastic Limit” 


G. H. HANDELMANt 
Carnegte Institute of Technology 


SUMMARY 


The torsion of thin-walled closed cylinders beyond the elastic 
limit is treated as a problem in ‘technical plasticity.” From 
assumptions on the nature of the wall thickness, the problem is 
reduced to one of a single stress and a single strain. Formulas 
for the distribution of the shear stress, the shear strain, and the 
warping function are found. An analog to the formula for the 
torque vs. angle of twist in the elastic state 1s developed. It is 
shown how the problem may be reduced to an equivalent elastic 
tube through the introduction of an artificial thickness function. 
Formulas are also given for multicell sections. Specific numeri- 
cal-graphical procedures are developed for sections of piecewise- 
constant wall thickness; two examples, single-cell and multi-cell, 


are given in detail. 


STATEMENT OF PROBLEM 


he NUMBER OF EXAMPLES that have been solved 
in the field of solids stressed beyond the elastic 
limit appears small when compared with the corre- 
sponding amount of material developed in the classical 
theory of linear elasticity. This is in large part due to 
the mathematical complexities that arise from the essen- 
tial nonlinearity of the theory. However, there appears 
to be one field of problems, not fully exploited, in 
which interesting results can be obtained without undue 
These questions might be classified as be- 
The cor- 


difficulty. 
longing to the field of ‘‘technical plasticity.”’ 
responding problems in the theory of linear elasticity 
are usually considered in texts on strength of mate- 
rials or applied elasticity. This paper is concerned 
with such a problem—the torsion of a thin-walled cyl- 
inder beyond the elastic limit. Various assumptions 
of the order of those used in the corresponding elas- 
ticity problem have been made, for the major question 
to be answered lies in the gross behavior of the section 
rather than a precise description of the stress-and-strain 
distribution. 

Torsion of thin-walled elastic cylinders has been dis- 
cussed in several books on applied elasticity’ and 
strength of materials,’ as well as in recent papers on 
applications of the theory and extensions to more com- 
plicated stress distributions.* * This paper is an at- 
tempt to extend these results to the case where the 
Stresses exceed the yield limit and the stress-strain rela- 


Received December 21, 1949 

* This paper is part of a thesis submitted in partial fulfillment 
of the requirements for the degree of Ph.D. in Applied Mathe 
matics at Brown University. 

t Assistant Professor of Mathematics 


tion is no longer linear. In particular, an analog to 
Bredt’s formula’ relating the angle of twist per unit 
length and the applied torque will be derived and applied 
to two sections, a single-cell section and a double-cell 
section. 

The cylinder in question is assumed to be a thin- 
walled, closed, prismatic tube subject to a distribution 
of end forces and couples that are statically equivalent 
to a single applied couple Q. In keeping with the semi- 
inverse method of elasticity, the exact distribution of 
these end forces can be ignored, and the stresses in the 
tube will be required only to have a resultant torque 
equal to Q. 
assumed to be free from external loads. 


The lateral surfaces of the cylinder are 
More precise 
conditions on the term “‘thin-walled”’ will be given later 
in the discussion. 

Consider a cross section of the tube taken normal to 
the axes of the cylinder and a closed curve C, the refer- 
ence line, lying in this cross section (see Fig. 1). If s 
denotes the are length along C, measured from some 
fixed point 0, then the inner wall may be defined by a 
function /;(s), which measures the distance of this wall 
from C taken along a normal to C at the point s. Simi- 
larly, the outer wall can be defined by a function f(s), 
the distance between the outer wall and C at s. The 
total thickness of the wall, ¢(s), is defined by the rela- 
tion 
t(s) = to(s) + t,(s) 


Two coordinate systems, related naturally to this 
description of the tube, will be used in the discussion to 
follow. The first coordinate system consists of the 
reference line C, the normals to C, and the parallels to 
the axis of the tube. Thus, at each point of C we have 
a basic trihedral of mutually orthogonal unit vectors, 
the unit vector f tangent to C, the unit vector 7 normal 
to C, and the ynit vector 2 parallel to the axis. These 
vectors are oriented so as to form a right-handed system 
in the order named (see Fig. 1). A point P; in the tube 
is then determined by three coordinates: 5, the dis- 
tance along C of the normal projection of P; on C; n, 
the normal distance of P; to C (positive if measured 
along the outward normal); and sz, the distance be- 
tween the cross section that contains P; and some fixed 
cross section. If the reference line is assumed to be a 
closed curve that admits a representation in terms of s 
and is twice differentiable with respect to that variable, 
then it can be shown’ that the element of are length in 
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t(s+ds,n) 
\ \ 
FIG.2 
these coordinates, d/, can be written as 
dl? = [1 + (n/R)]? ds? + dn? + dz? (1) 


where R is the radius of curvature of C. A detailed 
derivation of this formula is given in the. Appendix. 

A second coordinate system defined in the plane of 
the cross section (see Fig. 1) is particularly convenient 
for describing certain quantities that appear in the 
analysis. Let O bea fixed point within Cand P = P(s) 
be an arbitrary point on C itself. Then P c be 
described by giving r = r(s), the distance from O to 
the tangent to C at P, and p = p(s), the distance from 
O to the normal to Cat P. The functions 7 and p are 
not independent of one another but are related by two 


can 


where R is again the radius of curvature of C at p 
These relations will be found useful in discussing 
the properties of certain components of the strain 
tensor. 


ANALYSIS OF STRAIN 


The displacements arising in the tube are assumed ty 
be of the same form as those appearing in the § 
Venant solution of the elastic torsion problem 


Saint. 
that js, 
the displacements of a cross section consist of a smal] 
rotation about an axis parallel to the z- 
angle of @ radians per unit length of the cylinder and, 
displacement of each point of the cross section in the 


axis through ap 


direction of the z-axis. The vertical displacement or 
warping function w is supposed to be independent of ;: 
in other terms, w = w(s, 7). 

At a given point in the cross section, the displacement 
in the s-direction will be denoted u,, that in the n. 
direction by u,, and that in the z-direction by u,. The 
fixed point O, used in defining the functions r(s) and 
p(s), will be taken as the intersection of the axis of rota- 
tion with the plane of the cross section in question, 
Thus the components of the displacement vector will 


be given by 


us = (r+ n)0) 
Un = poz (4 
u, = w(s, n) j 


Now the physical components of the strain tensor, 


€ €,3, etc., in an arbitrary orthogonal curvilinear 


coordinate system, a, 6, y, are given by the following 


aa 


formulas :* 


OuUe o log hy Olog hk, ' 
Ene = h— —he— ty — he — u, 
Oa faye] Oy ' (3) 
hy 3 O 
és, = | iis al u,) +5 : (h2ttg) 


etc., where the element of arc length is given by 
dl? = hy”) + (dB? hy") + 


(5) can be evaluated for the case in question if the 


(da? (dy?/h3") (6) 


Eqs. 


coordinates a, 8, y are replaced by s, n, 2, respectively. 


differential equations that are derived in the Appen- : 
dix—namely, From Eqs. (1) and (6), it follows that 
hy = [1 R)|-", fAe=1, he =! 
(dr/ds) + (p/R) = 0 (2) " [1 + (#/R)] " 1, Me 
(dp/ds) — (r/R) + 1 = 0 (3) Therefore, 
] dr 
3 = 0s log [1 (n/R)| } pz 
7 1 + (n/R) ds ? on <u + Isp 
lf i alk l 
n= | R) Oz Oz) ¢ 
. 2)! + Ee Sree R) |+ize Ras ° f 
€nz = (1/2) [(Ow/On) + pb] (7) 
ly 1 ow eee || | ow | 
23 = l R) - Ze & (r+ n) 0 
_ 2\1 + ( (n/R) Os oo tren | + (n/R) f 2L1 + (n/R) Os r 
Enn = O 
é., = 0 
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TORSION OF 


Now €ss 


of the vector, €s, is one-half the change in the angle between the s 
As a consequence of the assumed displacements, all the strains except e,, and 


This has been shown to be the case for e,,, and €,, and will be verified directly for 
It is assumed, furthermore, that the cylinder wall is so thin that it will not support stresses 
that is, the shear stress 7, is zero. 


larly for the other components. 
¢,, should vanish identically. 

the strains €,, and és. 

in the direction of the vector n 


can be expected to vanish. Now 


-[ 1 * 4: l 
on i + (n/R) ds 1+ (n/R)R 


from Eq. (2). Similarly, 


(r+ n)/R l d 
om A — | —_ 
: 3 1 + (n/R) 1 + (n/R) ds 2[1 


‘ 


from Eq. (3). The condition that e,, vanishes implies 


that 


(Oow/dn) + pd = 0 


This equation can be integrated directly, since neither 


pnor 6 depends on n; hence, 


— pdn + f(s) (S) 
where f(s) is an arbitrary function of s which must be 
It may be assumed, without loss of gener- 
(S) shows that the warping 


determined. 
ality, that f(0) = 0. 
function w always varies linearly across a normal to the 
reference line; the slope depends only on the geometry 
of the cross section and the angle of twist, whereas the 
determined from the stress- 


Eq. 


constant term must be 
strain relations. 
, is one-half the shearing strain y as usually, 


Since the radius of 


Finally €,, 
employed in technical literature. 
curvature R is large compared to the wall thickness for 
most of the practical beam shapes, the shearing strain 
can be written in the following form if powers of (n/R) 


of the second degree and higher are ignored in Eq. 
(7) 
y = [1 R)| (Ow/Os) + (r + n)O (9) 


ANALYSIS OF STRESS 


It has been assumed in the discussion of the strains 
existing in the tube that only the 
Since the corresponding strain is independent of 
2, Ts: = t(s, m). Under this hypothesis, the boundary 
condition that the surface of the tube be free from stress 
The boundary condition 


shearing stress 7;, 


exists. 


will not be fulfilled exactly. 
would be satisfied if the walls were parallel to the refer- 
ence line. This assumption implies, then, that the 
walls are almost parallel to the reference line; the sur- 
face conditions will be met only to within this order of 


approximation. 


CYLINDERS BEYOND 


+ (n/R)| \ 
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measures the change in length of an infinitesimal vector in the s-direction divided by the original length 


and n directions after deformation, and simi- 


Consequently, the corresponding shear strain, 


Oz (< °) 
+ = 0 
1 + (n/R) \ds R 


|e = 


dpl 
ds$ 


62 (1 r 4 *) a 
2[1 + (n/R)| Rds} — 


jit o R) — [(r + n)/R) + 


Consider the equilibrium of an element of the cylinder 
obtained by taking two sections normal to the reference 
line C at a distance ds apart measured along C (see 
Fig. 2). If /7is the total height of the cylinder, equilib- 
rium of vertical forces requires that 


Hf r(s + ds, Hf r(s, n) dn 


where the first integral is taken over the thickness 
t(s + ds) at s + ds and the second is taken over the 
thickness ¢/(s) at s. It is customary to call the integral 


T(s) = S r(s, n)dn 


n) dn = (10) 


the shear flow because of its analogy to similar terms in 


hydrodynamics. Eq. (10) may then be rewritten as 


T(s + ds) = T(s) or T(s) = const. 
that is, the shear flow is constant throughout the whole 
cross section. This constant will be denoted simply by 
Bs 

In addition, the mean shearing stress, r 
defined by the relation 


*(s), will be 


r*(s) = [1/t(s)] S xs, n)dn = T/t(s) (11) 
Since the distribution of shearing stress across the wall 
may be expected to be well behaved, the shearing stress 
along the reference line, 7(s, 0), will approximate the 
mean shearing,stress provided the reference line is more 
or less centrally located between the outer and inner 
cylindrical surfaces. In the elastic case, the curve C 
could be determined by this requirement, the 
shear distribution across a normal will be a linear func- 
tion of the applied torque Q. Once C is determined 
for a single value of Q, it would remain the same for 
any other value. This will not be true beyond the 
elastic range, since the position of C will shift in general 
with Q because of the nonlinearity of the stress-strain 
Thus, it must be specifically assumed that, 


since 


relation. 
independent of the stress-strain law, 7*(s) = r(s, 0) = 
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r(s, n) on C to within the order of approximations con- 
sidered here. 

A closed, thin-walled, single-cell tube is statically 
determinate in the sense that the shear flow 7’ can be 
found from the applied torque Q without any reference 
to the stress-strain law. The distribution of the shear- 
ing stress along a normal to C, however, cannot be 
found without considering the strains involved. The 
resultant torque Q about a point O within the cross 
section (see Fig. 1) can be readily computed, for it is 


a *n,(s) 
Q = | / t(s, nm) [r(s) + n| dn ds (12) 
O JF ns) 


As in the previous section, r(s) is the distance from O to 


simply 


the tangent to Cats; mo(s) is the coordinate of the point 
on the inner wall corresponding to s; m,(s) is that of the 
corresponding point on the outer wall; and / is the total 
length of C. The integral in Eq. (12) can be split into 
two parts—namely, 


1 en (s) 
Q= / / t(s, n)r(s) dn ds + 
0 Ny (Ss) 
l ™n (s) 
IJ t(s,n)ndnds_ (13) 
0 No(s) 


The second integral in Eq. (13) is small compared to 
the first for most technically important sections and 
can be neglected. Since 


7 


/ t(s, n)r(s) dn = r(s)T 
J ny 


the torque may be expressed as 


l * nt, (s) "§ 
v= f [ t(s, n)r(s) dn ds = r | r(s) ds 
0 Sf (ss) 0 


(14) 


Eq. (14) is equivalent to assuming that the actual shear 
stresses may be replaced by the shear flow concentrated 
along C for the purpose of computing the torque. 
Since (1/2)r(s)ds = dA, the element of area in the 
r, s coordinate system 


Q = 2AT (15) 


where A is the total area enclosed by C. Given the 
applied torque Q, 7’can be found from Eq. (15) without 
reference to the stress-strain law. The relation be- 
tween Q and 7 is the same for both the elastic and the 
plastic range. The actual distribution of r(s, 2) cannot 
be determined from this result alone; the succeeding 
section will deal with this problem. Finally, the rela- 
tion between Q and 7 is independent of the position of 
the point O, and thus any point may be used as origin 


‘ 


within the curve C 


STRESS-STRAIN RELATIONS 


As has been pointed out in the previous sections, only 
one component of stress and one component of strain 
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(the shearing stress and corresponding shearing strajy 
need be considered. Accordingly, then, the only cha. 
acteristic of the material which is required is the shes, 
stress vs. shear strain curve. It is assumed that th 
functiont 


t(s, nm) = gly(s, n)] re) | 
has an inverse that will be denoted by 
y(s, n) = h[r(s, n)| 7 


Since the sign of 7 and y can be chosen arbitrarily, bot 
g and # must be odd functions of their arguments: oy 
sequently, g(0) = h(O) = 0. Within the elastic inj 
Eq. (16) can be written in the form 


: t(s, n) = G y(s, n) 
where the constant G is the modulus of rigidity of th 
material. Since the problem concerns only loading 
nothing will be said about the equation or conditions 
governing unloading and permanent set. 


Along the curve C, n = 0; the shear strain at points 





of Cis given by y(s, 0) and the shear stress by r(s, 0) = 


7*(s). Thus, from Eqs. (9) and (17), 


Ow(s, 0) 


a + r§? = h[r*(s)| 18 


y(s, 0) 


On the other hand, Eq. (8) states that 


df Ow(s, 0) 
ds —as 


f(s) = | h[r*(s)|ds — of r(s) ds 
0 0) 


Therefore, Eq. (S) becomes 


w(s,n) = —pén + f h [r*(s)] ds — 0 / r(s) ds 
0 e 0 


Now 7*(s) = 7'/t(s) by definition [Eq. (11)] and T = 


= h[r*(s)] — r@ 


or 


Q/2A |Eq. (15)|; consequently, the warping function | 
may be expressed in terms of the applied torque | 


and resultant angle of twist @ as 


es ) s 
w(s,n) = —pén + / a| ‘ |as — of r(s)d 
0 2At(s) 0 


] 
(19 


Once a relation is found between Q and 6, the warping 
function can be expressed solely in terms of one of thest 
parameters. For this purpose, consider Eq. (IS 


integrated along the total length of C, 


+ A more general stress-strain law of the form ; = g(y, 7/7 
where 7 and ¥ are the first derivatives of r and y with respect t 
some homogeneous parameter, might also be considered With 
reasonable assumptions on the function g, however, this equa 
tion can be integrated with respect to the parameter and rewrittet 
in the form of Eq. (16). 
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Ow(s, O) 


-_ 


The warping function must be a single-valued, periodic 
function of the total length of C; accordingly, the first 
It has been previously pointed out 
Therefore, 


ds + 0 $ ris = ¢ h [r*(s)|ds 


integral vanishes. 
that the second integral has the value 2A. 


l , Q 
2A p , laa | ” 


In the elastic case, h(t) = +/G; thus Eq. (20) reduces to 


the following expression : 


(Q/4A°G) § [ds /t(s)] 


= (20) 


§= (21) 


in the elastic range. Eq. (21) is one of the classical 


Bredt formulas.® It is natural then to introduce an 


artificial thickness functjon /’(s) defined by 
t'(s) = Q/2AGh[Q/2AXt(s)| 


(20) reduces to a generalized Bredt’s for- 


(1 - 
= (1- 


Then, Eq. 


y(s5,”) = 


- )} —né a 
R/\ ds 


from Eq. (3). Finally, 


nr ) 
R? 


( alal @ 4 a(2 es 
y(s,n) = _ — 
R}" (2An(s)) 7" R 
(25) 
In the elastic range, 
(: dL Q | — {, om 4 ~) 
VS, n) = — 1 a 
rR) Laci) > PRT Re 
(26) 


and with the introduction of the artificial thickness 


~.) 
R? 


(27) 


function, Eq. (25) becomes 
n QO n 
y\5,n) = {1 — = + Ot 2 ~ 
R/|2AGt'(s) R 


Comparison of Eq. (21) with Eq. (22), Eq. (23) with 
Eq. (24), and Eq. (26) with Eq. (27) shows that the 
angle of twist vs. applied torque relation, the warping 
function, and the shear strain for a material with a 


+ 


general stress strain law can be found by solving the 
elastic torsion problem for a tube with a thickness func- 
tion given by ?¢’(s). Unfortunately, the thickness of 
the equivalent elastic tube is a nonlinear function of the 
applied torque Q; thus, a different elastic tube must be 


studied for each value Q desired. The distribution of 


} Q 0 l 
bile 2AUs)| — S 


= no ) + A ? | ors +r + no 
, —_ — a i 7 
R)\" R | PAs) ' 
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mula—namely, 
0 = (Q/4A°G) $ [ds/t'(s)] (22) 


which is valid for a general shear stress-shear strain 
relation. 
Similarly, the warping function w(s, 2) within the 


elastic limit is given by 
Q f 
— pin 
anid 2AGS 0 


Introduction of the artificial thickness function ?¢’(s) 
reduces the general formula for the warping function 


[Eq. (19)] to 
7 QO f ds of j 
ae an (Ss) d 
aia 2AGSo t'(s) J0 a 


(24) 
The shearing strain can be evaluated by substituting 
Eq. (19) into Eq. (9), 


. 


— of r(s) ds 
J 0 


(23) 


ds 


w(s,n) = 
t(s) 


wis,n) = 


+ (r+ n)@ 


shear stress does not follow this pattern, however, for 


in the elastic case it 1s given by 


(: aT QO | 4 0G (2 n n =) 

s. = = NOW i 2 

—— R/ | 2At(s) R" R? 

whereas in the nonelastic case it is 

( a) QO | 

T(s =g<il — - 

‘5 8) = BY R) | 2At'(s) | > 
P (2 n 4 =) ion 
n a= Zs) 

R- R/S 


Finally, Eqs. (22), (24), (27), and (28) constitute a 
complete solution to the problem of finding the angle 
of twist, the warping function, and the distribution of 
shear stress and shear strain as a function of the 
applied torque. 

FOR MULTICELL 


TorQuE VERSUS ANGLE OF TWIST 


SECTIONS 


By means of the results obtained in the previous sec- 
tion, it is possible to solve the technical problem of 
finding torque vs. angle of twist curves for multicell 
sectious—that is, closed sections with several internal 
webs (see Fig. 5). The technique used by von Karman 
and Christensen’ for the elastic case will be generalized 


to stresses beyond the elastic limit. 
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TORSION OF CYLINDERS 


The reference line of the outer cell and each of the 
interior webs will be defined in the same manner as in 
the first section. Let p;,, denote the reference line 
between the 7th and kth sections and p; denote the refer- 
e of the remainder of the 7th cell not included 
in the webs. The thickness function for the 7-k web 
will be given by /;, and that of the 7th cell by ¢;. The 
shear flows, 7';, in each section will have the same posi- 
tive orientation (counterclockwise). The functions 
T, will not, in general, be continuous at the junction 
points, but the resulting stress concentrations will 
be ignored in the computation of the total applied 
The shear flow in the 


ence lin 


torque, as in the elastic case. 


2A,0 = 


where the integration is taken over the web or remainder 
of the cell designated by the subscript on the integral. 
Eqs. (30) constitute the nonlinear analogs of the 
linear, algebraic equations for the 7; found by von 
Karman and Christensen in the elastic case. 

Finally, the total applied torque can be found by 
adding the torque produced in each cell as computed 
from Eq. (15)—namely, 


(31) 


Eqs. (30) and (31) will give a complete solution of the 
problem of finding the angle of twist per unit length as 
a function of the applied torque. Normally, the com- 
putation is carried out by selecting a value of 6, finding 
the shear flows from Eqs. (30), and then finding the 


torque produced from Eq. (31). 


The next section contains an application of these 
methods. 
CONSTANT 


SECTIONS OF PIECEWISE 


THICKNESS 


APPLICATION TO 


Most technically important sections are ones in which 
the thickness function is piecewise constant. The com- 
putation problem is somewhat simplified for shapes of 


this type. Two examples, taken from the paper® of 
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i-k web will be denoted by 7;,. Finally, the area of 
the 7th cell will be given by A;. 

The proof, given in the third section, that the shear 
flow is constant holds equally well for each cell of the 
multicell section. If the positive sense of 7, is 
taken to be the same as 7), then the equilibrium condi- 
tions at the junction of the web yield 


ee ee (29) 


Each cell must have the same angle of twist 6 per unit 
length as the whole section; therefore, Eqs. (15), 
(26), and (29) may be applied to each cell, in turn, to 


give, for an n-cell section, 


=~ 7, , T, — Tiss 
ds + h — Ids 
i (S) Pi, i+] ti, i419) 


(30) 


Paul Kuhn,f will be considered. One section consists 
of a single cell; the other, of two cells. 

The single cell, shown in Fig. 4a, can be simply an- 
alyzed through a suitable modification of the scales on 
the stress-strain diagram, for in this case Eq. (20) can 
be written as 


l QO ) ( QO )] 
~ 24 E (= in CY TN 


6 = (32) 
where /; is the length of the reference line of the portion 
of thickness ¢; and /. is the length of the reference line of 
the portion of thickness 4. Eq. (32) states that the 
angle of twist can be found by adding together two 
shearing strains: one shearing strain is /;/2A times the 
shearing strain that corresponds to a shearing stress of 
magnitude Q/2At,, while the second shearing strain is 
l,/2A times the shearing strain that corresponds to a 
shearing stress of magnitude Q/2Aft.. Consider, for 
example, the shear stress vs. shear strain curve shown 
in Fig. 5. This curve has been modified by adding two 
scales to the vertical axis and two corresponding scales 
to the horizontal axis. The first vertical scale is found 
by replacing the 7 scale by a corresponding scale marked 


off in terms of Q by means of the relation 
a Q 2A 


, = Ti 


+ The first example is.simply the leading edge cell of the sec- 


ond exampie. 
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The corresponding horizontal scale, 6, scale, can be 


determined by the relation 
6; = yl, 2A 


The other pair of scales are found in exactly the same 
fashion from the equations 


T= TT: = Q 2A to, d= ¥ls 2A 


For the particular section studied, the constants have 
the values A = 392 sq. in., 4) = 51 in., 4, = 0.020 in., 
I, = 24 in., and fg = 0.051 in. The Q vs. 6 graph for 
the whole section can then be determined by picking 
a value of Q, finding this Q value on the 7 scale, and 
finding the corresponding 6, value from the stress-strain 
curve and the @ scale. Similarly, 6) is determined by 
reading off the value on the 6, scale which corresponds 
to the given Q reading on the 72 scale. Finally, 6 is 
found from Eq. (32) by adding @, and . In this 
manner, the Q vs. @ graph can be readily found; the re- 
sults of the computation are given in Table 1 and are 
shown in Fig. 6. 

The two-celled cross section shown in Fig. 4b cannot 
be handled so simply, for it reflects the essential non- 


linearity of the problem. Eq. (30) becomes, in this case, 


T td aT. —_ 7 
2A,0 = fi ( ') ds + / h ( ‘ ds 
p, ty Pi,2 ty 2 
T2 F T, — T2 
2A. = f»( ‘) ds — / i ( : ‘) ds 
P, ly pi,2 hie 


Since f, f, and tf,» are piecewise constant functions, 
these equations may be readily integrated; in terms of 
the constants shown in Fig. 4b, they become 


‘ Ty ) fi. = 
51 24 J 
‘  ( T ( 0.051 
‘oe T, 
44) : 20h : ) 
(<5) alla (—. r 
2 a 
sn ( : ) - un (4 *) 
0.030 0.051 


The computation problem is somewhat simplified if 


7840 


(33) 


1,9800 


these relations are written in the form 
y*(T2) = 2,7640 — 51h(T,/0.020) 
51h(71/0.020) = 7840 — 24h [(7, T»)/0.051] 


where 


TABLE 1 
Q (108 Ib. in.) 6 (rad. in.~!) 
0 0 
20 0.000186 
40 0.000435 
60 0.000661 
80 0.000929 
100 0.001694 
120 0.003430 
140 0.006230 
160 0.010750 
180 0.017080 
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y*(T2) = 44h(T2/0.073) + 20h(T2/0.036) + 


44h(T./0.30 


The functions 51h(7)/0.020), 24h[(7, — T>2) 0.051], and 
y*(T2) can be plotted as functions of the arguments 
T,, 7; — T2, and 7», respectively, by simple modif. 
cations of the stress-strain curve of Fig. 5 as in the 
previous example. Once these curves have been ob. 
tained, a rather simple iteration scheme can be employed 
to find 7, and 7, for a given value of 6. Let 7,"" anq 
T,‘” denote the mth approximation for 7; and 7). They 
the (x + 1) approximation can be obtained from 


51h(7\'"*? /0.020) = 7840 — | 
24h[(T" — T2")/0.051] \ (34 


+1) 9.020) 


y*(T2""*) = 2,7640 — 51n(7;" 
where the values 7\'"*", 7,°"*" can be found from the 
graphs just described. 

In practice, it has been found convenient to compute 
the 7, JT. vs. 6 curves in the following manner. The 
points are computed successively for increasing values 
of 6. For the first value of 0, the first approximation 
in Eqs. (34) for 7; — 72 was taken to be zero. 
ceeding values of 6, the first approximation taken for 
T, — 7T> was the final value of this difference obtained 
from the computation for the preceding value of 6. To 
within the order of accuracy used, the first stage of the 
iteration scheme required five steps; from this point on, 
If errors 


For suc- 


only three steps were required for each point. 
are committed in the computations, the process will 
still converge, although a little more slowly. A formal 
proof of the convergence can be given in the standard 
fashion. Finally, the torque vs. angle of twist per 
unit length curve can be found by evaluating Eq. (31 


for the section in question—that is, 


Q = 7847, + 1,980 72 
The results of the computation are given in Table 2 
and are plotted in Figs. 7 and 8. 


TABLE 2 


Angle of Twist 


per Unit 
Length 
(0 (rad. in.~' X Shear Flows Torque 
10-8) | [T; (Ib. in.~')] [7Z2 (Ib. in.~')] Q (Ib. in. X 10 
0 0 0 0 
25 43.0 77.0 186.2 
37 61.0 114.0 273.5 
50 78.0 153.0 364.1 
100 114.7 199.0 183.9 
150 129.0 227 .3 551.2 
200 140.0 250.9 606.5 
250 148.5 268.5 647.7 
300 155.8 282 .2 630.9 
350 162.0 293 . 4 707.9 
100 167.2 302.5 730.0 
$50 72.0 311.0 750.6 
500 176.3 319.0 769.8 
550 180.4 326.0 786.9 
600 184.2 333 .5 804.7 
650 187.7 340.0 820.4 
700 191.0 347 .2 837.2 
750 194.2 354.2 853.6 
800 197.0 361.5 870.2 
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Appendix 


(a) The Derivation of Eq. (1) 

The line element di, in the t, n, 2 coordinate system 
can be found by considering two points, one A with 
coordinates (s, m, 3) and the other B with coordinates 
(s + ds, n + dn, 2 + dz). The projection of B in the 
cross-section plane containing A will be denoted by 
B’ (see Fig. 9). This point has the coordinates (s + 
ds.n + dn, z). In terms of the auxiliary points shown 
in Fig. 9, 

dl? = AB” + dz? 
and 
AB” = AC* + dn? 
Let R(s) be the radius of curvature of the reference 
line at the point s. Then, by similar triangles, 
AC = [(R+ n)/R]ds = [1 + (n/R)|ds 


Therefore, the line element is given by Eq. (1)—that 
is, 


dl? = [1 + (n/R)|*ds* + dn? + dz? 


(b) The Proof of Identities (2) and (3). 

Consider a point P at s on the reference line C and 
a neighboring point on C,P’, at s + ds (see Fig. 10). 
The center of curvature of C corresponding to P will 
be denoted by A. The lines OA and OA’ are perpendic- 
ular to KP and KP’, respectively; therefore, 


OA =p, OA’ = p+dp 
AP=r, A'’P=r+d 
KP = KP’ =R, PP’ =ds 
ang. (PKP’) = ang. (AOA’) = ds/R 


Since triangle AOB is similar to triangle PAP’, 
OA/KP = AB/PP’, OA/KP = OB/KP’ (35) 
Now 


AB=AK — BK =r—R-—(r+dr—R) xX 
cos ang. (PKP’) = r — R — (r + dr — R) X 
cos (ds/R) = —dr 


since higher order terms in ds/R can be neglected. The 
first of Eqs. (35) can be applied to yield 


p/R = —dr/ds 


or Eq. (2). 
On the other hand, 
OB = OA’ — BA’ = p+ dp — 
(r + dr — R) tan (ds/R) 
= p+dp — (r + dr — R) (ds/R) 


after higher order terms in ds/R are neglected. The 
second of Eqs. (35) can thus be written as 
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+d ds 
emu = e-(r+dr—R) = 
R. R R? 
that is, 
(dp/ds) — (r/R) +1=0 (36) 


since the product drds is of higher order. Eq. (36) is 


identity (3). 
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Column Buckling in the Elastoplastic Range 


P. CICALA* 


Cordoba, Argentina 


SUMMARY 


A theory of inelastic buckling of a column of constant, H- 
shaped section is presented and an approximate procedure of 
solution indicated. The numerical application confirms qualita- 
tively the results found by Shanley on the idealized semirigid 
strut. On this last simplified scheme the effect of initial eccen- 
tricities is investigated and Southwell’s plot analyzed. Con- 
sidering the curvature of the stress-strain characteristic, it 
appears that, in accordance with experimental observations, the 
load decreases with increasing deflections after reaching a value 
that, for small initial eccentricities, may not attain the tangent- 
modulus value. 


INTRODUCTION 


HE STUDY OF COLUMN BUCKLING in the inelastic 
4g pwn has received an important contribution from 
Shanley’s paper,' which showed that the buckling 
process of a straight column starts at the tangent- 
modulus load. Shanley analyzed the problem on an 
extremely simplified scheme, considering the column 
as initially straight, supposing only a short central 
section to be deformable, and neglecting the curvature 
of the stress-strain diagram. In the present paper, 
this same problem is investigated in a more general 
form, and the three mentioned simplifications are 
separately removed. The results confirm and extend 
the conclusions reached by Shanley; 
they afford a fuller explanation of experimental ob- 


furthermore, 


servations and eliminate some doubts and objections? 
that the mentioned paper had aroused. 


ANALYSIS 


The Column of Constant Cross Section 


The section of the strut is supposed to be concen- 
trated in two equal flanges, at distance 2a apart, the 
web area being considered as negligible. Let A be the 
area of each flange and L be the strut length between 
end hinges (Fig. 1). Under the compressive load P, 
the axis of the strut, initially straight, deflects in the 
plane of symmetry parallel to the web. Quantities 
corresponding to flange 1 (on the concave side) are 
distinguished by subscript 1, while subscript 2 refers 
to the opposite flange. 


In a section whose deflection is denoted by ay, 
the flange loads are 
P, = (1/2)P(1 + 9), P, = (1/2)PQ —y) (1) 
Received November 28, 1949. 
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Let «1, € be the corresponding strains and x be the 
axial coordinate measured from the mid-section anq 
referred to the half-strut length. Variable quantitie 
are considered as functions of the parameters x and 
P,. Compressive stresses and strains are considered 
Thus, we write 


(4a/L?) (0°y/Ox?) 


as positive. 


= (€2 — €)/2a (2 


The stress-strain relations 
When P increases, the load in both flanges increases and 
the representative points lie on the stress-strain char. 
acteristic of the material (line cc). If the load on one 
flange is decreasing, the representative point departs 
from line ce along a line whose slope is /. We call 
“regression zone’ the part of the strut where OP, + 
oP <0. Section s-s (Fig. 1) is supposed to be taken 
in this zone, the corresponding state of stress being 
If the portion 


are represented in Fig, 2. 


represented by points 1,2 on Fig. 2a. 
of the characteristic here involved may be represented 
by -a straight segment of slope /,, the stresses in the 


flanges are expressed by 
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(a) Constant 
(c) Semirigid scheme. 


Scheme and notation. 
(b) Cross section. 


Fig. 1. 
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2 
| 
—~ 
me 
| 


=o,+ E,(a — &) : ym the value of y for x = 0, the solution may be ex- 
(3) - . 
or +t Ele — €) 


pressed as follows: 


Q 
| 
Ss 
~ 
N 
— 
Il 


are ¢C inates of the point X of ah . 
where Or and e, are the coordinates o | a “) —" 
From Eqs. (1), (2), and (3) we get , p Jo q(x’) sin ap(x e ex’ + 


departure. 
! h)cos apx — h (7) 
mt Aly bx J k+1 (tn + 
Mea’ Oy 5 (P — 2P,) + Py=0 (4) | 
[2 ox? 2k 2k according to symmetry with respect to section x = 0. 
We call ‘‘progression zone’ the parts of the strut 
where where OP,/O0P > 0. For these portions of the column, 
P, =Ao, k= E/E which lie near the hinges, the fundamental equations 
take the simpler form that can be deduced from Eqs. 
Putting (4) and (6) fork = 1. Therefore, in the upper portion 
p? = L’P SAEa’, q = L?P,/4AE@? | (5) of the column we have 
/ r (9 
oe Db = —_— : 
a? = (k + 1)/2k, h (k I)/(k + WS y = Csin p(l — x) (9) 
Eq. (4) may be written as follows: We denote by z the relative coordinate of the sections 
that limit the regression zone. For x = 3, Eqs. (7) 


(1/a?) (0°y/Ox?) + p*y + h(p? — g) = 0 (6) 2 
and (8) must give the same values of y and Oy/Ovw. 


where a, are constants and g is a function of x. This These conditions yield the quantities C and y,,. For 
equation is valid in the regression zone; denoting by _ the first we get 


: P ' arh “s 
Cleos p(1 — 2) cos apz — asin p(1 — 2) sin apz| = ah sin aps — g cos apx dx (9) 
7 
0 


From Eqs. (1), (8), and (9) the expression of P, may be deduced as a function of p and x, containing the quanti- 
ties z, function of p, and q, function of the variable of integration. For the limiting section x = s 


OF,/aPr = @ (10) 
P, = Pf, (11) 

According to Eqs. (5), Eq. (11) yieldst 
q = p? — Cp’ sin p(l — 2) (12) 


Hence, the following set of equations may be derived: 
g 1 ) 


(q — p*)p cos p(1 — 2) + (dg/dz) sin p(1 — z) = 0 

p? , dq : Z (13) 
ag p sin apz + ] cos apz + ahp g COS apx ax = o| 
a az 0 


which define in an integrodifferential form the relation between the quantities p, 2, and q. 


+ The argument of g may be called indifferently x or z. 


The problem may be given an approximate solution of gq at station 2; , ; can be calculated immediately, from 
by means of a step-by-step procedure based directly Eq. (12), as a function of p. Then, Eq. (14) defines 
on Eqs. (10) and (11). The expression of P:, as ob- p and, consequently, q. 


tained from Eqs. (1), (8), and (9), satisfies the condition The solution depends on one boundary condition 
>P,/d 0 for instance, the value of pp of p for which q = p? 
— (i.e., y = 0). This condition corresponds to the onset 
forx =z. Therefore, Eq. (10) reduces to of buckling; under the initial value p = po, the re- 
ression zone reaches as far as section 2 given by the 
g : 
0g/O0p = 0 (14) equation 
if this derivative is taken from Eq. (12), considering tan apo tan po (1 — 2%) = a (15) 
Z as a constant and p as an independent variable. " ie ald ; 
. f | : For po = 2/2—that is, if the buckling starts at the 
Thus, if the values Gi» Jz, .--» Gi Of g at stations 2; = 
tangent-modulus load represented by 


0, 2,..., 2; are known and if we express the integral 
that appears in Eq. (9) by Simpson’s rule, the value P, = 2x°AEa’/L? (16) 
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Fic. 2a. Stress-strain diagram. Points 1, 2 correspond to 
section s-s (Fig. 1) of the regression zone; points 1’, 2’ to section 


s’-s’ of the progression zone. Fic. 2b. Flange loads vs. strain 
by nonlinear characteristic. 
Eq. (15) yields 29 = 0; whereas, for po = m/2a 


that is, if the buckling starts at the double-modulus 
P,/a*®—the regression zone covers the entire 
strut atonce. For values of P < P,, no bent configura- 
tion of equilibrium exists. From the initial limits, the 
regression zone spreads out toward the ends of the 


load 


column, while the deflections increase. 
The most important case is that by which the 
buckling begins at P = P,. In this case, expanding 


the solution of Eqs. (13) into power series of z gives 


pb = (w/2) + [w3(k — 1)/96R] 2? +... 


g = (w?/4) + [1(k — 1)/3,840R]2° +... (17) 
The results of calculations carried out for k = 2 


are shown in Fig. 3. The dotted line b represents 
Eq. (15). Line a defines the variations of the regres- 
sion zone for the fundamental case pp = 7/2. Line 
c corresponds to z = 0.6. The solution corresponding 
to po = 1/2 is also represented by the dotted lines aj, a2 
of Fig. 4, which give the loads on the flanges as func- 
tions of the applied load at section x = 0. 


The Semirigid Column 


According to Fig. lb, the strut is supposed to be 
formed by two rigid parts of equal length L/2 and a 
short central deformable section of length 56. The 
deflection ay is supposed to be constant along the 
deformable portion; let yo be the value of y for P = 0 
(initial eccentricity). Eq. (2) modifies in this case as 
follows: 


(4a/Lb) (y — yo) = (4 — €)/2a (18) 
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while Eqs. (1) and (3) hold true. Thence, 


Py — yo) = a (P — 2P;) + : = ' Py (19) 
where 
P, = 8AEa*/Lb 
In the particular case yo = 0, P, = P,, Eq. (19 


agrees with Eq. (18) of Shanley’s paper. Neglecting 
the curvature of the stress-strain diagram in the region 
implied, the same Eq. (19) holds below the regression 
point with k = 1, thus giving 


y = wIP./(P. — P)| (20 
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Fic. 3. Length of the regression zone vs. load parameter 
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Fic. 4. Flange loads vs. applied load. Lines aj, az: constal 
section column. Lines aj’, a2’: semirigid column, yo = ' 
Lines ci, C2: semirigid column, y) = 1/400. 
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MN BUCKLING IN 


COt, Uv 


The value P,, of P, for which the reversal condition 


(namely, dP2/dP = 0) is verified, results: 
ff 
Pw = Pi _ V Yo) (21) 


Under this load, the stress a2 attains the value a, 


corresponding to point X where departure from the 
characteristic takes place (Fig. 2a). Thus, we obtain 
Pf. = P2/2F, (22) 


With this value of P,, the column deflections for any 
load P > Py» can be calculated from Eq. (19); thence, 
the flange loads can be computed according to Eqs. (1). 
The results of calculations carried out for k 2 are 
4, where the ratios P;/P, and P:/P, 
are plotted versus P/?.. and C2 give, re- 
spectively, the values of P; and P, for yo = 1/400. 
The line m, whose equation is P,/P, = P?/2P,”, repre- 
sents the locus of maximum P.. Its crossing divides 
line ¢: in two parts of different character: the ‘‘pro- 
(P < Pr»), where Eq. (20) holds; 
’ where Eq. (19) is valid. 
is established 


resumed in Fig. 
Lines ¢; 


gression branch” 


and the 
The corresponding division on line c 


“regression branch,’ 


by line m;, whose equation is 


P,/P, = (P/P,;) (P?/2P,,”) 


For any value of initial eccentricity, the regression 


branches have, as an asymptote, the vertical 
P = P, = 2kP,/(k + 1) 


For vanishing initial eccentricity, both progression 
branches tend jointly toward the straight line a (P; = 
P, = P/2), whereas the sugpeneiom branches approach, 
respectively, lines a;’ and ap’, which are obtained from 
Eqs. (1) and (19) for yo = 0, P,; = P;. For the initially 
straight column (yp = 0), from the same equations, 
b. may be obtained if P, is 
any intermediate P,. 
These lines represent the buckling process in the ideal 
maintained straight until the 
They have prac- 


other curves such as by, 
given value between P, and 
case that the column be 
tangent-modulus load is exceeded. 
tically no importance, since in no case are they ap- 
proached by equilibrium lines representing the natural 
behavior of actual columns.t 

It may be easily shown that, if the values of (y — 
o)/P are calculated using Eq. (19) and plotted against 
Y — ¥% according to Southwell’s representation, we 
obtain a curve whose slope tends asymptotically to 


es. & simpler representation is reached using the 


form of this plot as modified by Lundquist.’ In this 
case the parameters following are introduced: 
u=(y—y')/(P—P’), v=y-—y’ (23) 


where y’ is the deflection that corresponds to a load 
P’ arbitrarily chosen. Taking P’ > P,, we get from 

t In the load vs. column shortening graph, the lines corre- 
sponding to different initial buckling loads for yo = 0 do not cut 
each other; 1, the line starting at P = P, limits the 
region covered by equilibrium lines for imperfect columns. 


as in Fig. 
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Eq. (19) 


(P, — P’)u —v ‘+h (24) 


Except for the constant term /, the equation is the 
same as if the material were elastic and the critical 
load were P,. 


The Case of Nonlinear Stress-Strain Relation 


The buckling process can be easily analyzed if the 
idealization referred to in Fig. lb is adopted. The 
flange loads are plotted versus the relative shortening 
of the same fibers in Fig. 2b. For increasing strain, 
the relation is represented by line c, according to the 
properties of the material; line t is the tangent to c 
at point X, whence the representative point for flange 
2 departed along line e. Taking into account Eqs. 
(1), Eq. (18) may be written as follows: 


P,[{(Pi — P2)/(Pi + Pe) — yw] = Pi* P* (25) 
where P,* and P,* are the loads that, on the tangent 
t, correspond to the strains that actually exist in the 
flanges. For a given P;, the value of P;* may be read 


on the diagram, while P,* may be expressed by 











P* = Pf, — (P, — P.)/k (26) 
7 
P 
Vp 
fe 
] (6) 
+0°6 K=3 2 
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Fic. 5. applied load for various values of 


Flange loads vs. 
K = —(0/2E,)(dE./do). 
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Replacing these values in Eq. (25) yields a quadratic 
equation for P,. In order to perform the calculations 
in an analytical way, considering (P, — P;) as a small 
quantity, we may approximate the experimental rela- 


tion by writing 


(27) 


P,* _ P, + K{(P, oe P,,)?/ Fe) 
where 


K = —(1/2) [(¢/E,) (dE,/do)| (28) 
the values in the brackets being calculated at point XY. 
According to Osgood’s stress-strain law (Ke = o 


+ constant:o”), we obtain 


K = [(k — 1)/2k| (n — 1) (29) 
Using the stress-strain law given in R.A.S. Data 
Sheets 00.02 for an aluminum alloy having a yield 


stress of 55 ksi. and a ultimate stress of 73 ksi., at the 
stress for which /, = //2, we obtain the value K = 
2.80. 

Fig. 5 represents the values P; and P2 as functions of 
the applied load computed for k = 2, P, = P,, and 
the values of A indicated near each curve. 


CONCLUSIONS 


The investigation of the symmetrical buckling of a 
column of constant, H-shaped section leads to a rather 
complex analytical problem, even when the curvature 
of the stress-strain diagram is neglected. The equation 
implied may be solved with the desired accuracy by a 
step-by-step procedure. The analysis shows that 
the buckling can start at any load between the tangent- 
and the reduced-modulus value. At the onset of buck- 
ling, the stresses in one of the flanges begin to decrease; 
the zone where this occurs (we call it ‘‘regression zone’’) 
covers a central portion of the strut (shadowed in 
Fig. la). The more the initial buckling load exceeds 
the tangent-modulus the wider the initial 
regression zone; its amplitude successively grows with 
increasing axial load and lateral deflections. 


value, 


In the fundamental case, wherein the buckling starts 
at the tangent-modulus load, the regression zone pro- 
ceeds from the mid-section of the column, spreading 
rapidly over a wide part of the strut. Analyzing for 
this case the variation of the flange loads in the mid- 
section, the results obtained by Shanley on the idealized 
semirigid strut appear qualitatively confirmed. The 
lateral deflections and the maximum stress in the men- 
tioned region increase rapidly as the tangent-modulus 
load is exceeded. 
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Shanley’s semirigid scheme was used to investigate 
the effect of initial eccentricities. 
strain reversal occurs in one flange at a load less tha 
The increase of deflec 


In any case, thy 


the tangent-modulus value. 
tions, and, consequently, of the maximum Stress, js 
anticipated by 
supposed to vanish, the behavior of the strut approaches 
the ideal case of the straight column beginning ty 
buckle at the tangent-modulus load. Therefore, the 
possible buckling processes starting at higher loads, a 
indicated by both the rigorous and the simplified an 
alysis, appear as actually not attainable, unless ay 
artificial restraint prevents the unavoidable eccentricities 


initial imperfections. If these are 


from their natural growth. 

Furthermore, the analysis indicates that Southwell’s 
plot of lateral deflections shows the same asymptotic 
slope as if the critical load had the reduced-modulys 
value. The representation of the bucklinz process, in 
the form proposed by Lundquist for the elastic range, 
may yield a straight graph, parallel to the line that 
should be obtained if the column behaved elastically 
according to the reduced modulus. However, these 
diagrams do not have the same value here as in the 
elastic range, where they serve to detect the true critical 
load, the extraordinarily 
accurate specimens. 
they can conceal the real character of the phenomenon, 
disguising the bifurcation of the equilibrium line that 
exists at the tangent-modulus load for the straight 


construction of 
In the case under consideration, 


avoiding 


column. 

Another aspect of the problem which can be easily 
investigated on the ideal semirigid column is the effect 
of the curvature of the stress-strain diagram. When 
this factor is taken into account, the load resisted by 
the strut does not permanently increase with lateral 
deflections, as the previous analysis indicated. Even 
in the case of vanishing initial eccentricities, the 
maximum load resisted by the column is markedly 
less than the reduced-modulus value. This maximum 
load is considerably affected by the combined influence 


of previous imperfections of the column and curvature 
of the stress-strain diagram and may not attain the 


tangent-modulus value. 
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Hydrodynamic Longitudinal Stability 


GEORGE S. SCHAIRER* 


Boeing Airplane Company 


ABSTRACT 


[he various forms of hydrodynamic instability of flying boats 


are described, and the physical phenomena associated with these 
instabilities are examined. It 


tinuity in flow characteristics when the chine at the main step 


is shown that there is a discon 
becomes submerged, which changes the wetted area from a tri 
angle to a trapezoid. In cases where the afterbody or tail ex 
tensions are not in the water, this point of changing flow defines 
The lower limit is also partially de 
which the afterbody first contacts 


the lower porpoising limit 
termined by the point at 
the water. A form of upper limit porpoising which differs from 
skipping and closely resembles lower 
pointed out that hydrostabilizerst can mate 


limiting porpoising is 


described. It is 
rially change lower-limit porpoising 


INTRODUCTION 


we LONGITUDINAL STABILITY CHARACTERISTICS of 
seaplanes while riding on the water have been the 
subject of much study in recent years. This has come 
about because of numerous cases of oscillatory insta 
bility during take-off and landing in recent seaplanes. 
For many reasons these instabilities have become more 
important in recent seaplanes. As in the case of free- 
flight stability, a theory of oscillatory stability has ex 
isted for many years. Also similarly, this theory has 
been of little, if any, help in understanding and cor 
recting the faults in seaplanes. In nearly all cases these 
faults have been corrected, if at all, by a cut-and-try 
method without any fundamental understanding of 
what was going on. To a partial extent such funda 
mental understanding now exists for certain types of 
hydrodynamic instabilities. Certain facts concerning 
the physical phenomena of instability are presented in 
this paper in the hope that those working on the theory 
of stability may be able to better match theory to prac 


tice. 


REGIONS OF FLow 


The physical flows of water around a seaplane hull 
during take-off or landing cover a wide range of condi- 
When at rest and at low speeds, the hull be- 
haves very like a displacement boat, and virtually all 
In this régime re- 


tions. 


lit is derived from displacement. 
sistance is mostly surface friction and wave drag. 
Aerodynamic forces on the seaplane are normally neg- 
ligible except for propeller thrust. Longitudinal oscilla- 


Received November 18, 1949. 
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‘The hydrostabilizer is a lateral stability sponson attached 
on either side of the hull at the main step, as in Boeing 314 and 
Martin 130 seaplanes. 


tions are usually highly damped under these condi 
tions. At the opposite end of the speed scale near 
flying speed, entirely different conditions exist. In 
this régime displacement effects tend to be negligible, 
and lift is obtained from the planing action of the hull 
and aerodynamic lift on the wing. Resistance stems 
partially from surface friction and partially from a wave 
drag analogous to induced drag. Nearly all forms of 
hydrodynamic instability occur in this régime of planing 
flow. In the transition region between the displace 
ment régime and the planing régime mixed flows exist. 
Apparently, all cases of instability existing in this ré 
gime are identical to, or closely associated with, those 


in the planing régime. 


STABILITY LIMITS 


It has been found that the angle of trim of the sea 
plane hull to the free water surface is of primary im 
portance in describing stability limits. When operat 
ing at a given speed and weight, it will normally be 
found that within certain trim limits all longitudinal 
motions are damped and that within certain other limits 
longitudinal oscillations continue at finite magnitudes. 
The angle limits that divide these regions are sharp and 
definable in most experimental cases. In general, it 
has been found that the only effect of center of gravity 
is to change the trim angles that are obtainable because 
Center of gravity appears to 
For 


of control limitations. 
have no effect upon the angle limits for stability. 
this reason, it is convenient to use a trim angle versus 
speed plot to describe stability limits. Such a plot is 
shown in Fig. 1. On this plot are shown not only the 
trim limits for stability but also the free-to-trim angles 
for several elevator positions at one center-of-gravity 


location. 
TYPES OF LONGITUDINAL INSTABILITY 


The first type of longitudinal instability to attract 
much attention was a low-angle type of instability. In 
this type of instability, if the trim is reduced below cer- 
tain values, a combination pitching and heaving oscilla- 
tion starts. The motion is approximately sinusoidal 
with the two motions out of phase with each other and 
with a period of several seconds. This motion is com- 
monly described as “‘porpoising.”” The magnitude of 
the motion may be small or large depending upon the 
extent to which the seaplane is trimmed into the un- 
stable region of trims and upon the configuration of the 


hull. At low speeds the stability limit occurs at a high 
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Fic. 1. Trim limits for stability of a typical flying boat. 


trim angle, and at high speeds the trim limit is small. 
The stability limits for this type of instability for a con- 
ventional seaplane are shown in Fig. 1. In this type of 
motion the hull will practically never leave the water, 
and the motion will stop shortly after the trim has 
been shifted into the stable region. 

A high-angle type of instability, commonly described 
as “‘skipping,’’ has come into prominence in the last 
few years. In this type of instability the seaplane will 
commonly attempt to land at a fairly high angle, but 
immediately following contact with the water the angle 
of attack will increase and the craft becomes air-borne 
again. This cycle is repeated until the seaplane has 
lost flying speed or until flight is resumed. The action 
appears similar to the skipping of a stone on water. 
The same phenomenon may be experienced during take- 
off if the trim angle is increased into the skipping range. 
In this case the angle of attack increases suddenly, and 
the seaplane leaves the water. This maneuver is usually 
followed by a short pull-up, and the craft seldom con- 
tacts the water again. Skipping invariably involves 
leaving the water entirely, and, although the motion is 
repeated several times in some landings, it cannot be 
analyzed by any linear theory nor is the motion sinu- 
soidal. 

Trim limits for skipping are hard to measure full 
scale but are readily determined in constant-speed 
model tests. Two limits normally exist. There is a 
limit above which the trim must be raised to start the 
skipping, and there is a somewhat lower limit below 
which the trim must be reduced to stop the skipping. 
Typical limits are shown in Fig. 1. The severity of the 
motion at angles in the unstable region differs widely 
between hulls and is of more interest than the limits 
themselves. 

There is another type of high-angle instability. This 
motion is similar to low-angle porpoising and has a 
single upper limit as compared with the double upper 
limit of skipping. The hull does not leave the water, 
and the motion appears to be nearly sinusoidal. The 
appearance of this instability to the pilot is essentially 
identical to the low-angle type, except that the angle 
must be reduced to stop the motion. 


HYDRODYNAMIC STABILITY 
A theory to account for hydrodynamic stability Was 
given by Glauert in R. & M. 1493 in 1932. pyj. 
theory is similar in basis to most theories of stabjfity 


for complex systems and involves certain assumptions | 


and limitations that are usually considered admissible 
These assumptions include: that the motion is smaj 
that linear coefficients can be used, that no discontiny 
ties exist, and that terms of order other than the first 
may be neglected. The Glauert theory has heey 
difficult to apply and has seldom given results that agre 
with experiment. So far as is known, the Glauer 
theory has never been used to assist in understanding 
instabilities or correcting them. It is important 
that the Glauert theory be expanded into a useful and 
correct form or that a new theory be developed which 
adequately explains the physical phenomena. Th 
absence of a proper understanding of these stability 
problems has prevented the correction of serious faults 
in many seaplanes. 


INVESTIGATION OF WATER FLOWS 


In order to determine the overall action of a seaplane, 
a detailed investigation of the water flows is necessary 
Nearly all porpoising occurs in the planing region where 
the afterbody and tail extension are free of the water 
and only the straight portions of the main planing sur 
face and the afterbottom need be studied. Nearly all 


hulls consist of Vee bottoms with flat surfaces except | 


for a small amount of chine flare. To a first approxi- 
mation the wetted area in the planing régime will be 
the area of the bottom which would be wetted if the 
heave and trim were held fixed and the velocity reduced 
to zero—that is, the water is deflected but little before 
contacting the planing surface. Using this approxi 
mation, it can readily be seen how the wetted area varies 
with heave or displacement. A typical wetted area 
sketch is shown in Fig. 2. With flat surfaces, the area 
is triangular up to the point where the chine at the step 
is wetted. For greater displacements the areas are 
double trapezoids. The wetted area increases as the 
square of the heave until the chine is wetted and linearly 
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Fic. 2. Wetted area vs. heave. 
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from then on. The lift should be nearly proportional 
to wetted area. Also, the motion of the center of lift 
should bear some relation to the motion of the center 
of wetted area. A typical set of lift data is shown in 
These data are for a Vee planing surface with 
These data 


Fig. 3. 
chine flare as used in a conventional hull. 
show a definite change in the curves at a value of heave 
nearly identical to that for chine wetting when at rest. 
The data also show a parabolic shape followed by a 
linear shape and have been plotted in coefficient form, 
with the coefficients based on dynamic considerations 
rather than the conventional displacement coefficients. 
These data are not greatly affected by Froude Number 
or displacement effects in the regions where insta- 
bilities occur. 

In the case of a hull planing upon its forebottom, if 
the trim is increased sufficiently, the afterbottom will 
contact the water. As the trim is increased, still further 
considerable lift will be carried upon the afterbody. 
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This contacting of the afterbody will produce a marked 
change in the moment characteristics and some change 
in tit lift. The flow in the wake of the main step is 
complicated, and the water surface may be either above 
or below the undisturbed water surface and the angle 
may be inclined to the horizontal. Thus, it is difficult 
to predict either the trim angle at which the afterbody 
will contact or the forces following contact. Since the 
water surface will normally be curved, the forces may 
depend more on the curvature than upon the angle, and 
so-called ‘“‘sticking 


discontinuities may result from 


flow,’ as over a spoon when held with the convex side 


under a water faucet. 


STABILITY THEORY 


An examination of the basic flows as given above is 
necessary before the application of planing surface data 
to the stability theory. The examination above shows 
a number of different flows over the hull. The stability 
theory must be applied separately to each region of 
flow and must not be applied at all in the immediate 
vicinity of the changes between types of flow. The ap- 
plication of the linear theory in the region where the 
the displacement will be 
The next 


lift varies as the square of 
difficult at best and not necessarily possible. 
consideration is to determine in what flow regions in- 
stabilities actually exist on seaplanes. These regions 
are sharply defined for normal hulls and can be plotted 
on the trim versus speed plot. If the aerodynamic lift 
is assumed constant with angle, the hull will operate at a 
given weight and speed at constant lift coefficient. For 
a constant water lift coefficient, the variation of trim 
with displacement can be determined from Fig. 3 and 
The chine is wetted at low angles 
and free at high angles. A typical afterbody has been 
shown in Fig. 4 but was not considered in the lift 
At some trim angle the sternpost will 


plotted as in Fig. 4. 


calculations. 
touch. 

The regions of flow are plotted on a speed versus trim 
plot in Fig. 5. Not only is the lower envelope of these 
curves similar in shape to the lower limit but, for several 
cases, it is identical to the lower limit. 

This is a most important fact. For several cases, it 
shows that stability was not dependent upon some 
fine gradation of damping but was determined by 
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changing physical phenomena. It has frequently been 
observed that wetting of the sternpost is never present 
in lower-limit porpoising, but the chine limit is be- 
lieved to be a new discovery. Model tests of the XP- 
BB1 confirm the coincidence of the lower limit with 
both the force break in the lift curves and the point of 
visual wetting of either the sternpost or chine. In the 
full-scale tests it was not possible to observe the stern- 
post wetting, but in the speed region where the stern- 
post was probably free the lower limit was identical to 
the point at which the chine was wet. If the chine were 
visually wet more than '/, in. forward of the step, a defi- 
nitely perceptible heaving oscillation could be felt. 
This was checked at a number of speeds. 

Additional tests were run on the model at various 
weights and without the wing, and in all cases the lower 
limit was defined as described above. 

Extensive dynamic stability tests were run on a 
model of the Boeing 314 which had hydrostabilizers. 
The full-scale tests had never indicated the existence of 
a lower limit to stability. The model tests showed that 
a mild type of lower-limit porpoising was possible in 
certain limited regions of speed and trim. These are 
shown in Fig. 6. When the hydrostabilizers were re- 
moved from the model, conventional lower limits were 
obtained as shown in Fig. 6. The lower-limit porpoising 
It is obvious that the hydro- 
Here 


was violent in this case. 
stabilizers have a major effect on the lower limit. 
again, the physical flows were investigated, and it was 
found that there was no instability when the hydro- 
stabilizers were contacting the water. Below certain 
trim angles, the hydrostabilizers planed on the water 
with a triangular wetted area. Since the step of the 
hydrostabilizers intersected the hull proper slightly 
above the chine, it was possible to wet the chine without 
wetting the hydrostabilizers in certain regions, and it 
was only in these regions that instability could be found. 
Since the angle limits for this instability were narrow, 
only a small oscillation could exist without exceeding 


the limits and a mild motion resulted. 

It might be presumed that a complete explanation of 
lower-limit porpoising has been found. This is cer- 
tainly not the case. Lower-limit porpoising is essen- 
tially a coupled oscillation of a system with at least two 
degrees of freedom. Without going into the mathe- 
matics of the cases, it would appear that for certain com- 
binations of physical phenomena stability normally 
exists and that for certain other combinations insta- 
bility normally exists. Nothing has been presented 
herewith to prove that an unstable case might not be 
found under physical conditions normally giving sta- 
bility or vice versa. Damping of sufficient magnitude 
and properly applied should make any case stable. 
Possibly, such damping can be supplied by an aero- 
dynamic tail surface. In any case, data conflicting with 
the above findings may exist. 
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DISCUSSION 


In examining data it is important to know the basis 
of the data. 
of oscillation (such as a double amplitude of 3°) as the 
definition of the limit of stability, and such data will not 
agree with the data of this paper wherein the limits are 
Much broader limits wil 


Some investigators use a finite amplitude 


for no perceptible oscillation. 
be obtained if a finite oscillation is used as the basis 0! 
stability limits. 

It is not difficult to find common hull types to whieh 
the physical limits described above cannot be applied 
A hull with a long, pointed, main step without any break 
This 


exactly the form of surface used for the afterbody 


between the step and the chines is such a case. 


most recent flying boats. 
surface will vary widely and for certain conditions m4) 
be stable and for certain other conditions may be ut 
stable. It will probably be difficult to find any physical 


The characteristics of such é 
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na in connection with the flow over such sur- 


phenome on 
faces which will tie in with the stability limits. Other 


hull forms not covered include fluted bottoms, flat 
bottoms, Vee steps, sea sleds, and all forms of hull dif- 
fering much from those investigated. 

Most investigators have attempted to check theory 
against practice in the flow region where the hull char- 
acteristics most nearly approximate those of a flat plan 
ing surface. This is, of necessity, the region where the 
afterbody is free and the chines are wet. It is not sur 
prising that instability was indicated for most cases. 
The fact that the hull may operate under widely dif- 
ferent conditions has usually been neglected. It is 
these different conditions that supply the stable regions 
for most seaplanes. 

For conventional hulls, it would appear that the sea- 
plane designer has no control over the lower stability 
limits. The sternpost angle is the only important vari- 
able. The designer has considerable latitude in the ap- 
plication of the hull in order to provide that the seaplane 
will normally be operated within its stable region. 
This control is exerted mostly by properly locating the 
center of gravity, but it is also influenced to some extent 
by the overall configuration of the seaplane and to some 
extent by its loading. Highly loaded hulls are much 
more susceptible to trouble than lightly loaded hulls. 
It has been found that, if the speed scale is divided by 
the square root of the hydrodynamic lift, the trim limits 
are nearly independent of the weight. Thus, increasing 
the weight moves the trim limits out to higher speeds. 
In some cases this can be serious. Since take-off and 
landing are accelerated conditions, the time spent in any 
given region will be short. If the instability is not great, 
it is perfectly feasible to accelerate through an unstable 
region without ever encountering an oscillation. How- 
ever, at high weights the accelerations are reduced, and 
sufficient time may be spent in a condition to realize 
fully the instability characteristics of the region. Un- 
stable regions should be avoided if at all possible. 

Skipping has not yielded to any simple physical 
reasoning. The hull-water combination is analogous to 
the landing gear-tire-oleo combination of a landplane. 
If a seaplane lands with a vertical velocity, the hull dis- 
placement will exceed the statically required displace- 
ment, and an upward velocity will be imparted to the 
hull. This is similar to a spring effect. Of course, there 
will be some energy absorption as by an oleo, but the 
hull may leave the water and make a second contact. 
This can all occur without the trim angle changing ap- 
preciably. The result is a simple bouncing. If the total 
water reaction results in an increase in trim angle, both 
the hydrodynamic and aerodynamic lift will increase, 
and the vertical velocity leaving the water can be 
higher than that when entering. Such a condition is 
similar to that occurring on an airplane with the main 
wheels ahead of the center of gravity, except that it may 

On the Boeing 314 in its original con- 
figuration, a high-angle landing would put the after- 


be more severe. 
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body and a portion of the main bottom in the water. 
Although all of this area was behind the center of grav- 
ity, the resultant forces on the seaplane resulted in an 
increase in trim, and the ship usually left the water im 
mediately. A number of structural failures occurred on 
this and other seaplanes in which the afterbottom struc 
ture failed downward. Apparently, the increase in trim 
results from downward forces upon at least a portion 
of the afterbottom. On the Boeing 314, it was found, 
by the cut-and-try method, that after deepening the 
step this tendency to increase trim was sufficiently re- 
duced to prevent leaving the water on most landings. 
Some semblance of the phenomena remained in the 
scale model after changing the step, and the angles at 
which it occurred were essentially unchanged but the 
magnitude of the motion was greatly reduced. 

Sudden trim increases during take-off resulting in un- 
intentional take-offs also occurred in the Boeing 314. 
These were analogous to the skipping on landing but 
proved that the difficulty was not primarily associated 
with the vertical velocity on landing and was more a 
function of the hydrodynamic moments than anything 
else. Apparently, the moments are discontinuous with 
trim angle. This may not be a static moment condition 
but does appear to be related to the effect of aerody- 
namic ventilation upon the hydrodynamic flow over the 
afterbottom. Deep steps and ventilation nearly always 
reduce the severity of skipping but have little effect 
upon the trim limits for stability. In the case of the 
Boeing 314, the step depth was increased by moving the 
step aft, and all skipping was eliminated in the full- 
scale case both during take-off and landing. Such a 
change would not have been possible in a normal sea- 
plane, but in this particular case there was no lower- 
limit porpoising, and the trim change resulting from 
moving the step had unimportant effects. Without 
the hydrostabilizers this seaplane would have trimmed 
in the lower-limit porpoising range both before and 
after the step change. In the case of the XPBBI, the 
hull was initially designed with an extremely deep step, 
and no evidence of skipping was ever seen in the full- 
Double upper limits of skipping could be 
The de- 


celeration during a landing is sufficiently great to pre- 


scale article. 
found in the model, but the motion was mild. 


vent mild motions from becoming evident to the pilot. 
So far as is known to the author, the XPBBI never left 
the water after its first contact in landing. In landings 
with high vertical velocities, the ship would initially 
sink noticeably below the static displacement and would 
then rise nearly clear of the water, followed by a normal 
slowing down in the water. 

In the XPBBI1 definite evidence was found at full 
scale of an upper-limit instability that closely resembled 
lower-limit porpoising and which had none of the char- 
acteristics of skipping. The motion was mild and had a 
single upper limit as distinguished from the double 
upper limit of skipping. There was no tendency to 
leave the water, and the condition could be found from 
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hump speed to near take-off speed. There is no assur- 
ance that this instability was encountered on the model, 
possibly because the model did not aerodynamically re- 
semble the full-scale article at the extremely high angles 
at which this motion was encountered. During this 
motion, the afterbottom was wet to a great extent, and 
the main step was wet with a triangular area. It oc- 
curs to the author that this instability may result from 
the same combination of events applied to the after- 
body which causes lower-limit porpoising when applied 
to the forebody alone. 


AERONAUTICAL 


SCIENCES—AUGUST, 1950 
end that the longitudinal hydrodynamic Stability prob. 
lems will become well understood. It appears necesgyp 
to make the theory agree with the physically observabj; 
facts if success is to be had. In the case of lower-linit 
porpoising, there appears to be excellent agreement jy 
tween the regions of stabilities and the regions of cert, 
basic flow types. Skipping and a new type of Upper 
limit porpoising have been described. 

Since preparation of this paper, a paper “‘Some Ax 


pects of the Flow Round Planing Seaplane Hull ¢; | 


Floats and Improvement in Step and Afterbody p, 





sign,’ by K. M. Tomaszewski and A. G. Smith yas 
given at the VII International Congress of Applied 


Mechanics. The latter contains information concerning 


CONCLUSION 
It is the hope of the author that the material of this upper-limit porpoising which complements the informa 


paper will stimulate further investigations toward the — tion on lower-limit porpoising given above. 
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A Mechanical Analogy for Hypersonic Flow 


H. REESE IVEY* 
Langley Aeronautical Laboratory, N.A.C.A. 


SUMMAR\ 


\ mechanical analogy for hypersonic flow is described which 
wy be useful for the Mach Number range of 5 to infinity 
The basic characteristics of the various flow types are reviewed: 


ll 


gas dynamics, slip flow, and free molecule flow. From that dis 
cussion the requirements are determined for an analogy to simu 
late hypersonic flows. The theory of the analogy is developed, 
and a simple apparatus is described for demonstrating it. The 
functions of the various parts of the device are explained, and 
typical photographs of the resulting flow are presented. The 
pictures show the gradual degeneration of shock waves as the 
effective operating altitude is increased from the gas dynamics 
region to the free molecule region. Possible uses of the analogy 
are suggested, and the operating limitations are discussed. 


INTRODUCTION 


Ss“ OF THE FEATURES that characterize modern 
aerodynamics may be conveniently illustrated by 
considering missiles. For example, missiles designed 
for operating at high supersonic speeds present several 
complex aerodynamic problems related to the large 
changes in atmospheric density encountered along the 
flight paths. At the lower altitudes the mean free paths 
of the molecules are short, and the air can be considered 
as a continuous medium. In this case the laws of gas 
dynamics apply: strong, sharply defined shock waves 
occur and fully developed boundary layers are formed. 

As the missile rises into the upper atmosphere, the 
mean free paths of the molecules become appreciable 
but are still small by comparison with the characteristic 
missile dimension. Since shock waves are at least one 
mean-free-path length thick, the shock waves on the 
missile are no longer sharply defined but are noticeably 
thick and diffuse. Similarly, the boundary layers begin 
to behave differently from the boundary layers of con- 
tinuous fluids. For example, a molecule can carry 
momentum from a relatively high speed level of the 
boundary layer to the surface of the missile without al- 
ways hitting intermediate molecules or losing momen- 
tum to them. Thus, the flow velocity at the missile 
surface is no longer zero, as is the case for continuous 
fluids, but is higher, being determined from statistical 
considerations of the velocities of the molecules hitting 
the surface. Since the velocity at the surface is not zero, 
the flow can be said to “‘slip,”’ and the missile is operat- 
ing in the ‘‘slip-flow”’ region of flight. At a still higher 
altitude, the molecular mean free path becomes large 
with respect to the missile dimensions. There, the 

Presented at the Aerodynamics Session, Eighteenth Annual 
Meeting, I.A.S., New York, January 23-26, 1950. 

* Aeronautical Research Scientist, Theoretical Aerodynamics 
Division 


molecules in the free stream can hit the missile surface at 
full speed and are re-emitted without direct inter- 
ference from other molecules. Shock waves and 
boundary layers are no longer recognizable as such in 
this free molecule flow. 

Relatively little information is available concerning 
the transition of the flow from one type to another as 
the flow density is varied while the Mach Number is 
kept constant at some high supersonic value. There- 
fore, it seems desirable to develop an apparatus that 
can provide visual information that may give a better 
understanding of the phenomena and can give qualita 
tive force measurements. The preceding discussion 
suggests that the flow in such an apparatus must be- 
have as though it consists of finite particles or molecules 
rather than acting as a continuous medium. The 
present paper presents a mechanical analogy for hyper- 
sonic flows which uses hard elastic pellets to represent 


the molecules of a real gas. 


BASIS FOR THE MECHANICAL ANALOGY 


The Mach Number .1/ of a flow is usually defined as 
the ratio of the directed macroscopic speed I’ of the flow 


to the speed of sound a 
M = V/a (1) 


From the kinetic theory of gases, with the usual 
limitations of equipartition of energy and _ perfect 
elasticity, the speed of sound is determined by the 
number of excited degrees of freedom mn of the gas 
molecules and the root-mean-square random speed 


Vrms OL the molecules—that is, 
a = [(n + 2)/3n] O,ms (2) 
Combining Eqs. (1) and (2) yields 
M = (V/d;ms)[3n/(n + 2)] (3) 


Eq. (3) serves as the basis for a hypersonic flow analogy. 
For a given type of ideal molecule, the Mach Number 
of the flow is determined by the ratio of the free-stream 
speed to the root-mean-square random speed of the 
molecules. 

In addition to securing the proper Mach Number, one 
must obtain the desired mean free path of the molecules 
in order to have the right type of flow. The mean free 
path L is a function of the effective diameter d and the 
number of molecules V in a unit volume 


L 


0.2251/ Nd? (4) 
Eq. (4) shows that a few large molecules can give the 
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same mean free path as many small molecules. As long 
as the mean free path is large compared with the par- 
ticle size, the molecules of a real gas can be represented 
by relatively few small elastic spheres in order to form a 
mechanical analogy of the real flow. The effective 
Mach Number [Eq. (3)] is determined by the ratio of 
the directed speed of a stream of pellets to their random 
speed. 

The Reynolds Number of the flow is proportional to 
the ratio of the Mach Number and mean free path for 


given dimensions and gas properties. 


A SIMPLE MECHANICAL ANALOGY APPARATUS 


In a simple piece of equipment for demonstrating 
the mechanical analogy, the elastic pellets are simply 
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dropped onto deflectors and then onto a model to be 
tested (Fig. 1). 

The various parts of the apparatus and their function 
are as follows: 

(A) Pellet storage box stores the pellets and supplies 
them in a steady stream. 

(B) The flow regulation valve distributes and regy 
lates the mass flow. 

(C) The acceleration section gives directed motion ty 
the stream of pellets by allowing them to fall. 


(D) The deflectors give small disturbance velocities ty § 


the pellets. 

(E) The smoothing section furnishes the pellets with 
a space in which the disturbed pellets can collide with 
each other and develop a random disturbance velocity 
superimposed on the general stream velocity. 

(F) The model is placed at the lower end of th 
smoothing section so that the pellets will stream over 
and around it. 

Gravity can be used to furnish the required partick 
speeds, but in that case the length of the acceleratioy 
section must be long when compared with the length oj 
the model in order to approximate unaccelerated flow, 

The flow density at constant Mach Number may be 
At high 
flow densities, the flow corresponds to the gas dy- 
namics type of flow, and regular shock waves are 


varied by means of a flow-regulation valve. 


formed as evidenced by a sharp demarcation between 
high-density and low-density regions of the flow. At 
lower flow densities slip flow occurs, and at sufficiently 
low densities the flow corresponds to free molecule 
flow. With free molecule flow, the pellets that hit the 
model bounce off and leave the vicinity of the model 
without running into other pellets. In this case no 
phenomenon that can be identified as a shock wave 
occurs, and only a gradual increase in density exists 
near the front portion of the model. Weak disturbances 
from the model are visible considerably ahead of the 
position a shock wave would have if the flow were 
sufficiently dense. 

One parameter that is of some importance in simulat- 
ing a supersonic flow is the ratio of the specific heats 0! 
the gas y. If the energy is equally divided among the 
various degrees of the then 
kinetic theory gives the ratio of specific heats simply asa 
function of the number of degrees of freedom 


(n + 2)/n (5 


freedom of molecules, 


7 = 

Air molecules at room temperature possess many 0 
the features of rigid dumbbell molecules, which have 
three translational degrees of freedom and two rota: 
tional degrees. Then, for air, 


y= 7/9 = 14 


In order for a hypersonic analogy to give the prope! 
shock-wave angles, the effective value of y should be 
close to 1.4. If spherical pellets are used in the analogy 
to represent molecules, the value of y analogous to the 
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MECHANICAL ANALOGY 


ratio of specific heats in a real gas is found from Eq. (5). 
A spherical pellet has three translational and three rota- 
tional degrees of freedom; however, the rotational de- 
grees may not be fully excited. Thus, one might expect 
a flow of spherical pellets to behave somewhat like a 
gas with fewer than six degrees of freedom. The values 
of yfora real gas having three and six degrees of molecu- 


lar freedom, respectively, are 


n+ 2 5 = a) = 
“<* =? = 1.67; % =~ = 1.38 
‘ n 3 18) 


From these considerations, a flow of spherical pellets 
should behave somewhat like a gas having a value of y 
between 1.33 and 1.67. In a real gas having nearly 
spherical effective boundaries for the molecules (for 
example Ne), the excitation of molecular rotations may 
take ten collisions, whereas only one or two are required 
for approaching translational equilibrium. The ratio of 
specific heats is discussed in more detail subsequently in 
the present paper. 

When the lower density flows are being considered, 
the proper type of molecular reflections from the model 
surface must be secured. When an air molecule hits a 
surface, it seems to stick momentarily and then leave 
again in some random direction that is not related to its 


original direction of motion. In order to simulate such 


diffuse reflections with the mechanical analogy, the 
model surface must be machined with a particular type 
of surface roughness designed to make the pellets 
bounce in random directions with a Maxwell distribu- 
tion of speeds. The surface pressures are not critically 
affected by slight deviations from the desired velocity 
distribution. 

The pressure exerted on a surface by a low-density 
airflow depends on the momentum change of the 
molecules hitting the surface. The momentum change 
depends partly on the velocity with which molecules are 
re-emitted from the surface. For a translational ac- 
commodation coefficient of 1, the average velocity of 
emission is proportional to the square root of the sur- 
face temperature. If the surface temperature is less 
than the stagnation temperature (as is the usual case for 
hypersonic missiles at high altitudes), the molecules 
lose part of their energy to the surface. It thus be- 
comes apparent that an inelastic surface in the mechani- 
cal analogy is analogous to a cold surface in air. Simi- 
larly, a perfectly elastic surface in the analogy corre- 
sponds to a surface at stagnation temperature. A vibrat- 
ing rough surface may have to be used in the analogy to 
simulate a surface hotter than stagnation temperature 
by increasing the energy of the particles that hit the 
surface. 

The Mach Number of the flow can be altered by 
changing the angle of the deflectors that are responsible 
for the random motion of the pellets. This change 
essentially corresponds to securing a different speed of 
sound in the real gas. Possibly the simplest and most 
interesting case is that in which no deflectors are in- 


FOR 
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stalled. The speed of sound is then effectively zero in 
the corresponding real gas and the Mach Number in- 


finite. 


TYPICAL VISUAL RESULTS 


Interesting optical studies of the flow can be made 
with ordinary photographic equipment. In order to 
secure information on the instantaneous distribution of 
pellets, short exposure times are used. A series of simi- 
lar pictures taken with a high-speed motion-picture 
camera is sufficient to determine the paths and speeds 
of the pellets. Time exposures are useful in giving 
smooth flow patterns. 

Fig. 2 is a typical picture of the flow around a wedge 
airfoil at a Mach Number estimated from high-speed 
motion pictures of the flow to be about 10. The white 
region represents high-density flow, and the dark region 
represents low-density flow. The flow rate in the 
apparatus is such that the mean free path of the flow is 
small compared with the model dimensions but is not 

This flow rate is analogous to slip flow in a 
Note that the shock wave represented by the 
this 


negligible. 
real gas. 
leading edge of the white flow region is diffuse; 
diffuseness indicates that the mean free path of the 
pellets is appreciable. In fact the mean free path of the 
particles is so large in the dark region of the flow that 
no trailing edge shock wave is noticeable. The blunt 
edge of the test airfoil causes the shock wave to be 
slightly detached from the nose. For this particular 
picture, the ratio of specific heats for the corresponding 
real flow was approximately 1, as is explained later in 
the paper, and, hence, the flow separates from the air- 
foil at the midchord of the lower surface as predicted by 
the author.! 

Fig. 3 shows the same airfoil as Fig. 2, except that the 


flow density is lower. The shock wave is becoming 





flow. 


in slip 


Fic. 2. Hypersonic flow around wedge airfoil 
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more diffuse and the high-density region is incre 
extent. 


ASing jp 


; 





Fig. 4 shows the airfoil in a still lower density floy 
This time the mean free path is so long that the le 
shock wave has broken down and the molecules are re 
flecting from the surface in quasi-specular fashion free 
molecule flow with the ratio of specific heats approxi 
mately 1.33). The influence of the airfoil can be seen 
well ahead of the usual gas dynamics shock-wave log 
tion. 

Fig. 5 shows the flow around a circular cylinder wher 
the mean free path is in the slip-flow region with th 
ratio of specific heats in the corresponding real gas equal 
to approximately 1.0. Again the diffuse shock wave ca 
be seen. This type of picture should serve as a partial 
check on the hypersonic method of reference 1, which 
used the assumption of y = | because of resulting sim 
plifications. 





; é' ae. Fig. 6 pictures the flow around the cylinder of Fig, 5 
Fic: 3. Hypersonic flow around wedge airfoil in low-density , acs Se a eg ; 

slip flow. at a lower density. The density increase in front of th 
cylinder has become gradual, and the density increas 


begins ahead of the gas dynamics shock-wave location 


DETERMINATION OF OPERATING CONDITIONS 


The determination of the operating conditions in an 
analogy apparatus rests primarily on the development 
of instrumentation, but a possible approach to the prob- 
lem can be discussed briefly. 


ading § 





When w is the weight of pellets in the pellet storage 
box 6 and ¢ is the time, the weight flow rate dw/dt in the | 
tunnel can be found from the rate of change of weight oi ' 
pellets (dw/dt) in the pellet storage box and can be | 
written as 

dw/dt = —(dw/dt), 6 


The directed stream velocity |’ in the test section is 





Fic. 4. Hypersonic flow around wedge airfoil in free molecule 
flow. 








Fic. 6. Hypersonic flow around circular cylinder in low-density 
Fic. 5. Hypersonic flow around circular cylinder in slip flow, slip flow. 


found from the acceleration of gravity g, the height j 
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MECHANICAL ANALOGY 


that the pellets fall from the flow regulation valve to the 
flectors, and the distance / from the valve to the test 


de 
section. This velocity must be corrected for the drag 
force D on the deflectors. Air resistance is eliminated 


by operating in a vacuum. The approximate directed 


D )+ a( D ) I/s 
dt 2h\dw/ dt 


l , 
: —-¥V Be < 
h dn 5 
(7) 


speed is 


V = V 2gh 


The density p of the flow is found from the test-section 
cross section A, the velocity, and the mass flow. 


p = (dw/di)/AVg (8) 


The speed of sound a in the analogous real gas is found 
from the ratio of specific heats y, the pressure p exerted 
on the test-section walls by the pellets in the corre- 
sponding real gas, and the density p. Thus, with all the 
assumptions previously made, 


/ 
a= WV yp/p (9) 


Eqs. (1), (7), and (9) determine the Mach Number. 
Eqs. (7) and (8) serve as a basis for determining the 
dynamic pressure. 

The mean free path of the pellets is important in de- 
termining the type of flow over the model: free mole- 
cule flow, slip flow, or gas dynamic flow. Since the mean 
free path is essentially determined by the cross section 
of particles in a given volume, numerous small pellets 
are better than a few large ones if the mean free path is 
In the 
present investigation, solid glass beads approximately 
0.006 in. 
and are readily obtainable. 

The mean free path Z of the pellets is a function of 
the number .V in a unit volume and of their diameter d 


L= Nd? 


to remain larger than the molecular diameter. 


in diameter were found to be satisfactory 


0.2251 


OPERATING LIMITATIONS 


The mechanical analogy apparatus possesses certain 
basic limitations that must be considered. 

It has been pointed out that the spherical pellets 
have six degrees of freedom instead of the five degrees 
that cool air has. The three rotational degrees may not 
be fully excited. Hence, the equivalent ratio of specific 
This 


will give shock-wave angles and 


heats y is between 1.33 and 1.67 instead of 1.4. 
analogous value of ¥ 
density ratios that are slightly in error and pressures 
that are approximately correct except at extremely high 
supersonic speeds. 

Other values of y can be simulated by means of special 
are attainable by using 


designs. Lower values of ¥ 


more complicated ‘‘molecules,’’ and higher values are 
obtainable by eliminating some of the degrees of free 
dom. For example, if the apparatus is only one “mole 


cule” wide, the spheres have only two translational de 
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grees of freedom and one rotational degree (the rotation 
may not be fully excited). The simulated value of y is 
then between 1.67 and 2.0. 

Whereas air molecules essentially are perfectly elas- 
tic, the small pellets used in the mechanical analogy 
possess a definite internal damping that is related to the 
This damping limits the 
For low 


coefficient of restitution. 
number of collisions allowable in the flow. 
flow densities where free molecule flow exists, the damp- 
ing is not important, but, as the flow density is increased 
until gas dynamic flow begins, the damping becomes 
important and sets the limit on the usable flow densi- 
ties. If the flow density is high, the random motions 
will damp out as though the molecule represented by 
the pellet had an infinite number of degrees of freedom 
which shared equally in energy. The value of y simu- 
lated by high-density flow in the mechanical analogy is 
then 


lim 
—> 


7 = [(1 + 2)/n] = 1 

Thus, some control can be exercised over the value of y 
represented by the flow of pellets; however, except for 
the case of y = 1, the flow density must be kept low 
because of internal damping. If the Mach Number is 
known for simulated gas dynamic flow over a wedge, the 
effective value of y can be determined by measuring the 
shock-wave angle and applying the classical shock-wave 
equations. In order to minimize the effective varia- 
tion of y with density, it is necessary to use elastic 
particles (glass beads were used in the present investi- 
gation) and a massive model (no thin skin) made 
either of an elastic material or a material that has a 
high modulus of elasticity. 

The importance of y in determining the pressures and 
shock-wave angle at supersonic speeds can be esti- 
mated from the following table for the flow over a 4 
slope at a Mach Number of 7: 

Ratio of Pressure 
Coefficients 
to That of Air 


Pressure 


1 Shock Angle Coefficient 


1.0 10° 29’ 0.02555 0.9538 
1.4 (air 11 0.02682 1.0 
1.667 Li? 23’ 0.02775 1.035 


At higher spéeds the per cent error is more for 
straight-sided shapes but is small for many curved 
shapes. 

It is probably well to take into consideration the fact 
that solid pellets do not necessarily reflect specularly 
from smooth surfaces. For deflections at small angles 
of incidence the reflection angle may be less than the 
angle of incidence because the coefficient of restitution 
is less than 1. At high angles of incidence, the angle of 
reflection may be larger than the angle of incidence be 
cause of the drag exerted on the pellet by the surface 
friction which causes the pellet to rotate. 

One of the most important problems in developing 


the mechanical analogy is the overcoming of the effects 
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of inelasticity of the pellets. One approach to the 
problem would be to add continually random energy to 
the flow by means of vibrating surfaces or by means of 
electrical or magnetic excitation of some kind. Another 
approach is to increase effectively the efficiency of the 
collisions; for example, possibly the mechanical colli- 
sions can be replaced by some more efficient electrical 
repulsion between pellets. 

USES THE MECHANICAL ANALOGY 


POSSIBLE OF 


The successful application of the mechanical analogy 
requires a certain understanding of the basic behavior 
of air molecules under similar circumstances. For ex- 
ample, in order to secure true pictures or force data for 
an airfoil in free molecule flow, it is necessary to make 
the proper adjustments to the surface finish so that the 
pellets will undergo diffuse reflection. n- effect, the 
mechanical analogy apparatus is a device for carrying 
out involved computations, provided it is given the 
proper basic functions with which to work. 


One of the most important characteristics of the 
mechanical analogy is that the flow pattern, the density 
distribution, and even the paths of the individual 
“molecules” are readily apparent to the eye or to a 
camera. 
attaining a quick understanding of the mechanism of 


This flow visibility is vitally important in 


various phenomena. 
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The mechanical analogy may be useful in sever 
fields of aerodynamics as a means of supplementin, 
other equipment, particularly in preliminary wo, 
where the ease of changing the model, the Mach Nyy, 
ber, and the flow density may be important. Probab) 
the 
either free molecule flow or slip flow, since the intern 


most accurate information will be obtained {, 
damping of the particles is not so large. 


CONCLUDING REMARKS 


A mechanical analogy has been described which js 
suitable for studying hypersonic flows at speeds frop 
Mach Number 5 to infinity. Photographs have beep 
presented which show the gradual flow evolution g 
high altitudes where discontinuous shock waves po 
Possible uses and limitations of the 
The 


limitation of the analogy is the inelasticity of the par. 


longer occur. 
apparatus have been discussed. most serious 
ticles used to simulate air molecules. At the present 
stage of development the analogy can be used for 
visual studies of various flow phenomena and fo 
qualitative force measurements. 
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First U. S. National Congress of Applied Mechanics 


Plans have now been fo-mulated to hold the First U.S. National Congress of Applied Mechanics on 
The Illinois Institute of Technology will bz host to the Congress. 
Members of the Institute, one of the sponsoring societies, are cordially invited to submit papers, 
which should constitute original research in applied mechanics. 
words or the equivalent in equations, tables, and diagrams.) Abstracts must be submitted to the Chair- 
man of the Editorial Committee before Ap-il 14, 1951; complete papzrs must be submitted prior to June 


Additional information about the Congress may be obtained by writing to N. M. Newmark, 
Secretary of the Congress, University of Illinois, Urbana, III. 


(Their length must not exceed 5,000 


All members are invited to attend. 
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RIEF REPORTS of investigations in the aeronautical sciences 

and discussions of papers published in the JOURNAL are 

presented im this special department. The publication will be 

completed approximately 8 weeks after receipt of the material. The 

Editorial Committee does not hold itself responsible for the opinions 

expressed by the corre spondents. Contributions should not exceed 
800 words in length. 


On the Stability for a Certain Degenerate Type 
of Disturbance of Viscous Fluid Flow 
Between Parallel Walls 


S. |. Pai : 
Institute for Fluid Dynamics and Applied Mathematics, University of 
Maryland, College Park, Md. 


May 23, 1950 


N THE STUDY OF HYDRODYNAMICAL STABILITY of two-dimen- 
I sional parallel flows, H. B. Squire made a great contribution 
by showing that three dimensional wavy disturbances are more 
' It simplifies the mathe 


However, the analysis of a 


stable than two-dimensional ones 
matical analysis to a large extent. 
certain degenerate type of disturbance treated by Squire in refer- 
ence 1 is incomplete. This note gives the complete analysis of 
the stability for this degenerate type of disturbance. 

For convenience of comparison, we use the same notations as 
those used in reference 1 

We consider the steady two-dimensional flow of an incom 
pressible viscous fluid upon which a disturbance is superimposed. 
The origin O is taken at a point midway between the walls, the 
axis of x is in the direction of the flow and parallel to the walls, 
the axis of y is perpendicular to the walls, and the axis of z is 
parallel to the walls and perpendicular to the other axes. If 26 is 
the distance between the walls, these will coincide with the plane 


y = +b. The main flow is defined by 


u= U(y),, v=0, w=0 (1) 


where u, v, w are the velocities in the direction Ox, Oy, Oz, re 
spectively 

If now a small disturbance be superposed on the main flow, the 
velocities will be given by 


») 
) 


u= Uy +m, ven, w= w (2 


where 1, 2, %,, are small quantities whose squares may be 
neglected 

It will be assumed that any disturbance can be broken up into 
a set of disturbances which is periodic with respect to x and to z 


Thus we may take 


uy = u(y) exp i(jx + kz — ot) / 
%, = v(y) exp 1(jx + ks — at) > (3) 
w, = wy) exp 7(jx + kz — ot) \ 


where j and & are real and positive and o is in general complex. 
For the degenerate type of disturbance v(y) = 0, the equation 
of motion for the disturbance in nondimensional form will be 


A(u) = [D2 — (a? + B?) — tRa(V — c) “u = Q (4) 
with the boundary conditions of uw. 
u=Q0 for y= +1 (5) 


where 


D =d/dy 
a =jb 
= kb 
V = U/Uo 
cUy = o/j 
y =y/b 
R = Uob/v = Reynolds Number 
Uo = some reference velocity 
Squire has shown that if ¢ is real, no real value of Ra can be 


found which satisfies Eqs. (4) and (5). This only means that, 
for finite Reynolds Number R, there is no neutral stable disturbance 
of this degenerate type. 

In order to determine whether the flow is stable or not for this 
degenerate type of disturbance, we have to find a complex ¢ that 
satisfies Eqs. (4) and (5) with positive real a, 8, and R. Reynolds 
Number must be real and positive in order to have some physical 
The imaginary part of this complex ¢ determines 


significance. 
the stability of this type of disturbance. 
Write 
C = CRT IC; (6) 
Eq. (4) becomes 
A(u) = [D? — (a? + B?) — tRal(V — cr — ici)|u = 0 (7) 


Let @ be the function conjugate to u which satisfies the equation 


A(ui) = [D? — (a? + B*?) + iRal(V — cr + ici)]i = 0 (8) 
and the same boundary conditions as u. 

Multiply Eq. (7) by # and Eq. (8) by u and subtract and in- 
tegrate across the field. Then 


- bd | 
0 = f [a Au — uAu| dy = —2iRa | (V — cr)uti dy 


(9) 


performing the same operations and adding, we have, after in- 
tegrating by parts and making use of the boundary conditions, 


eo} na | ao | 
—Rati { , ee dy = 2f _ un’ dy + (a? + B?) { _ uit dy 


(10) 

where u’ = du/dy, and &’ = du/dy 
All the integrals in Eq. (10) are positive real numbers. By 
Hence, for positive 
Hence, all pro 


definition, a and 8 are positive real numbers. 
real Reynolds Number R, c; must be negative. 
files V(y) are stable for disturbances of this degenerate type 

Eq. (9) determines the velocity of propagation of this degen- 
erate type of disturbance. 


REFERENCE 
Stability for Three-Dimensional Disturbance of 


1 Squire, H: B., On the 
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Base Pressures at Supersonic Velocities 


A. Gabeaud 
Centre d’Essais de Méchanique des Fluides, Paris, France 


May 1, 1950 


HAVE JUST NOTICED in the Readers’ Forum (JOURNAL OF THE 
I AERONAUTICAL SCIENCES, Vol. 17, No. 3, P. 185, March, 1950) 
the note by Freeman K. Hill on ‘Base Pressures at Supersonic 
Velocities,” in which the author, after recalling the results ob- 
served in free flight and in a*wind tunnel for Mach Numbers 
up to 2.5 and mentioning the theoretical relationship that | 
proposed some time ago,' drew attention to a recent result ob 
The author suggested a 


tained at VW = 5.1 in a wind tunnel. 


525 
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semiempirical expression for the base pressure for Mach Numbers 
between 2 and 6 as follows: 
9 


, Y 
Cp =f -(1 — 


1 +) (1 
VE —— et Pe (1) 
vo ~ 


where f and \ are two adjustable coefficients whose values are 
determined from the experimental results observed for J = 
2.5and M = 5.1 (f = 0.5394; » = 0.055). 

It has not escaped Hill’s notice that this expression leads to a 
negative value of the base pressure for values of ./ greater than 
9.5, and for this reason he recommended its use only up to J = 6. 

Having recalled this fact, I wish to make the following two 
remarks: 

(1) The theoretical law that I proposed in 1936 gives Cz = 
— 0.034 at M = 5.1; 
in 1949! and which takes the simple form 


the most recent law that I communicated 


2 2yM? — (vy -— 1) 2 3 2 
Cz = _ (2) 
> M? y¥ + 1 ¥+1 M?-8 > M? 


gives Cg = —0.0415. The theoretical result is then confirmed 
very well by experiment (Cg = —0.0457). 

(2) The same theoretical law does not give a perfect vacuum 
until MW = o, 

From these two remarks it would appear that the theoretical 
law, Eq. (2), which is well verified by experiment in the regions 
from JJ = 1.8 to M = 5.1 and which does not predict a perfect 
vacuum until JJ = ©, is preferable to the semiempirical law, 
Eq. (1) [Hill’s Eq. (2)], the use of which presents the incon- 
veniences indicated by its author above J = 5.1. 
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Note on the Calculation of Deflections of 
Indeterminate Structures 


Bruno A. Boley* and Roland H. Moore? 
Goodyear Aircraft Corporation, Akron, Ohio 
May 10, 1950 
FS arpaeeeitenn IS PRESENTED HERE Of the usual method of 
calculating deflections of indeterminate structures by means 

of the principle of virtual displacements. A review of the litera- 
ture, including some of the best-known texts on the subject, was 
carried out; none of the works consulted appeared to mention the 
possibility of this simplification (however well known it may be to 
practicing engineers). It was therefore decided to present it in 
the hope that, once the advisability of mentioning it were pointed 
out, it would find its way in new books and in new editions of old 
ones 

The usual procedure in calculating deflections is as follows. 
Let the deflection 6,4 of a point A in the structure in question be 
required in a given direction. The internal strains in the struc- 
ture caused by the external loads will be denoted by e,, ¢ a 


Y,2 and y,,, following the usual notation. If the structure is in 
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determinate, the evaluation of these strains will require use of vir 
tual work or some equivalent concept related to the distortions of 
the structure. This calculation may be considered the first ste] 
in the determination of deflections, and no simplification of it 1s 
suggested here. 

In the next step, a dummy system of self-equilibrating loads ts 
applied to the structure. ‘he dummy systein is chosen in such 
manner that it does work only when a deflection 6,4 occurs at 
point A in the desired direction and has such a magnitude that 


this work is equal toé,. The internal stresses due to the dummy 


system will be denoted by a;’, ¢,,’,¢,', 7,,', 7,2’, and 7: 

The equilibrium of the structure under the dummy system & 
then tested (though it was established at the outset) by means 0 
the principle of virtual displacements. This principle states that 
a system of forces is in equilibrium if, and only if, the work dom 
by it during any arbitrary virtual displacement is zero. The work 
done by the dummy system must then also be zero when the a 
tual displacements of the structure are used. This work is th 
sum of the external and internal work; or 


= (), from which the deflection 6,4 is easily obtained 


\dx dy d 
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The only requirements that the dummy system is expected to 
satisfy in this derivation are those of equilibrium; at no time is it 
required that the stresses due to it satisfy the compatibility con- 
ditions for the structure. Hence, the following simplification 
may be employed: that, in obtaining the stresses due to the 
dummy load, it is not necessary to ensure continuity of displace- 
ments throughout the structure and that, therefore, the redun- 
dant members of the structure may be neglected in this part of the 
calculations. 

It may seem at first that, if this is done, the displacements used 
to calculate the virtual work violate the conditions of restraint of 
the structure considered. This is not the case, however, since 
the principle of virtual displacements is in reality applied only to 
the simplified structure acted upon by the dummy loads. The 
displacements of the given structure do not violate any of the re- 
straints of the simplified one. 

Three examples of the use of this simplification are indicated 
below 

Let the deflection at point A in the beam of Fig. 1 be required. 
Then the dummy system consists of a 1-lb. load downward at A 
and a '/s-lb. load upward at each end, the redundant end mo- 
ments being neglected. Hence, for this part of the calculations 
the beam reduces to a simple beam. 

Fig. 2b shows the dummy system to be used in finding the 
downward deflection of point A of the truss of Fig. 2a. In cal 
culating the forces due to the dummy system the redundant 
bars may be disregarded; thus, the simplified truss of Fig. 2b 
results. 

In obtaining the shortening of diameter AA’ in Fig. 3a, the 
dummy system of Fig. 3b must be used. The determination of 
the bending moments due to the dummy loads may be simplified 
by the insertion of imaginary hinges (Fig. 2b) which eliminate 
two internal redundants (the third redundant being eliminated 


by considerations of symmetry). 
° 


The Effects of Weight Variation on the Range 
of Air-Borne Craft 


Ralph W. Allen 
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URING A PRELIMINARY ANALYSIS of a long-range aircraft 
design, it is frequently necessary to adjust the final weights 
from the original estimations. The results of this adjustment, 
among other things, usually affects the range of the craft. In 


order to achieve a first approximation as to the magnitude of 
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this effect, without entering into exact calculations, the equa- 
tion and curves derived herein can be applied. 

The converse can also follow—that is, a change in the original 
estimated range of the craft introduces a change in the weight. 
These changes can be visualized immediately from the curves 
shown in Fig. 1. 

The range of an air-borne craft traveling at a constant speed 
and maintaining a constant lift to drag ratio, powered by a system 
that maintains thrust to drag equalization, is expressed by the 
Breguet formula 


R = I(L/D) V loge (We /We) (1) 

where: 

R = range 

I = specific impulse 

L = lift 

D = drag 

We = gross weight of craft 

Wer = weight empty of craft 


The weight of the craft can be broken up into three major 


parts—namely, 
We = Wr + Ws + W; (2) 
in which 
Wr = weight of Fuel 
Ws = structural weight 


Wp = pay load 


The weight empty is defined as the sum of the pay load and 
structural weights 


We = Wp + Ws (3) 


Let We/We be denoted by u; then, if Ry represents the initial 
range based on yo and R; is the new range based on yw, the ratio 
of R, /Ry becomes, from Eq. (1), 


R,/Ro = log wu: /log uo (4) 


This reduces to 


S$ 
+ 
WY 
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From Eqs. (2), (3), and (4), the following relations are de- Under a given load criteria, the structural ratio Ws/We 
usually remains fairly constant with respect to range changes, 
so that for a given structural weight ratio the ratio of Bros 

weight to pay load is determined by Eq. (5). 
From the results of Eq. (5), it is noticed that the gross Weight 
; : — — P of the craft increases exponentially with the range iNcreage 
Wo/Wep = 1/wo~®'/® — (Ws/We)] ‘ Thus, given a gross weight and pay load, a new gross weight jg 

found resulting from the designated range ratio. 
The curves of Fig. 1 quickly aid in determining these Various 
(Wp/We) + (Ws/We) = po ®/o changes. 


duced: 


wm = [Woe/Welh = [We/(We + Ws)h = po Ro 


reducing, 
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